Scheduled Downtime
On Tuesday 24 October 2023 @ 5pm MT the forums will be in read only mode in preparation for the downtime. On Wednesday 25 October 2023 @ 5am MT, this website will be down for maintenance and expected to return online later in the morning.
Normal Operations
The forums are back online with normal operations. If you notice any issues or errors related to the forums, please reach out to help@ucar.edu

MARBL ERROR in running CESM2.1.3 with --compset BHISTcmip6 --res f19_g16

E

Eun-Jo SONG

EJSONG
New Member
Dear all

I have some problem in running CESM 2.1.3 BHISTcmip6 compset withe resolution f19_g16.
cesm.log as followings:

(Task 132, block 1) Message from (lon, lat) ( 52.017, 36.614), which is global (i,j) (83, 288)
(Task 132, block 1) MARBL WARNING (marbl_co2calc_mod:drtsafe): (marbl_co2calc_mod:drtsafe) it = 1
(Task 132, block 1) MARBL WARNING (marbl_co2calc_mod:drtsafe): (marbl_co2calc_mod:drtsafe) x1,f = 0.6309573E+001-0.6290290E+001
(Task 132, block 1) MARBL WARNING (marbl_co2calc_mod:drtsafe): (marbl_co2calc_mod:drtsafe) x2,f = 0.1584893E+002-0.1580036E+002
(Task 132, block 1) Message from (lon, lat) ( 53.141, 36.671), which is global (i,j) (84, 288)
(Task 132, block 1) MARBL WARNING (marbl_co2calc_mod:drtsafe): (marbl_co2calc_mod:drtsafe) it = 1
(Task 132, block 1) MARBL WARNING (marbl_co2calc_mod:drtsafe): (marbl_co2calc_mod:drtsafe) x1,f = 0.6309573E+001-0.6290290E+001
(Task 132, block 1) MARBL WARNING (marbl_co2calc_mod:drtsafe): (marbl_co2calc_mod:drtsafe) x2,f = 0.1584893E+002-0.1580036E+002
(Task 132, block 1) Message from (lon, lat) ( 52.017, 36.614), which is global (i,j) (83, 288)
(Task 132, block 1) MARBL WARNING (marbl_co2calc_mod:drtsafe): (marbl_co2calc_mod:drtsafe) it = 2
(Task 132, block 1) MARBL WARNING (marbl_co2calc_mod:drtsafe): (marbl_co2calc_mod:drtsafe) x1,f = 0.3981072E+001-0.3968936E+001
(Task 132, block 1) MARBL WARNING (marbl_co2calc_mod:drtsafe): (marbl_co2calc_mod:drtsafe) x2,f = 0.2511886E+002-0.2504184E+002
(Task 132, block 1) Message from (lon, lat) ( 53.141, 36.671), which is global (i,j) (84, 288)
(Task 132, block 1) MARBL WARNING (marbl_co2calc_mod:drtsafe): (marbl_co2calc_mod:drtsafe) it = 2
(Task 132, block 1) MARBL WARNING (marbl_co2calc_mod:drtsafe): (marbl_co2calc_mod:drtsafe) x1,f = 0.3981072E+001-0.3968936E+001
(Task 132, block 1) MARBL WARNING (marbl_co2calc_mod:drtsafe): (marbl_co2calc_mod:drtsafe) x2,f = 0.2511886E+002-0.2504184E+002
(Task 132, block 1) Message from (lon, lat) ( 52.017, 36.614), which is global (i,j) (83, 288)
(Task 132, block 1) MARBL WARNING (marbl_co2calc_mod:drtsafe): (marbl_co2calc_mod:drtsafe) it = 3
(Task 132, block 1) MARBL WARNING (marbl_co2calc_mod:drtsafe): (marbl_co2calc_mod:drtsafe) x1,f = 0.1584893E+001-0.1580113E+001
(Task 132, block 1) MARBL WARNING (marbl_co2calc_mod:drtsafe): (marbl_co2calc_mod:drtsafe) x2,f = 0.6309573E+002-0.6290212E+002
(Task 132, block 1) Message from (lon, lat) ( 53.141, 36.671), which is global (i,j) (84, 288)
(Task 132, block 1) MARBL WARNING (marbl_co2calc_mod:drtsafe): (marbl_co2calc_mod:drtsafe) it = 3
(Task 132, block 1) MARBL WARNING (marbl_co2calc_mod:drtsafe): (marbl_co2calc_mod:drtsafe) x1,f = 0.1584893E+001-0.1580113E+001
(Task 132, block 1) MARBL WARNING (marbl_co2calc_mod:drtsafe): (marbl_co2calc_mod:drtsafe) x2,f = 0.6309573E+002-0.6290212E+002
(Task 132, block 1) Message from (lon, lat) ( 52.017, 36.614), which is global (i,j) (83, 288)
(Task 132, block 1) MARBL WARNING (marbl_co2calc_mod:drtsafe): (marbl_co2calc_mod:drtsafe) it = 4
(Task 132, block 1) MARBL WARNING (marbl_co2calc_mod:drtsafe): (marbl_co2calc_mod:drtsafe) x1,f = 0.2511886E+000-0.2504969E+000
(Task 132, block 1) MARBL WARNING (marbl_co2calc_mod:drtsafe): (marbl_co2calc_mod:drtsafe) x2,f = 0.3981072E+003-0.3968851E+003
(Task 132, block 1) MARBL ERROR (marbl_co2calc_mod:drtsafe): bounding bracket for pH solution not found
(Task 132, block 1) MARBL ERROR (marbl_co2calc_mod:drtsafe): (marbl_co2calc_mod:drtsafe) dic = -0.1000000E+001
(Task 132, block 1) MARBL ERROR (marbl_co2calc_mod:drtsafe): (marbl_co2calc_mod:drtsafe) ta = -0.1000000E+001
(Task 132, block 1) MARBL ERROR (marbl_co2calc_mod:drtsafe): (marbl_co2calc_mod:drtsafe) pt = -0.1000000E+001
(Task 132, block 1) MARBL ERROR (marbl_co2calc_mod:drtsafe): (marbl_co2calc_mod:drtsafe) sit = -0.1000000E+001
(Task 132, block 1) MARBL ERROR (marbl_co2calc_mod:drtsafe): (marbl_co2calc_mod:drtsafe) temp = -0.1000000E+001
(Task 132, block 1) MARBL ERROR (marbl_co2calc_mod:drtsafe): (marbl_co2calc_mod:drtsafe) salt = -0.1000000E+004
(Task 132, block 1) Message from (lon, lat) ( 53.141, 36.671), which is global (i,j) (84, 288)
(Task 132, block 1) MARBL WARNING (marbl_co2calc_mod:drtsafe): (marbl_co2calc_mod:drtsafe) it = 4
(Task 132, block 1) MARBL WARNING (marbl_co2calc_mod:drtsafe): (marbl_co2calc_mod:drtsafe) x1,f = 0.2511886E+000-0.2504969E+000
(Task 132, block 1) MARBL WARNING (marbl_co2calc_mod:drtsafe): (marbl_co2calc_mod:drtsafe) x2,f = 0.3981072E+003-0.3968851E+003
(Task 132, block 1) MARBL ERROR (marbl_co2calc_mod:drtsafe): bounding bracket for pH solution not found
(Task 132, block 1) MARBL ERROR (marbl_co2calc_mod:drtsafe): (marbl_co2calc_mod:drtsafe) dic = -0.1000000E+001
(Task 132, block 1) MARBL ERROR (marbl_co2calc_mod:drtsafe): (marbl_co2calc_mod:drtsafe) ta = -0.1000000E+001
(Task 132, block 1) MARBL ERROR (marbl_co2calc_mod:drtsafe): (marbl_co2calc_mod:drtsafe) pt = -0.1000000E+001
(Task 132, block 1) MARBL ERROR (marbl_co2calc_mod:drtsafe): (marbl_co2calc_mod:drtsafe) sit = -0.1000000E+001
(Task 132, block 1) MARBL ERROR (marbl_co2calc_mod:drtsafe): (marbl_co2calc_mod:drtsafe) temp = -0.1000000E+001
(Task 132, block 1) MARBL ERROR (marbl_co2calc_mod:drtsafe): (marbl_co2calc_mod:drtsafe) salt = -0.1000000E+004
(Task 132, block 1) MARBL ERROR (marbl_co2calc_mod:comp_htotal): Error reported from drtsafe
(Task 132, block 1) MARBL ERROR (marbl_co2calc_mod:marbl_co2calc_surface): Error reported from comp_htotal()
(Task 132, block 1) MARBL ERROR (marbl_surface_flux_mod:marbl_surface_flux_compute): Error reported from marbl_co2calc_surface() with flux_co2
(Task 132, block 1) MARBL ERROR (marbl_interface:surface_flux_compute): Error reported from marbl_surface_flux_compute()
(Task 132, block 1) MARBL ERROR (ecosys_driver:ecosys_driver_set_sflux): Error reported from marbl_instances(1)%surface_flux_compute()
ERROR reported from MARBL library
------------------------------------------------------------------------

POP aborting...
Stopping in ecosys_driver:print_marbl_log

;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;

I found that reported grid is near Caspian Sea,
but I do not know how to fix.

Any help and comments are welcome!

Thanks!

EJSONG
 
M

mlevy

Michael Levy
CSEG and Liaisons
Staff member
This is often not a problem in MARBL itself, but rather the first part of the code to trigger an error when your run goes off the rails. Have you made any changes from the standard BHISTcmip6 compset? How long did your case run for before this error cropped up? You can see from the

Code: 
(Task 132, block 1) MARBL ERROR (marbl_co2calc_mod:drtsafe): (marbl_co2calc_mod:drtsafe) dic = -0.1000000E+001
(Task 132, block 1) MARBL ERROR (marbl_co2calc_mod:drtsafe): (marbl_co2calc_mod:drtsafe) ta = -0.1000000E+001
(Task 132, block 1) MARBL ERROR (marbl_co2calc_mod:drtsafe): (marbl_co2calc_mod:drtsafe) pt = -0.1000000E+001
(Task 132, block 1) MARBL ERROR (marbl_co2calc_mod:drtsafe): (marbl_co2calc_mod:drtsafe) sit = -0.1000000E+001
(Task 132, block 1) MARBL ERROR (marbl_co2calc_mod:drtsafe): (marbl_co2calc_mod:drtsafe) temp = -0.1000000E+001
(Task 132, block 1) MARBL ERROR (marbl_co2calc_mod:drtsafe): (marbl_co2calc_mod:drtsafe) salt = -0.1000000E+004

block that something is very wrong -- with the exception of temperature, it is not physical for any of these values to be negative (and temperature probably shouldn't be -1 celsius at this location).

Knowing nothing about your run, I would suggest at the very least switching to the f19_g17 resolution; this is a slightly updated land mask that removes the Caspian Sea from the ocean domain (instead, the land model treats it as a lake). It may just be the case that you are trying to initialize some aspect of your run with data on the gx1v7 grid rather than the gx1v6 grid, and weird things are happening in the Caspian because some data is masked out / missing.
 
E

Eun-Jo SONG

EJSONG
New Member
Thank you for your reply.

I do not change standard BHISTcmip6 compset.
The error immediately occur after job submitting.

Sincerely

EJSONG
 
M

mlevy

Michael Levy
CSEG and Liaisons
Staff member
Have you tried using the f19_g17 resolution instead of f19_g16? I would recommend making a new case to test this, rather than modifying your existing case.
 
E

Eun-Jo SONG

EJSONG
New Member
mlevy said:
Have you tried using the f19_g17 resolution instead of f19_g16? I would recommend making a new case to test this, rather than modifying your existing case.
Click to expand...
Thank you for your comments.
I try to run cesm2 with both resolution f19_g17 and f09_g17 in BHISTcmip6 compset,
but this problem still remains.

Sincerely
EJSONG
 
You must log in or register to reply here.
Top