Merge pull request NCAR#120 from pjpegion/pjpegion

Pjpegion

README.md

@@ -83,3 +83,15 @@ compiler environment variables.
   * When invoking the run script, the only required argument is the name of the case to run. The case name used must match one of the case configuration files located in ../etc/case_config (without the .nml extension!). If specifying a suite other than the default, the suite name used must match the value of the suite name in one of the suite definition files located in ../../ccpp/suites, (e.g. `SCM_GFS_v15`). If specifying a namelist other than the default, the value must be an entire filename that exists in ../../ccpp/physics_namelists. The -g flag can be used to run the executable through the gdb debugger (assuming it is installed on the system).
 2. A netcdf output file is generated in the location specified in the case
 configuration file (is present), or in an output directory created by default in `bin` with the case name and suite name appended.
+
+## Running the SCM with FV3GFS initial conditions
+model initial conditions are needed to initialize the land surface in order to run with an interactive land model
+1. Prepare model initial conditions.
+  * `cd to scm/etc/scripts/`
+  modify path to files you can read in extract_FV3GFS_column_ic.py, this is set up for C96
+2. run extract_FV3GFS_column_ic.py, it will create fv3_model_point.nc in ../../data/processed_case_input/
+  * `./extract_FV3GFS_column_ic.py`
+3. cd to bin directory
+  * `cd ../../bin/`
+4. Run the SCM with the fv3_model_point case and C96 namelist
+  * `./run_gmtb_scm.py -c fv3_model_point -n input_GFS_v15_C96.nml`

ccpp/physics_namelists/input_GFS_v15_C96.nml

@@ -0,0 +1,104 @@
+&gfs_physics_nml
+       fhzero         = 6.
+       ldiag3d        = .true.
+       fhcyc          = 24.
+       nst_anl        = .true.
+       use_ufo        = .true.
+       pre_rad        = .false.
+       ncld           = 5
+       imp_physics    = 11
+       pdfcld         = .false.
+       fhswr          = 3600.
+       fhlwr          = 3600.
+       ialb           = 1
+       iems           = 1
+       IAER           = 111
+       ico2           = 2
+       isubc_sw       = 2
+       isubc_lw       = 2
+       isol           = 2
+       lwhtr          = .true.
+       swhtr          = .true.
+       cnvgwd         = .true.
+       shal_cnv       = .true.
+       cal_pre        = .false.
+       redrag         = .true.
+       dspheat        = .true.
+       hybedmf        = .true.
+       satmedmf       = .false.
+       shinhong       = .false.
+       do_ysu         = .false.
+       lheatstrg      = .false.
+       lgfdlmprad     = .false.
+       effr_in        = .false.
+       random_clds    = .false.
+       trans_trac     = .false.
+       cnvcld         = .true.
+       imfshalcnv     = 2
+       imfdeepcnv     = 2
+       cdmbgwd        = 0.125,3.0
+       prslrd0        = 0.
+       ivegsrc        = 1
+       isot           = 1
+       debug          = .false.
+       oz_phys        = .false.
+       oz_phys_2015   = .true.
+       h2o_phys       = .true.
+       nstf_name      = 2,1,1,0,5
+       xkzminv        = 0.3
+       xkzm_m         = 1.0
+       xkzm_h         = 1.0
+       do_sppt        = .false.
+       do_shum        = .false.
+       do_skeb        = .false.
+       do_sfcperts    = .false.
+/
+
+&gfdl_cloud_microphysics_nml
+       sedi_transport = .true.
+       do_sedi_heat = .false.
+       rad_snow = .true.
+       rad_graupel = .true.
+       rad_rain = .true.
+       const_vi = .F.
+       const_vs = .F.
+       const_vg = .F.
+       const_vr = .F.
+       vi_max = 1.
+       vs_max = 2.
+       vg_max = 12.
+       vr_max = 12.
+       qi_lim = 1.
+       prog_ccn = .false.
+       do_qa = .false.
+       fast_sat_adj = .false.
+       tau_l2v = 225.
+       tau_v2l = 150.
+       tau_g2v = 900.
+       rthresh = 10.e-6  ! This is a key parameter for cloud water
+       dw_land  = 0.16
+       dw_ocean = 0.10
+       ql_gen = 1.0e-3
+       ql_mlt = 1.0e-3
+       qi0_crt = 8.0E-5
+       qs0_crt = 1.0e-3
+       tau_i2s = 1000.
+       c_psaci = 0.05
+       c_pgacs = 0.01
+       rh_inc = 0.30
+       rh_inr = 0.30
+       rh_ins = 0.30
+       ccn_l = 300.
+       ccn_o = 100.
+       c_paut = 0.5
+       c_cracw = 0.8
+       use_ppm = .false.
+       use_ccn = .true.
+       mono_prof = .true.
+       z_slope_liq  = .true.
+       z_slope_ice  = .true.
+       de_ice = .false.
+       fix_negative = .true.
+       icloud_f = 1
+       mp_time = 150.
+/

scm/etc/Cheyenne_setup_gnu.csh

@@ -10,6 +10,7 @@ module load gnu/8.3.0
 module load ncarcompilers/0.5.0
 module load mpt/2.19
 module load netcdf/4.6.3
+module load python/2.7.16
 
 echo "Setting CC/CXX/FC environment variables"
 setenv CC gcc

scm/etc/Cheyenne_setup_gnu.sh

@@ -10,6 +10,7 @@ module load gnu/8.3.0
 module load ncarcompilers/0.5.0
 module load mpt/2.19
 module load netcdf/4.6.3
+module load python/2.7.16
 
 echo "Setting CC/CXX/FC environment variables"
 export CC=gcc

scm/etc/Cheyenne_setup_intel.csh

@@ -10,6 +10,7 @@ module load intel/19.0.2
 module load ncarcompilers/0.5.0
 module load mpt/2.19
 module load netcdf/4.6.3
+module load python/2.7.16
 
 echo "Setting CC/CXX/FC environment variables"
 setenv CC icc

scm/etc/Cheyenne_setup_intel.sh

@@ -10,6 +10,7 @@ module load intel/19.0.2
 module load ncarcompilers/0.5.0
 module load mpt/2.19
 module load netcdf/4.6.3
+module load python/2.7.16
 
 echo "Setting CC/CXX/FC environment variables"
 export CC=icc

scm/etc/Cheyenne_setup_pgi.csh

@@ -9,6 +9,7 @@ module load ncarenv/1.3
 module load pgi/19.3
 module load ncarcompilers/0.5.0
 module load netcdf/4.6.3
+module load python/2.7.16
 
 echo "Setting CC/CXX/FC environment variables"
 setenv CC pgcc

scm/etc/Cheyenne_setup_pgi.sh

@@ -9,6 +9,7 @@ module load ncarenv/1.3
 module load pgi/19.3
 module load ncarcompilers/0.5.0
 module load netcdf/4.6.3
+module load python/2.7.16
 
 echo "Setting CC/CXX/FC environment variables"
 export CC=pgcc

scm/etc/case_config/fv3_model_point.nml

@@ -0,0 +1,26 @@
+$case_config
+model_name = 'FV3',
+n_columns = 1,
+case_name = 'fv3_model_point',
+dt = 900.0,
+time_scheme = 1,
+runtime = 86400,
+output_frequency = 900.0,
+n_levels = 64,
+output_file = 'output',
+case_data_dir = '../data/processed_case_input',
+vert_coord_data_dir = '../data/vert_coord_data',
+thermo_forcing_type = 0,
+model_ics=.true.,
+mom_forcing_type = 0,
+relax_time = 7200.0,
+sfc_flux_spec = .false.,
+sfc_type=1,
+reference_profile_choice = 2,
+C_RES = 96,
+year = 2014,
+month = 8,
+day = 1,
+hour = 0,
+$end
+

scm/etc/pjpegion_setup.sh

@@ -0,0 +1,6 @@
+#!/bin/bash
+
+echo "Setting environment variables for SCM-CCPP on MACOSX with clang/gfortran"
+
+export CFLAGS="-fmax-stack-var-size=30000 -I /Volumes/Cluster/opt/local/include/"
+export FFLAGS="-fmax-stack-var-size=30000 -fbounds-check -fcheck=all -gdwarf-4 -fvar-tracking-assignments -fbacktrace -fcheck=bounds -fdec"