[cam-users] CGD Forums
phpbb at cgd.ucar.edu
phpbb at cgd.ucar.edu
Tue Nov 2 09:00:01 MST 2004
Dear cam-users,
As you requested, here is the latest digest of messages posted on CGD Forums forums. Please come and join the discussion!
<< radiation calculation, http://bb.cgd.ucar.edu//viewtopic.php?t=31 >>
Postercathryn.meyer at yale.eduPostedMon Nov 01, 2004 2:26 pm MST, http://bb.cgd.ucar.edu//viewtopic.php?p=119#119
Message: The SWCF and LWCF were nonzero after I replaced cld with zero in the code tphysbc.F90. Would setting cldmin>1.0 in the codes RADCSWMX and RADCLWMX have the effect of not allowing the radiation code to see the clouds? Would this affect the way the model sees the clouds at all?
<br />
<br />Thanks,
<br />Cathy
------------------------------
PosterpjrPostedMon Nov 01, 2004 4:09 pm MST, http://bb.cgd.ucar.edu//viewtopic.php?p=120#120
Message: The SWCF and LWCF were nonzero after I replaced cld with zero in the code tphysbc.F90. Would setting cldmin>1.0 in the codes RADCSWMX and RADCLWMX have the effect of not allowing the radiation code to see the clouds? Would this affect the way the model sees the clouds at all?
<br />
<br />Thanks,
<br />Cathy
<br />
<br />How big were the numbers you saw in SWCF and LWCF?
------------------------------
Postercathryn.meyer at yale.eduPostedMon Nov 01, 2004 5:37 pm MST, http://bb.cgd.ucar.edu//viewtopic.php?p=121#121
Message: Most were zero or on the order of 10^-15, but occasionally there were areas where the numbers got as large as 50. or 75.
<br />
<br />Cathy
------------------------------
PosterpjrPostedMon Nov 01, 2004 6:06 pm MST, http://bb.cgd.ucar.edu//viewtopic.php?p=122#122
Message: Most were zero or on the order of 10^-15, but occasionally there were areas where the numbers got as large as 50. or 75.
<br />
<br />Cathy
<br />
<br />That is not good. Before I start worrying, or poking around in the code
<br />can you assure me that you have done the test "exactly" the way I
<br />suggested --- neither more, nor less. You must realize that we
<br />are juggling many different tasks here, and I dont want to be spending
<br />any time on sorting out other issues.
<br />
<br />Thanks
<br />
<br />Phil
------------------------------
Postercathryn.meyer at yale.eduPostedMon Nov 01, 2004 6:22 pm MST, http://bb.cgd.ucar.edu//viewtopic.php?p=123#123
Message: Here is the change that I made to tphysbc.F90 (replaced "cld" with "zero" in the radctl call):
<br />-----------------------------------------------------------------------------------------------
<br />! Complete radiation calculations
<br />!
<br /> call t_startf ('radctl')
<br /> call radctl (lchnk, ncol, lwup, emis, state%pmid, &
<br /> state%pint, state%lnpmid, state%lnpint, state%t, state%q, & zero, cicewp, cliqwp, coszrs, asdir, asdif, &
<br /> aldir, aldif, pmxrgn, nmxrgn, fsns, fsnt ,flns ,flnt , &
<br /> qrs, qrl, flwds, rel, rei, &
<br /> sols, soll, solsd, solld, &
<br /> landfrac, state%zm, state, fsds)
<br /> call t_stopf ('radctl')
<br />--------------------------------------------------------------------------------------------
<br /!
>
<br />I ran the model for 100 days in aquaplanet mode on my linux pc, although I don't think that would make any difference. I can send you ascii files of the SWCF and LWCF output if you'd like to see the numbers for yourself. Here are the maxes and mins:
<br />
<br />(0) SWCF max and min
<br />(0) 0.0001462979
<br />(0) -1.450159
<br />(0) LWCF max and min
<br />(0) 67.38367
<br />(0) -1.61599e-10
<br />
<br />Although most of the values are zero or near zero (10^-15 or less).
<br />
<br />I will try the model run again just in case I screwed anything up, but the output I just showed you was from the second try of the run already.
<br />
<br />Cathy
------------------------------
PosterpjrPostedMon Nov 01, 2004 8:40 pm MST, http://bb.cgd.ucar.edu//viewtopic.php?p=125#125
Message: Here is the change that I made to tphysbc.F90 (replaced "cld" with "zero" in the radctl call):
<br />-----------------------------------------------------------------------------------------------
<br />! Complete radiation calculations
<br />!
<br /> call t_startf ('radctl')
<br /> call radctl (lchnk, ncol, lwup, emis, state%pmid, &
<br /> state%pint, state%lnpmid, state%lnpint, state%t, state%q, & zero, cicewp, cliqwp, coszrs, asdir, asdif, &
<br /> aldir, aldif, pmxrgn, nmxrgn, fsns, fsnt ,flns ,flnt , &
<br /> qrs, qrl, flwds, rel, rei, &
<br /> sols, soll, solsd, solld, &
<br /> landfrac, state%zm, state, fsds)
<br /> call t_stopf ('radctl')
<br />--------------------------------------------------------------------------------------------
<br /!
>
<br />I ran the model for 100 days in aquaplanet mode on my linux pc, although I don't think that would make any difference. I can send you ascii files of the SWCF and LWCF output if you'd like to see the numbers for yourself. Here are the maxes and mins:
<br />
<br />(0) SWCF max and min
<br />(0) 0.0001462979
<br />(0) -1.450159
<br />(0) LWCF max and min
<br />(0) 67.38367
<br />(0) -1.61599e-10
<br />
<br />Although most of the values are zero or near zero (10^-15 or less).
<br />
<br />I will try the model run again just in case I screwed anything up, but the output I just showed you was from the second try of the run already.
<br />
<br />Cathy
<br />Stay tuned. We will get back to you.
<br />Phil
------------------------------
PosterpjrPostedMon Nov 01, 2004 8:47 pm MST, http://bb.cgd.ucar.edu//viewtopic.php?p=126#126
Message: My fault. The array "zero" is not a 2-D array, eg it was dimensioned
<br />(pcols), rather than (pcols,pver).
<br />
<br />Declare a 2-D array, like
<br />
<br />real(r8) zero2d(pcols,pver)
<br />
<br />set it to zero
<br />
<br />zero2d = 0
<br />
<br />and use that array instead of cld in the call to radctl. It worked for us.
<br />
<br />Phil
------------------------------
<< Using SCAM with specified surface properties, http://bb.cgd.ucar.edu//viewtopic.php?t=42 >>
Posterrneale_ca,PostedMon Nov 01, 2004 8:00 pm MST, http://bb.cgd.ucar.edu//viewtopic.php?p=124#124
Message: Is the option to use specified surface properties actually implemented in SCAM? Selecting the 'use specified surface properties' in the UI doesn't seem to have any effect and when I try and trace the logical variable use_srfprop it doesn't seem to get used anywhere after being initialized. I am using SCAM from the CCSM3.0 tarball. Thanks Rich.
------------------------------
Format TEXT
Show Message Text YES
Show My Messages YES
Digest Frequency DAY
Show only new messages since last time I logged in FALSE
Send empty digests NO
Time of day to send digest 9 AM
Maximum characters per message in digest 32000
More information about the cam-users
mailing list