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

Starting w/ SOM

B

brunke

New Member
I've been trying to run CAM with the SOM, and I can't get it to run. I've followed the directions to add the needed fields to the input dataset and am using that dataset. I've tried running it initially and as a branch run from a run using climatological SSTs as previous posts seem to indicate that I can. Both has not worked.

For the initial run, I get several QNEG4 warnings about the max possible LH flux being exceeded:

QNEG4 WARNING from TPHYSAC , lchnk = 3; Max possible LH flx exceeded at 2 points. Worst excess = -3.2043E-06 at i = 114

And I'm setting reset_csim_iceprops to .true. in my namelist. Any thoughts on what else I need to do?
 
S

suvarchal_kumar@gmail_com

Member
QNEG4 warnings are not unusual so dont worry about it(you get that if you get -ve moisture values check subroutine qneg4.f in physics dir) but does the model stop giving some error when you start initially? what error do you get?

brunke said:
I've been trying to run CAM with the SOM, and I can't get it to run. I've followed the directions to add the needed fields to the input dataset and am using that dataset. I've tried running it initially and as a branch run from a run using climatological SSTs as previous posts seem to indicate that I can. Both has not worked.

For the initial run, I get several QNEG4 warnings about the max possible LH flux being exceeded:

QNEG4 WARNING from TPHYSAC , lchnk = 3; Max possible LH flx exceeded at 2 points. Worst excess = -3.2043E-06 at i = 114

And I'm setting reset_csim_iceprops to .true. in my namelist. Any thoughts on what else I need to do?
Click to expand...
 
B

brunke

New Member
I know, I've noticed them all the time before. But this time, they're paired with the model stopping on the FIRST time step. There's also an indication given that all the ice is sublimated away.

suvarchal said:
QNEG4 warnings are not unusual so dont worry about it(you get that if you get -ve moisture values check subroutine qneg4.f in physics dir) but does the model stop giving some error when you start initially? what error do you get?
Click to expand...
 
S

suvarchal_kumar@gmail_com

Member
i am not sure but check for proper values of icecov in your boundary sst dataset.
can you show log file at point of error.
brunke said:
I know, I've noticed them all the time before. But this time, they're paired with the model stopping on the FIRST time step. There's also an indication given that all the ice is sublimated away.
Click to expand...
 
B

brunke

New Member
Here's what's reported in the log file:

SOMINI: Read som data for dates 816 43200 and 916
0
INITEXT: ocnfrac= 0.000000000000000E+000
INITEXT: ocnfrac= 0.000000000000000E+000
INITEXT: ocnfrac= 0.000000000000000E+000
INITEXT: ocnfrac= 0.872788838702892
INITEXT: ocnfrac= 8.582871866770614E-002
INITEXT: ocnfrac= 4.800755222529873E-002
INITEXT: ocnfrac= 6.636520427649995E-002
INITEXT: ocnfrac= 5.828821236287018E-002
INITEXT: ocnfrac= 0.000000000000000E+000
INITEXT: ocnfrac= 0.000000000000000E+000
INITEXT: ocnfrac= 1.00000000000000
INITEXT: ocnfrac= 1.00000000000000
INITEXT: ocnfrac= 1.00000000000000
INITEXT: ocnfrac= 1.00000000000000
INITEXT: ocnfrac= 1.00000000000000
INITEXT: ocnfrac= 1.00000000000000
INITEXT: ocnfrac= 1.00000000000000
INITEXT: ocnfrac= 1.00000000000000
INITEXT: ocnfrac= 1.00000000000000
INITEXT: ocnfrac= 0.000000000000000E+000
INITEXT: ocnfrac= 0.000000000000000E+000
INITEXT: ocnfrac= 0.136661332801539
INITEXT: ocnfrac= 0.296736121277573
INITEXT: ocnfrac= 7.368663635916983E-002
INITEXT: ocnfrac= 0.278036514714450
INITEXT: ocnfrac= 1.00000000000000
INITEXT: ocnfrac= 1.00000000000000
INITEXT: ocnfrac= 1.00000000000000
INITEXT: ocnfrac= 1.00000000000000
INITEXT: ocnfrac= 0.746351819132566
INITEXT: ocnfrac= 8.622286955950287E-002
INITEXT: ocnfrac= 1.00000000000000
nstep, te 0 3340167301.58743 0.000000000000000E+000
0.000000000000000E+000 98466.2417923964
icefrac shcoverage= 16.2616 million km^2
icefrac nhcoverage= 5.7440 million km^2
icevol shcoverage= 1.6262 10^13m^3
icevol nhcoverage= 1.1488 10^13m^3
snowvol shcoverage= 0.2420 10^13m^3
snowvol nhcoverage= 0.0031 10^13m^3
QNEG4 WARNING from TPHYSAC , lchnk = 3; Max possible LH flx exceeded at 2 points. Worst excess = -3.6576E-06 at i = 113
-43.7080217391730 1200.00000000000 -11557.7595944038
-2.341661977199045E-004 -1542.11293491346 52449.6260870076
sublimate away all sea ice
something is probably seriously wrong
ice state at dh stop nstep = 0
lon,lat= 115 62
aice 2.367103600791906E-013
nc= 1
ain 2.367103600791906E-013
vin 2.416898670440375E-026
-203.701822772566 1200.00000000000 -3.664691026906504E-006
-1.168408854360669E-005 -1.956657223769444E-006 244442.187327079
hin 1.021036286553664E-013
vsn 1.492887638742491E-028
hsn 6.306811574461930E-016


atm states and fluxes
zlvl = 58.6107045525678
potT = 267.165192884373
Tair = 266.595456659027
Qa = 2.395227224182845E-003
swvdr = 3.266638044879295E-009
swvdf = 5.60537896977331
swidr = 2.719906579292941E-009
swidf = 2.64645762046137
Flw = 268.339688830285
Fsnow = 3.851046112497146E-009

ocn states and fluxes
frzmlt = -28799.7912952528
sst = -1.17204983880527
Tf = -1.80000000000000

srf states and fluxes

ENDRUN:DH
sublimate away all sea ice
something is probably seriously wrong
ice state at dh stop nstep = 0
lon,lat= 106 6
aice 1.606790929024821E-013
nc= 1
ain 1.606790929024821E-013
vin 8.185989502035633E-028
hin 5.094620186214145E-015
vsn 1.350696212187494E-028
hsn 8.406172749601261E-016


atm states and fluxes
zlvl = 55.5999282632879
potT = 253.724164517893
Tair = 253.183091610114
Qa = 5.574551676275658E-004
swvdr = 0.000000000000000E+000
swvdf = 0.000000000000000E+000
swidr = 0.000000000000000E+000
swidf = 0.000000000000000E+000
Flw = 171.264553444071
Fsnow = 5.205608497371813E-019

ocn states and fluxes
frzmlt = -44312.5553359473
sst = -1.35372041769112
Tf = -1.80000000000000

srf states and fluxes

ENDRUN:DH
-31.2052145292305 1200.00000000000 -1.476558076868132E-044
-7.009501694178016E-005 -8.804943979032775E-006 37446.2574350766
sublimate away all sea ice
something is probably seriously wrong
ice state at dh stop nstep = 0
lon,lat= 110 59
aice 2.017814643624049E-013
nc= 1
ain 2.017814643624049E-013
vin 6.167162431873481E-027
hin 3.056357258264858E-014
vsn 0.000000000000000E+000
hsn 0.000000000000000E+000


atm states and fluxes
zlvl = 59.5578481243232
potT = 271.403227672318
Tair = 270.824453735447
Qa = 2.876840852798750E-003
swvdr = 0.000000000000000E+000
swvdf = 0.000000000000000E+000
swidr = 0.000000000000000E+000
swidf = 0.000000000000000E+000
Flw = 207.245659811671
Fsnow = 4.919892299307382E-057

ocn states and fluxes
frzmlt = -128640.306429803
sst = 2.02208176219236
Tf = -1.80000000000000

srf states and fluxes

ENDRUN:DH

Any hints would be very helpful. Thanks.

suvarchal said:
i am not sure but check for proper values of icecov in your boundary sst dataset.
can you show log file at point of error.
Click to expand...
 
You must log in or register to reply here.
Top