Which variables need to be output from CAM

[MarionBWeinzierl]

Need to find out which variables we need to output from CAM for Cloud Liquid and Cloud Ice.

[dorchard]

(was that CLUBB variables to be output?)
AR by @paogorman @judithberner

[paogorman]

Expanding on this: I had previously asked Jack to also output the Clubb tendencies for comparison and he asked which exactly.

…

On Wed, May 1, 2024 at 10:44 AM Dominic Orchard ***@***.***> wrote: AR by @paogorman <https://github.com/paogorman> @judithberner <https://github.com/judithberner> — Reply to this email directly, view it on GitHub <#11 (comment)>, or unsubscribe <https://github.com/notifications/unsubscribe-auth/ARPLW62ZHP3FRKIDDNTJSLDZAD5THAVCNFSM6AAAAABHCCGK2KVHI2DSMVQWIX3LMV43OSLTON2WKQ3PNVWWK3TUHMZDAOBYGU3DSMBRGU> . You are receiving this because you were mentioned.Message ID: ***@***.***>

[paogorman]

These are the additional tendencies I think we need to output to better compare to the NN scheme. There are CLUBB and microphysics which could be added together. Also I think we might have been missing part of ZM due to evaporation of rain (this could be added to the other ZM tendencies). ZM tendencies: In addition to: *call* addfld ('ZMDT', (/ 'lev' /), 'A', 'K/s','T tendency - Zhang-McFarlane moist convection') *call* addfld ('ZMDQ', (/ 'lev' /), 'A', 'kg/kg/s','Q tendency - Zhang-McFarlane moist convection') *call* addfld ('ZMDICE', (/ 'lev' /), 'A', 'kg/kg/s','Cloud ice tendency - Zhang-McFarlane convection') *call* addfld ('ZMDLIQ', (/ 'lev' /), 'A', 'kg/kg/s','Cloud liq tendency - Zhang-McFarlane convection') I think we also need to add in tendencies related to evaporation of precipitation (this is added as a separate tendency in ZM for some reason): *call* addfld ('EVAPTZM', (/ 'lev' /), 'A', 'K/s','T tendency - Evaporation/snow prod from Zhang convection') *call* addfld ('EVAPQZM', (/ 'lev' /), 'A', 'kg/kg/s','Q tendency - Evaporation from Zhang-McFarlane moist convection') CLUBB tendencies: call addfld ('RVMTEND_CLUBB', (/ 'lev' /), 'A', 'g/kg /s', 'Water vapor tendency') call addfld ('STEND_CLUBB', (/ 'lev' /), 'A', 'k/s', 'Temperature tendency') call addfld ('RCMTEND_CLUBB', (/ 'lev' /), 'A', 'g/kg /s', 'Cloud Liquid Water Tendency') call addfld ('RIMTEND_CLUBB', (/ 'lev' /), 'A', 'g/kg /s', 'Cloud Ice Tendency') *call* addfld ('TTENDICE', (/ 'lev' /), 'A', 'K/s', 'T tendency from Ice Saturation Adjustment') *call* addfld ('QVTENDICE', (/ 'lev' /), 'A', 'kg/kg/s', 'Q tendency from Ice Saturation Adjustment') *call* addfld ('QCTENDICE', (/ 'lev' /), 'A', 'kg/kg/s', 'CLDICE tendency from Ice Saturation Adjustment') Microphysics tendencies *call* addfld ('MPDT', (/ 'lev' /), 'A', 'W/kg', 'Heating tendency - Morrison microphysics' ) *call* addfld ('MPDQ', (/ 'lev' /), 'A', 'kg/kg/s', 'Q tendency - Morrison microphysics' ) *call* addfld ('MPDLIQ', (/ 'lev' /), 'A', 'kg/kg/s', 'CLDLIQ tendency - Morrison microphysics' ) *call* addfld ('MPDICE', (/ 'lev' /), 'A', 'kg/kg/s', 'CLDICE tendency - Morrison microphysics'

…

On Wed, May 1, 2024 at 10:45 AM Paul O'Gorman ***@***.***> wrote: Expanding on this: I had previously asked Jack to also output the Clubb tendencies for comparison and he asked which exactly. On Wed, May 1, 2024 at 10:44 AM Dominic Orchard ***@***.***> wrote: > AR by @paogorman <https://github.com/paogorman> @judithberner > <https://github.com/judithberner> > > — > Reply to this email directly, view it on GitHub > <#11 (comment)>, > or unsubscribe > <https://github.com/notifications/unsubscribe-auth/ARPLW62ZHP3FRKIDDNTJSLDZAD5THAVCNFSM6AAAAABHCCGK2KVHI2DSMVQWIX3LMV43OSLTON2WKQ3PNVWWK3TUHMZDAOBYGU3DSMBRGU> > . > You are receiving this because you were mentioned.Message ID: > ***@***.***> >

[jatkinson1000]

To follow up here.
The variables currently output are as follows:

mfilt=1440,1440,1440
nhtfrq=1,1,1
fincl2 = 'ZMDT:I', 'ZMDQ:I', 'ZMDICE:I', 'ZMDLIQ:I', 'EVAPTZM:I', 'EVAPQZM:I', 'PRECCDZM:I',
         'RVMTEND_CLUBB:I', 'STEND_CLUBB:I', 'RCMTEND_CLUBB:I', 'RIMTEND_CLUBB:I',
         'TTENDICE:I', 'QVTENDICE:I', 'QCTENDICE:I',
         'MPDT:I', 'MPDQ:I', 'MPDLIQ:I', 'MPDICE:I',
         'YOGDT:I', 'YOGDQ:I', 'YOGDICE:I', 'YOGDLIQ:I', 'PRECC:I', 'YOGPREC:I',
         'RELHUM:I'
fincl3 = 'PMID:I','PS:I','T:I','TBP:I','Q:I','QBP:I','CLDICE&IC','CLDLIQ&IC','ICWMR:I','ICIMR:I'

With recent changes to run YOG as a separate scheme with deep_scheme='off' rather than as deep_scheme='YOG' the ZM setup routines will no longer be called so we will not be able to output the ZM variables from a YOG run.

With reference to @MarionBWeinzierl's original query on this issue.
This arose because cloud liquid and ice were being written as zero values from CLDLIQ and CLDICE.
Judith was going to look up why this might be and if these were not the correct variables to look at.
After some digging myself I am now using:

        float ICIMR(time, lev, lat, lon) ;
                ICIMR:mdims = 1 ;
                ICIMR:units = "kg/kg" ;
                ICIMR:long_name = "Prognostic in-cloud ice mixing ratio" ;
                ICIMR:cell_methods = "time: mean" ;
        float ICWMR(time, lev, lat, lon) ;
                ICWMR:mdims = 1 ;
                ICWMR:units = "kg/kg" ;
                ICWMR:long_name = "Prognostic in-cloud water mixing ratio" ;
                ICWMR:cell_methods = "time: mean" ;

instead.
I have asked @judithberner to check that these were the appropriate variables on 2024-06-04 but haven't yet received confirmation.