KeyClim Cloud2 tunings

Externals_CAM.cfg

@@ -2,7 +2,7 @@
 local_path = src/chemistry/oslo_aero
 protocol = git
 repo_url = https://github.com/NorESMhub/OSLO_AERO.git
-branch = noresm2_3_develop
+tag = oslo_aero_2_3a1
 required = True
 
 [externals_description]

bld/build-namelist

@@ -2266,7 +2266,7 @@ if ($chem eq 'trop_mam_oslo' ) {
     unless (defined $nl->get_value('srf_emis_type')) {
         add_default($nl, 'srf_emis_type',      'val'=>'CYCLICAL');
         add_default($nl, 'srf_emis_cycle_yr',  'val'=>'2000');
-    } 
+    }
 
     # Vertical emission datasets:
     %species = ();
@@ -2350,7 +2350,7 @@ if ($chem =~ /_mam4/ and $phys =~ /cam6/) {
                 'pom_a4_an_srf_file'  => 'pom_a4',
                 'pom_a4_bb_srf_file'  => 'pom_a4' );
     }
-	    
+
     # for mechanism missing full tropospheric chemistry
     if ($chem =~ /trop_mam/ or $chem =~ /waccm_ma/ or $chem =~ /waccm_sc/) {
         %species = (%species,
@@ -2905,7 +2905,9 @@ unless (defined $nl->get_value('dms_source_type')) {
 #add_default($nl, 'opom_cycle_year');
 #--djlo--
 
-if ($chem =~ /_mam_oslo/) {$nl->set_variable_value('phys_ctl_nl','use_hetfrz_classnuc','.true.');}
+if ($chem =~ /_mam_oslo/) {
+  add_default($nl, 'hetfrz_aer_scalfac');
+}
 
 if ($waccmx) {
   my $wmx_opt = $nl->get_value('waccmx_opt');
@@ -2944,8 +2946,8 @@ if ($waccmx) {
     add_default($nl,'wei05_coefs_file');
     add_default($nl,'solar_wind_data_file');
   }
-  add_default($nl,'ionos_xport_nsplit');  
-  add_default($nl,'steady_state_ion_elec_temp', 'val'=>'.false.');  
+  add_default($nl,'ionos_xport_nsplit');
+  add_default($nl,'steady_state_ion_elec_temp', 'val'=>'.false.');
 }
 
 # Chemistry options
@@ -3135,6 +3137,10 @@ if ($clubb_sgs  =~ /$TRUE/io) {
    add_default($nl, 'clubb_l_lscale_plume_centered');
    add_default($nl, 'clubb_l_use_ice_latent');
    add_default($nl, 'clubb_do_liqsupersat');
+   add_default($nl, "clubb_meltpt_temp");
+
+   add_default($nl, "clubb_meltpt_temp");
+   add_default($nl, "clubb_dt_low");
 }
 
 # Tuning for wet scavenging of modal aerosols

bld/namelist_files/namelist_defaults_cam.xml

@@ -1025,7 +1025,7 @@
 
 <!-- prescribed oxidants for prognostic aerosol chemistry -->
 <tracer_cnst_file     ver="fixed_ox">oxid_1.9x2.5_L26_1850-2005_c091123.nc</tracer_cnst_file>
-                                              
+
 <!-- prescribed oxidants for prognostic aerosol chemistry -lightweight set for running SCAM and producing SCAM IOP boundary data -->
 <tracer_cnst_file     phys="cam6"  ver="fixed_ox" scam="1">oxid_1.9x2.5_L26_1850clim_c091123.nc</tracer_cnst_file>
 <tracer_cnst_cycle_yr phys="cam6"  ver="fixed_ox" scam="1">1850</tracer_cnst_cycle_yr>
@@ -1196,9 +1196,14 @@
 <clubb_mult_coef                                  > 1.0D0   </clubb_mult_coef>
 <clubb_lambda0_stability_coef                     > 0.04    </clubb_lambda0_stability_coef>
 <clubb_l_lscale_plume_centered                    > .false. </clubb_l_lscale_plume_centered>
-<clubb_l_use_ice_latent                           > .false.  </clubb_l_use_ice_latent>
+<clubb_l_use_ice_latent                           > .false. </clubb_l_use_ice_latent>
 <clubb_do_liqsupersat                             > .false. </clubb_do_liqsupersat>
 
+<!-- CLUBB NorESM -->
+<clubb_meltpt_temp           >268.15D0</clubb_meltpt_temp>
+<clubb_meltpt_temp camnor="1">243.15D0</clubb_meltpt_temp>
+<clubb_dt_low                >238.15D0</clubb_dt_low>
+<clubb_dt_low camnor="1"     >238.15D0</clubb_dt_low>
 
 <!-- Microphysics scheme -->
 <microp_scheme                          >NONE         </microp_scheme>
@@ -1244,8 +1249,12 @@
 <nucleate_ice_subgrid_strat microphys="mg1"            >1.2D0</nucleate_ice_subgrid_strat>
 <nucleate_ice_subgrid_strat microphys="mg2"            >1.2D0</nucleate_ice_subgrid_strat>
 
-<use_hetfrz_classnuc            >.false.</use_hetfrz_classnuc>
-<use_hetfrz_classnuc phys="cam6">.true.</use_hetfrz_classnuc>
+<use_hetfrz_classnuc                     >.false.</use_hetfrz_classnuc>
+<use_hetfrz_classnuc phys="cam6"         >.true.</use_hetfrz_classnuc>
+<use_hetfrz_classnuc chem="trop_mam_oslo">.true.</use_hetfrz_classnuc>
+
+<hetfrz_aer_scalfac           >1.0D0</hetfrz_aer_scalfac>
+<hetfrz_aer_scalfac camnor="1">0.001D0</hetfrz_aer_scalfac>
 
 <use_preexisting_ice            >.false.</use_preexisting_ice>
 <use_preexisting_ice phys="cam6">.true.</use_preexisting_ice>

bld/namelist_files/namelist_definition.xml

@@ -65,41 +65,41 @@
         Nudging tendencies are applied as a relaxation force between the current
         model state values and target state values derived from the avalilable
         analyses. The form of the target values is selected by the 'Nudge_Force_Opt'
-        option, the timescale of the forcing is determined from the given 
-        'Nudge_TimeScale_Opt', and the nudging strength Alpha=[0.,1.] for each 
+        option, the timescale of the forcing is determined from the given
+        'Nudge_TimeScale_Opt', and the nudging strength Alpha=[0.,1.] for each
         variable is specified by the 'Nudge_Xcoef' values. Where X={U,V,T,Q,PS}
- 
+
                F_nudge = Alpha*((Target-Model(t_curr))/TimeScale
- 
+
         WINDOWING:
         ----------
-        The region of applied nudging can be limited using Horizontal/Vertical 
-        window functions that are constructed using a parameterization of the 
-        Heaviside step function. 
-    
-        The Heaviside window function is the product of separate horizonal and vertical 
+        The region of applied nudging can be limited using Horizontal/Vertical
+        window functions that are constructed using a parameterization of the
+        Heaviside step function.
+
+        The Heaviside window function is the product of separate horizonal and vertical
         windows that are controled via 12 parameters:
-    
-            Nudge_Hwin_lat0:     Specify the horizontal center of the window in degrees. 
-            Nudge_Hwin_lon0:     The longitude must be in the range [0,360] and the 
+
+            Nudge_Hwin_lat0:     Specify the horizontal center of the window in degrees.
+            Nudge_Hwin_lon0:     The longitude must be in the range [0,360] and the
                                  latitude should be [-90,+90].
-            Nudge_Hwin_latWidth: Specify the lat and lon widths of the window as positive 
-            Nudge_Hwin_lonWidth: values in degrees.Setting a width to a large value (e.g. 999) 
+            Nudge_Hwin_latWidth: Specify the lat and lon widths of the window as positive
+            Nudge_Hwin_lonWidth: values in degrees.Setting a width to a large value (e.g. 999)
                                  renders the window a constant in that direction.
-            Nudge_Hwin_latDelta: Controls the sharpness of the window transition with a 
-            Nudge_Hwin_lonDelta: length in degrees. Small non-zero values yeild a step 
+            Nudge_Hwin_latDelta: Controls the sharpness of the window transition with a
+            Nudge_Hwin_lonDelta: length in degrees. Small non-zero values yeild a step
                                  function while a large value yeilds a smoother transition.
-            Nudge_Hwin_Invert  : A logical flag used to invert the horizontal window function 
+            Nudge_Hwin_Invert  : A logical flag used to invert the horizontal window function
                                  to get its compliment.(e.g. to nudge outside a given window).
-    
-            Nudge_Vwin_Lindex:   In the vertical, the window is specified in terms of model 
-            Nudge_Vwin_Ldelta:   level indcies. The High and Low transition levels should 
-            Nudge_Vwin_Hindex:   range from [0,(NLEV+1)]. The transition lengths are also 
-            Nudge_Vwin_Hdelta:   specified in terms of model indices. For a window function 
+
+            Nudge_Vwin_Lindex:   In the vertical, the window is specified in terms of model
+            Nudge_Vwin_Ldelta:   level indcies. The High and Low transition levels should
+            Nudge_Vwin_Hindex:   range from [0,(NLEV+1)]. The transition lengths are also
+            Nudge_Vwin_Hdelta:   specified in terms of model indices. For a window function
                                  constant in the vertical, the Low index should be set to 0,
-                                 the High index should be set to (NLEV+1), and the transition 
-                                 lengths should be set to 0.001 
-            Nudge_Vwin_Invert  : A logical flag used to invert the vertical window function 
+                                 the High index should be set to (NLEV+1), and the transition
+                                 lengths should be set to 0.001
+            Nudge_Vwin_Invert  : A logical flag used to invert the vertical window function
                                  to get its compliment.
        Default: FALSE
        </entry>
@@ -171,16 +171,16 @@
 <entry id="Nudge_Force_Opt" type="integer" category="nudging"
        group="nudging_nl" valid_values="" >
        Select the form of nudging forcing, where (t'==Analysis times ; t==Model Times)
-         0 -> NEXT-OBS: Target=Anal(t'_next)                 
+         0 -> NEXT-OBS: Target=Anal(t'_next)
          1 -> LINEAR:   Target=(F*Anal(t'_curr) +(1-F)*Anal(t'_next))
                             F =(t'_next - t_curr )/Tdlt_Anal
        Default: 0
        </entry>
 
 <entry id="Nudge_TimeScale_Opt" type="integer" category="nudging"
        group="nudging_nl" valid_values="" >
-       Select the timescale of nudging force, where (t'==Analysis times ; t==Model Times) 
-         0 -->  TimeScale = 1/Tdlt_Anal                     
+       Select the timescale of nudging force, where (t'==Analysis times ; t==Model Times)
+         0 -->  TimeScale = 1/Tdlt_Anal
          1 -->  TimeScale = 1/(t'_next - t_curr )
        Default: 0
        </entry>
@@ -2647,6 +2647,12 @@ Add diagnostic output for heterogeneous freezing code.
 Default: .false.
 </entry>
 
+<entry id="hetfrz_aer_scalfac" type="real" category="microphys"
+       group="hetfrz_classnuc_nl" valid_values="" >
+Scaling factor for aerosols
+Default: 0.001
+</entry>
+
 <entry id="use_preexisting_ice" type="logical" category="microphys"
        group="nucleate_ice_nl" valid_values="" >
 Switch to turn on treatment of pre-existing ice in the ice nucleation code.
@@ -3325,6 +3331,20 @@ energy.
 Default: true
 </entry>
 
+<!-- CLUBB NorESM options -->
+<entry id="clubb_meltpt_temp" type="real" category="NorESM"
+       group="clubb_noresm_nl" valid_values="">
+Temperature used for the melting temp of ice crystals [K]
+Default: 268.15K
+</entry>
+
+<entry id="clubb_dt_low" type="real" category="NorESM"
+       group="clubb_noresm_nl" valid_values="">
+Temperature at which detrained water is classified as entirely ice (no
+liquid) in the CLUBB parameterization in units of (K).
+Default: 238.15K
+</entry>
+
 <!-- CARMA Sectional Microphysics -->
 
 <entry id="carma_model" type="char*32" category="carma"
@@ -4113,7 +4133,7 @@ Default: .false.
 <entry id="dme_energy_adjust" type="logical"  category="NorESM"
        group="phys_ctl_nl" valid_values="" >
 Switch to use appropriate energy adjustment in dry-mass adjustment at the
-end of tphysac. 
+end of tphysac.
 Default: .false.
 </entry>
 
@@ -7724,7 +7744,7 @@ Default: none
 
 <entry id="volc_fraction_coarse" type="real" category="cam_oslo"
        group="oslo_ctl_nl" valid_values="" >
-Fraction of volcanic aerosols which will end up in coarse mode 
+Fraction of volcanic aerosols which will end up in coarse mode
 Default: 0.0
 </entry>
 <entry id="aerotab_table_dir" type="char*256" input_pathname="abs" category="cam_oslo"

cime_config/config_component.xml

@@ -8,7 +8,7 @@
       CAM
     ===============
    -->
-    <desc atm="CAM60[%1PCT][%4xCO2][%2xCO2][%CCTS][%CFIRE][%CVBSX][%PORT][%RCO2][%SCAM][%SDYN][%WCCM][%WCMD][%WCSC][%WCTS][%NORESM][%FRC2][%GHGONLY][%NATONLY][%AEROXIDONLY][%AERONLY][%OXIDONLY][%OZONEONLY][%LUONLY][%PINTCF][%PIAER][%PIAEROXID][%NORBC][%GHG2014][%GHGNOH2O2014][%CO22014][%N2O2014][%CH42014][%CH4NOH2O2014][%BC2014][%OC2014][%SO22014][%AER2014][%AEROXID2014][%SO2OXID2014][%NTCF2014][%ANTHRO2014][%GHGOZONELU2014][%OZONE2014][%H2O2014][%OXID2014][%FSST][%NORPDDMSBC][%NORPIBC][%SSP245][%AERLOW][%FRC2EXT][%COVBASLIN][%COVFOSFUE][%COVMODGRE][%COVSTRGRE][%COVTWOBLI]">CAM cam6 physics:</desc>
+    <desc atm="CAM60[%1PCT][%4xCO2][%2xCO2][%CCTS][%CFIRE][%CVBSX][%PORT][%RCO2][%SCAM][%SDYN][%WCCM][%WCMD][%WCSC][%WCTS][%NORESM][%FRC2][%CLOUD2][%GHGONLY][%NATONLY][%AEROXIDONLY][%AERONLY][%OXIDONLY][%OZONEONLY][%LUONLY][%PINTCF][%PIAER][%PIAEROXID][%NORBC][%GHG2014][%GHGNOH2O2014][%CO22014][%N2O2014][%CH42014][%CH4NOH2O2014][%BC2014][%OC2014][%SO22014][%AER2014][%AEROXID2014][%SO2OXID2014][%NTCF2014][%ANTHRO2014][%GHGOZONELU2014][%OZONE2014][%H2O2014][%OXID2014][%FSST][%NORPDDMSBC][%NORPIBC][%SSP245][%AERLOW][%FRC2EXT][%COVBASLIN][%COVFOSFUE][%COVMODGRE][%COVSTRGRE][%COVTWOBLI]">CAM cam6 physics:</desc>
     <desc atm="CAM50[%CCTS][%CLB][%PORT][%RCO2][%SCAM][%SDYN][%WCSC][%WCTS]"   >CAM cam5 physics:</desc>
     <desc atm="CAM40[%PORT][%RCO2][%SCAM][%SDYN][%TMOZ][%WXIE][%WXIED][%WCMD]" >CAM cam4 physics:</desc>
     <desc atm="CAM[%ADIAB][%DABIP04][%TJ16][%HS94][%KESSLER][%RCO2]">CAM simplified and non-versioned physics :</desc>
@@ -22,7 +22,7 @@
     <desc option="2xCO2" >abrupt doubling of CO2 with other forcings maintained at 1850 piControl levels :</desc>
     <desc option="1PCT"  >ramped CO2 increasing by 1% per year with other forcings maintained at 1850 piControl levels (CMIP6 DECK 1pctCO2 experiment) :</desc>
 
-    <desc atm="CAM60%NORESM" >cam 6 and general NorESM  changes + Production tagged aerosols (OSLO_AERO)</desc>
+    <desc atm="CAM60%NORESM" >CAM 6, general NorESM changes, and production tagged aerosols (OSLO_AERO)</desc>
 
   <!--
     ===============
@@ -43,6 +43,13 @@
     <desc option="RCO2"  >CAM CO2 ramp: </desc>
     <desc option="TMOZ"  >CAM tropospheric chemistry with bulk aerosols:</desc>
 
+  <!--
+    ================
+     CAM-Nor Options
+    ================
+  -->
+  <desc option="CLOUD2">Runtime settings for the keyClim Cloud2 experiment"</desc>
+
   <!--
     ===============
      WACCM Options
@@ -532,6 +539,7 @@
     <value compset="DOCN%SOMAQP">$SRCROOT/components/cam/cime_config/usermods_dirs/aquap</value>
     <value compset="_CAM40%WX.*_CLM40" grid="a%1.9x2.5_l%1.9x2.5">$SRCROOT/components/cam/cime_config/usermods_dirs/waccmx</value>
     <value compset="CAM[456]0%SCAM" >$SRCROOT/components/cam/cime_config/usermods_dirs/scam_mandatory</value>
+    <value compset="CLOUD2">$SRCROOT/components/cam/cime_config/usermods_dirs/cmip6_noresm_keyCLIM_cloud2</value>
     </values>
     <group>run_component_cam</group>
     <file>env_case.xml</file>

cime_config/config_compsets.xml

@@ -60,9 +60,8 @@
 
   <!-- comment on the naming of NorESM compsets
 
-(1) CAM60%NORESM versus CAM60%PTAERO :
+(1) CAM60%NORESM
 CAM60%NORESM : the compsets used and created for CMIP6 use "NORESM" to activate typically NorESM2 settings
-CAM60%PTAERO : older compsets use "PTAERO" to activate typically NorESM2 settings
 
 (2) presence of frc2 (FRC2) in the compset name :
 WITHOUT "frc2" in name : uses multiple emission files per species (leads on fram HPC to irreprudicibility for fully-coupled compsets)

cime_config/usermods_dirs/cmip6_noresm_fsst_hifreq_xaer/user_nl_cam

@@ -1,9 +1,9 @@
 
 history_aerosol=.true.
 
-nhtfrq = 0, -24, -6, -6, -3,  -3 , -1,   1, -24 
+nhtfrq = 0, -24, -6, -6, -3,  -3 , -1,   1, -24
 
-mfilt  = 1,   73, 292, 292, 584, 584, 584, 240, 365 
+mfilt  = 1,   73, 292, 292, 584, 584, 584, 240, 365
 
 ndens  = 2,   2,  2,  2,   2, 2,   2,   1, 1
 
@@ -30,6 +30,9 @@ fincl6='UBOT:I','VBOT:I','TREFHT:I','QREFHT:I','TS:I','SST:I','PS:I','ICEFRAC:I'
 fincl7='PS:A','TREFHT:A','MMRPM2P5_SRF:A'
 
 
-do_circulation_diags           = .true. 
-
+do_circulation_diags           = .true.
 
+! CMIP6 settings for NorESM physics
+hetfrz_aer_scalfac = 1.0D0
+clubb_meltpt_temp = 268.15D0
+clubb_dt_low      = 238.15D0

cime_config/usermods_dirs/cmip6_noresm_fsst_xaer/user_nl_cam

@@ -1,9 +1,9 @@
 
 history_aerosol=.true.
 
-nhtfrq = 0, -24, -6, -3,  -1,   1, -24,-120,-240 
+nhtfrq = 0, -24, -6, -3,  -1,   1, -24,-120,-240
 
-mfilt  = 1,   5, 20, 40, 120, 240, 365,  73, 365 
+mfilt  = 1,   5, 20, 40, 120, 240, 365,  73, 365
 
 ndens  = 2,   2,  2,  2,   2,   2,   1,   1,   1
 
@@ -12,12 +12,15 @@ fincl1 = 'SST','TAUX','TAUY','TAUBLJX','TAUBLJY','BTAUNET','PRECC','PRECL','PREC
 
 
 fincl2 = 'ABSVIS:A', 'ACTNL:A', 'ACTREL:A', 'AOD_VIS:A', 'cb_BC:A', 'cb_DMS:A', 'cb_DUST:A', 'cb_OM:A', 'cb_SALT:A', 'cb_SO2:A', 'cb_SULFATE:A',
-   'CDNUMC:A', 'CLDICE:A', 'CLDLIQ:A', 'CLDTOT:A', 'CLOUD:A', 'CMFMC:A', 'CMFMCDZM:A', 'DAYFOC:A', 'FCTL:A', 'FLDS:A', 'FLDSC:A', 'FLNR:A', 'FLNS:A', 'FLNSC:A', 
-   'FLNT:A', 'FLNTC:A', 'FLUT:A', 'FLUTC:A', 'FSDS:A', 'FSDSC:A', 'FSNR:A', 'FSNS:A', 'FSNSC:A', 'FSNTOA:A', 'FSNTOAC:A', 'ICEFRAC:A' , 'LHFLX:A', 'MASS:A', 'OMEGA:A', 
-   'OMEGA500:A', 'PBLH:A', 'PDELDRY:A', 'PRECC:A', 'PRECT:A', 'PS:A', 'PSL:A', 'Q:A', 'QREFHT:A', 'QSNOW:A', 'RHREFHT:A', 'SHFLX:A', 
-   'SOLIN:A', 'SOLLD:A', 'SOLSD:A', 'SST:A' ,'T:A', 'T500:A', 'T700:A', 'T850:A', 'TAUBLJX:A', 'TAUBLJY:A', 'TAUGWX:A', 'TAUGWY:A', 'TAUX:A', 'TAUY:A', 
+   'CDNUMC:A', 'CLDICE:A', 'CLDLIQ:A', 'CLDTOT:A', 'CLOUD:A', 'CMFMC:A', 'CMFMCDZM:A', 'DAYFOC:A', 'FCTL:A', 'FLDS:A', 'FLDSC:A', 'FLNR:A', 'FLNS:A', 'FLNSC:A',
+   'FLNT:A', 'FLNTC:A', 'FLUT:A', 'FLUTC:A', 'FSDS:A', 'FSDSC:A', 'FSNR:A', 'FSNS:A', 'FSNSC:A', 'FSNTOA:A', 'FSNTOAC:A', 'ICEFRAC:A' , 'LHFLX:A', 'MASS:A', 'OMEGA:A',
+   'OMEGA500:A', 'PBLH:A', 'PDELDRY:A', 'PRECC:A', 'PRECT:A', 'PS:A', 'PSL:A', 'Q:A', 'QREFHT:A', 'QSNOW:A', 'RHREFHT:A', 'SHFLX:A',
+   'SOLIN:A', 'SOLLD:A', 'SOLSD:A', 'SST:A' ,'T:A', 'T500:A', 'T700:A', 'T850:A', 'TAUBLJX:A', 'TAUBLJY:A', 'TAUGWX:A', 'TAUGWY:A', 'TAUX:A', 'TAUY:A',
    'TGCLDIWP:A', 'TGCLDLWP:A', 'TMQ:A', 'TREFHT:A', 'TREFHTMN:M', 'TREFHTMX:X', 'TS:A', 'TSMN:M', 'TSMX:X', 'U10:A', 'UTGWORO:A',
     'Z500:A',
     'Z100:A','Z1000:A','DOD550:A','TAUTMODIS:A','CLTMODIS:A'
 
-
+! CMIP6 settings for NorESM physics
+hetfrz_aer_scalfac = 1.0D0
+clubb_meltpt_temp = 268.15D0
+clubb_dt_low      = 238.15D0

cime_config/usermods_dirs/cmip6_noresm_keyCLIM_cloud2/user_nl_cam

@@ -0,0 +1,30 @@
+
+history_aerosol=.true.
+
+nhtfrq = 0, -24, -6, -3,  -1,   1, -24,-120,-240
+
+mfilt  = 1,   5, 20, 40, 120, 240, 365,  73, 365
+
+ndens  = 2,   2,  2,  2,   2,   2,   1,   1,   1
+
+fincl1 = 'SST','TAUX','TAUY','TAUBLJX','TAUBLJY','BTAUNET','PRECC','PRECL','PRECT','FREQZM','PCONVB','PCONVT','PRECCDZM','Z700','Z500','Z200','Z300','Z100','Z050','U200','U850','V200','V850','T200','T500',   'T700','T1000','OMEGA500','OMEGA850','VTHzm','WTHzm','UVzm','UWzm','Uzm','Vzm','THzm','Wzm','dUzm','dVzm','dUazm','dVazm','dUfzm','U','V','T','Q','Z3','dU','dV','dUa','dVa','dUf','EFLX','PTTEND','IETEND_DME', 'PTTEND_DME','TFIX','EFIX','EP','QFLX','MEANPTOP','MEANTTOP','MEANTAU','TCLDAREA','RHREFHT','TREFMXAV','TREFMNAV','ozone','O3','TROP_P','TROP_T','TROP_Z','VT100'
+
+
+fincl2 = 'ABSVIS:A', 'ACTNL:A', 'ACTREL:A', 'AOD_VIS:A', 'cb_BC:A', 'cb_DMS:A', 'cb_DUST:A', 'cb_OM:A', 'cb_SALT:A', 'cb_SO2:A', 'cb_SULFATE:A',
+   'CDNUMC:A', 'CLDICE:A', 'CLDLIQ:A', 'CLDTOT:A', 'CLOUD:A', 'CMFMC:A', 'CMFMCDZM:A', 'DAYFOC:A', 'FCTL:A', 'FLDS:A', 'FLDSC:A', 'FLNR:A', 'FLNS:A', 'FLNSC:A',
+   'FLNT:A', 'FLNTC:A', 'FLUT:A', 'FLUTC:A', 'FSDS:A', 'FSDSC:A', 'FSNR:A', 'FSNS:A', 'FSNSC:A', 'FSNTOA:A', 'FSNTOAC:A', 'ICEFRAC:A' , 'LHFLX:A', 'MASS:A', 'OMEGA:A',
+   'OMEGA500:A', 'PBLH:A', 'PDELDRY:A', 'PRECC:A', 'PRECT:A', 'PS:A', 'PSL:A', 'Q:A', 'QREFHT:A', 'QSNOW:A', 'RHREFHT:A', 'SHFLX:A',
+   'SOLIN:A', 'SOLLD:A', 'SOLSD:A', 'SST:A' ,'T:A', 'T500:A', 'T700:A', 'T850:A', 'TAUBLJX:A', 'TAUBLJY:A', 'TAUGWX:A', 'TAUGWY:A', 'TAUX:A', 'TAUY:A',
+   'TGCLDIWP:A', 'TGCLDLWP:A', 'TMQ:A', 'TREFHT:A', 'TREFHTMN:M', 'TREFHTMX:X', 'TS:A', 'TSMN:M', 'TSMX:X', 'U10:A', 'UTGWORO:A',
+    'Z500:A'
+
+fincl4 = 'Q:I','QREFHT:I','Q850:I','PSL:I', 'PS:I','T:I', 'TREFHT:I','U:I','V:I','Z500:I'
+
+clubb_gamma_coef  = 0.330
+
+clubb_meltpt_temp = 243.15D0
+clubb_dt_low      = 238.15D0
+
+hetfrz_aer_scalfac = 0.001D0
+
+micro_mg_berg_eff_factor = 0.5D0