search

NCAR / amwg_dev

Repo to store model sandboxes and cases used for CAM development
9 stars 2 forks source link

f.e21.FWscHIST.ne30_L48_BL10_cam6_3_041_control.hf.001 #53

Open cecilehannay opened 2 years ago

cecilehannay commented 2 years ago

Description: New control with cam6_3_041 as requested in: https://github.com/NCAR/amwg_dev/issues/52

Case directory: Locally (if still available): /glade/p/cesmdata/cseg/runs/cesm2_0/f.e21.FWscHIST.ne30_L48_BL10_cam6_3_041_control.hf.001

On github: https://github.com/NCAR/amwg_dev/tree/f.e21.FWscHIST.ne30_L48_BL10_cam6_3_041_control.hf.001

Sandbox: Locally (if still available): /glade/work/hannay/cesm_tags/cam6_3_041

On github: https://github.com/ESCOMP/CAM/tree/cam6_3_041 hash: 963ebdc

Diagnostics: AMWG diags (if available) https://webext.cgd.ucar.edu/FWscHIST/f.e21.FWscHIST.ne30_L48_BL10_cam6_3_041_control.hf.001/atm/

Contacts: @cecilehannay @JulioTBacmeister @adamrher

Settings

This run uses:

CAM_CONFIG_OPTS="-phys cam_dev -age_of_air_trcs -chem waccm_sc_mam4 -nlev 58"

se_rsplit = 3 se_nsplit = 2 se_hypervis_subcycle = 4

ncdata = '/glade/p/cesm/amwg_dev/juliob/FWsc_ne30pg3_58L_GRID_48_taperstart10km_lowtop_BL10_v3_beta1p75_Top_43km.nc'

For full setting, check: https://github.com/NCAR/amwg_dev/tree/f.e21.FWscHIST.ne30_L48_BL10_cam6_3_041_control.hf.001

islasimpson commented 2 years ago

@cecilehannay I'm looking for the higher frequency data for this run as it is the control for @dan800 's comparison of the TEM diagnostics but it seems like only the h0 files are in

/glade/scratch/hannay/archive/f.e21.FWscHIST.ne30_L48_BL10_cam6_3_041_control.hf.001/atm/hist/

Don't suppose they have been moved to somewhere else? I'm guessing maybe the higher frequency data has been scrubbed if it hasn't been touched.

cecilehannay commented 2 years ago

Unfortunately, this run was done in January and the higher frequency has been scrubbed.

+++++++++++++++++++++++++++++++++++++++++++++++++++++++++ Cecile Hannay National Center for Atmospheric Research email: @.*** webpage: http://www.cgd.ucar.edu/staff/hannay/

"Live as if you were to die tomorrow. Dream as if you were to live forever." +++++++++++++++++++++++++++++++++++++++++++++++++++++++++

On Fri, Jul 1, 2022 at 10:35 AM islasimpson @.***> wrote:

@cecilehannay https://github.com/cecilehannay I'm looking for the higher frequency data for this run as it is the control for @dan800 https://github.com/dan800 's comparison of the TEM diagnostics but it seems like only the h0 files are in

/glade/scratch/hannay/archive/f.e21.FWscHIST.ne30_L48_BL10_cam6_3_041_control.hf.001/atm/hist/

Don't suppose they have been moved to somewhere else? I'm guessing maybe the higher frequency data has been scrubbed if it hasn't been touched.

— Reply to this email directly, view it on GitHub https://github.com/NCAR/amwg_dev/issues/53#issuecomment-1172521930, or unsubscribe https://github.com/notifications/unsubscribe-auth/ACKF2VA6FLOWFA7GVPNHPNTVR4M5JANCNFSM5LXEDGZQ . You are receiving this because you were mentioned.Message ID: @.***>

islasimpson commented 2 years ago

Ok, thanks

dan800 commented 2 years ago

@cecilehannay Is f.e21.FWscHIST.ne30_L48_BL10_cam6_3_041_control still our baseline from which all improvements are being based or have we moved forward? If it is we really need to redo it so we can diagnose the baseline TEM. If it's not then is there high frequency output from the new baseline to analyze?

cecilehannay commented 2 years ago

Here is the new baseline. /glade/scratch/hannay/archive/f.cesm3_cam058_mom_b.FWscHIST.ne30_L58.001

with high frequency output.

dan800 commented 2 years ago

@islasimpson @JulioTBacmeister What would be the best comparison? To run a L92 with this exact configuration (identical physics) but dumping GW momentum at lid. Presumably this would have a better residual circulation. We wanted to diagnose if we can get away with zero flux UBC in that configuration anyway.

islasimpson commented 2 years ago

Here is the new baseline. /glade/scratch/hannay/archive/f.cesm3_cam058_mom_b.FWscHIST.ne30_L58.001 with high frequency output.

Thanks. I'm looking at this though and there are no .h2. files. Were they not saved for this one? If so, I think I can calculate the TEM diagnostics using the daily means for this one and the 6 hourly for the other one. That's ok.

cecilehannay commented 2 years ago

Yes, this correct. We dropped the 6-hourly output based on this document. We don't output h2 files anymore.

dan800 commented 2 years ago

Are we still OK with calculating TEM from daily or do we need to revisit the default output for this stage of model evaluation? If revisiting, then @adamrher -- can you see if we still need all those microphysical quantities on every run. What about N2O.

adamrher commented 2 years ago

I don't believe N2O is a dynamic tracer in the compsets that Cecile has been running. We are evaluating N2O in CAM-CHEM runs, separately.

I assume you are referring to the monthly microphysics/aerosol fields... this stuff?

'ACTREL','ACTNL','ACTREI','FCTL','FCTI','ACTNI','BURDENSO4','BURDENBC','BURDENDUST',
'BURDENSOA','BURDENSEASALT','BURDENPOM','CCN3','CDNUMC','FICE','PBLH','RELHUM',

Doesn't seem like too much trouble to keep these around. Especially as we start to evaluate MAM5.

dan800 commented 2 years ago

re: N2O -- This is where things get confusing and I still think we need to consolidate. I think in SC (so L92) configurations N2O is a dynamic tracer. It sounds like it might be useful to have it in both L58 and L92 if it can help diagnose tracer transport issues.

That CAMP7 doc has a lot of output under monthly (some repeated) and it's not clear if the ISSCP stuff is defaulted to be output or not ((that's more of what I was referring to).

adamrher commented 2 years ago

@cecilehannay are we turning on the isscp diagnostics by default?

you're right, Dan. N2O is advected by the dycore in Cecile's L58 runs. I just double checked the log file. So we should ass N2O to the monthly files. We could do this for all chemical species. I think the BURDEN variables are enough to summarize the behavior advected aerosol species ... we should then consider whether to add the chemical species.

dan800 commented 2 years ago

So I looked at the google sheet, which is supposed to supercede the google doc - N2O isn't listed, even as an optional output. CH4 is listed as optional but I think is not optional - it's hugely valuable for understanding strat transport and water vapor. If CFCs are still being transported then they need to be added to the CAM7 output options.

adamrher commented 2 years ago

not familiar with the spreadsheet, but just looking at user_nl_cam for her latest run here, it does appear CH4 is included in her default i/o. That is the only chemical species included. The rest of the advected chemical species:

    11  N2O       N2O                                                                                                                               dry
    12  CFC11     CFC11                                                                                                                             dry
    13  CFC12     CFC12                                                                                                                             dry
    14  H2O2      H2O2                                                                                                                              dry
    15  H2SO4     H2SO4                                                                                                                             dry
    16  SO2       SO2                                                                                                                               dry
    17  DMS       DMS                                                                                                                               dry
    18  SOAG      SOAG                                                                                                                              dry

Should be added to the monthly fields.

cecilehannay commented 2 years ago

Yes, the google sheet supersedes the google doc. I was referring to the google doc because this is where we had the discussion to forgo the 6-hourly output.

The latest version of the fincl is:

nhtfrq         =                  0, -24,  -6,  -3,  -1,   1, -24,-24

fincl1 = 'ACTNI', 'ACTNL', 'ACTREI', 'ACTREL', 'AODDUST', 'AODVIS', 'BURDENBC', 'BURDENDUST', 'BURDENPOM', 'BURDENSEASALT', 
'BURDENSO4', 'BURDENSOA', 'CAPE', 'CCN3', 'CDNUMC', 'CH4', 'CLDHGH', 'CLDICE', 'CLDLIQ', 'CLDLOW', 'CLDMED', 'CLDTOT', 'CLOUD', 'CMFMC_DP', 
'CT_H2O', 'DCQ', 'DQCORE', 'DTCOND', 'DTCORE', 'DTV', 'EVAPPREC', 'EVAPSNOW', 'FCTI', 'FCTL', 'FICE', 'FLDS', 'FLNS', 'FLNSC', 'FLNT', 'FLNTC', 'FLUT', 
'FREQZM', 'FSDS', 'FSDSC', 'FSNS', 'FSNSC', 'FSNT', 'FSNTC', 'FSNTOA', 'ICEFRAC', 'LANDFRAC', 'LHFLX', 'LWCF', 'MPDICE', 'MPDLIQ', 'MPDQ', 'MPDT', 
'OCNFRAC', 'OMEGA', 'OMEGA500', 'PBLH', 'PHIS', 'PINT', 'PMID', 'PRECC', 'PRECL', 'PRECSC', 'PRECSL', 'PRECT', 'PS', 'PSL', 'PTEQ', 'PTTEND', 'Q', 
'QFLX', 'QRL', 'QRS', 'QTGW', 'RCMTEND_CLUBB', 'RELHUM', 'RVMTEND_CLUBB', 'SHFLX', 'SOLIN', 'SST', 'STEND_CLUBB', 'SWCF', 
'T', 'TAUX', 'TAUY', 'TFIX', 'TGCLDIWP', 'TGCLDLWP', 'TMQ', 'TREFHT', 'TS', 'TTGW', 'U', 'U10', 'UBOT', 'UTGWORO', 'UTGW_TOTAL', 
'V', 'VBOT', 'VTGWORO', 'VTGW_TOTAL', 'WPRTP_CLUBB', 'WPTHLP_CLUBB', 'Z3', 'ZMDQ', 'ZMDT'

fincl2 = 'OMEGA', 'PMID', 'PS', 'Q', 'QRL', 'QRS', 'T', 'TROP_P', 'TROP_T', 'U', 'V', 'Z3'

fincl4 =  'PRECC','PRECL'

I will add the variables Adam suggests in the next run.