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 https://github.com/m2lines/CAM-ML/issues/11#issuecomment-2088569015, 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 https://github.com/m2lines/CAM-ML/issues/11#issuecomment-2088569015, 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.