Ca develop (#96)

Updates to cellular automata stochastic physics. Added new regression test. Added possibility to run global and sub-grid CA, updates to coupling to convection.
ufs-community · May 5, 2020 · 6585dd5 · 6585dd5
1 parent d202a02
commit 6585dd5
Show file tree

Hide file tree

Showing 36 changed files with 51,256 additions and 49,716 deletions.
diff --git a/FV3 b/FV3
diff --git a/compsets/fv3.input b/compsets/fv3.input
@@ -266,6 +266,10 @@ fv3_defaults = {
  IAU_INC_FILES="''"
  IAU_DRYMASSFIXER='.false.'
 
+ DO_CA='.F.'
+ CA_SGS='.F.'
+ CA_GLOBAL='.F.'
+
  CPLFLX='.F.'
  CPLWAV='.F.'
  CPLWAV2ATM='.F.'

diff --git a/log/report-hera.intel-log/build_fv3.exe.log b/log/report-hera.intel-log/build_fv3.exe.log
diff --git a/log/report-hera.intel-log/build_fv3_32bit.exe.log b/log/report-hera.intel-log/build_fv3_32bit.exe.log
diff --git a/log/report-hera.intel-log/build_fv3_32bit_debug.exe.log b/log/report-hera.intel-log/build_fv3_32bit_debug.exe.log
diff --git a/log/report-hera.intel-log/build_fv3_appbuilder.exe.log b/log/report-hera.intel-log/build_fv3_appbuilder.exe.log
diff --git a/log/report-hera.intel-log/build_fv3_ww3.exe.log b/log/report-hera.intel-log/build_fv3_ww3.exe.log
diff --git a/log/report-hera.intel-log/rtreport.txt b/log/report-hera.intel-log/rtreport.txt
diff --git a/log/report-wcoss.cray-log/build_fv3.exe.log b/log/report-wcoss.cray-log/build_fv3.exe.log
diff --git a/log/report-wcoss.cray-log/build_fv3_32bit.exe.log b/log/report-wcoss.cray-log/build_fv3_32bit.exe.log
diff --git a/log/report-wcoss.cray-log/build_fv3_appbuilder.exe.log b/log/report-wcoss.cray-log/build_fv3_appbuilder.exe.log
diff --git a/log/report-wcoss.cray-log/build_fv3_ww3.exe.log b/log/report-wcoss.cray-log/build_fv3_ww3.exe.log
diff --git a/log/report-wcoss.cray-log/rtreport.txt b/log/report-wcoss.cray-log/rtreport.txt
diff --git a/log/report-wcoss_dell_p3-log/build_fv3.exe.log b/log/report-wcoss_dell_p3-log/build_fv3.exe.log
diff --git a/log/report-wcoss_dell_p3-log/build_fv3_32bit.exe.log b/log/report-wcoss_dell_p3-log/build_fv3_32bit.exe.log
diff --git a/log/report-wcoss_dell_p3-log/build_fv3_appbuilder.exe.log b/log/report-wcoss_dell_p3-log/build_fv3_appbuilder.exe.log
diff --git a/log/report-wcoss_dell_p3-log/build_fv3_ww3.exe.log b/log/report-wcoss_dell_p3-log/build_fv3_ww3.exe.log
diff --git a/log/report-wcoss_dell_p3-log/rtreport.txt b/log/report-wcoss_dell_p3-log/rtreport.txt
diff --git a/modulefiles/jet.intel/fv3 b/modulefiles/jet.intel/fv3
@@ -35,6 +35,8 @@ module load esmf/8.0.0
 ##
 ## load cmake
 ##
+module use /apps/modules/modulefiles
+module load cmake/3.16.1
 setenv CMAKE_C_COMPILER mpiicc
 setenv CMAKE_CXX_COMPILER mpiicpc
 setenv CMAKE_Fortran_COMPILER mpiifort

diff --git a/parm/ccpp_ca.nml.IN b/parm/ccpp_ca.nml.IN
@@ -0,0 +1,311 @@
+&amip_interp_nml
+ interp_oi_sst = .true.
+ use_ncep_sst = .true.
+ use_ncep_ice = .false.
+ no_anom_sst = .false.
+ data_set = 'reynolds_oi',
+ date_out_of_range = 'climo',
+/
+
+&atmos_model_nml
+ blocksize = 32
+ chksum_debug = .false.
+ dycore_only = .false.
+ fdiag = @[FDIAG]
+ ccpp_suite = '@[CCPP_SUITE]'
+/
+
+&diag_manager_nml
+ prepend_date = .F.
+/
+
+&fms_io_nml
+ checksum_required = .false.
+ max_files_r = 100,
+ max_files_w = 100,
+/
+
+&fms_nml
+ clock_grain = 'ROUTINE',
+ domains_stack_size = 3000000,
+ print_memory_usage = .false.
+/
+
+&fv_grid_nml
+ grid_file = 'INPUT/grid_spec.nc'
+/
+
+&fv_core_nml
+ layout = @[INPES],@[JNPES]
+ io_layout = 1,1
+ npx = @[NPX]
+ npy = @[NPY]
+ ntiles = 6,
+ npz = @[NPZ]
+ grid_type = -1
+ make_nh = @[MAKE_NH]
+ fv_debug = .F.
+ range_warn = .F.
+ reset_eta = .F.
+ n_sponge = 24
+ nudge_qv = .T.
+ tau = 5.
+ rf_cutoff = 7.5e2
+ d2_bg_k1 = 0.15
+ d2_bg_k2 = 0.02
+ kord_tm = -9
+ kord_mt = 9
+ kord_wz = 9
+ kord_tr = 9
+ hydrostatic = .F.
+ phys_hydrostatic = .F.
+ use_hydro_pressure = .F.
+ beta = 0.
+ a_imp = 1.
+ p_fac = 0.1
+ k_split = 2
+ n_split = 6
+ nwat = 2
+ na_init = @[NA_INIT]
+ d_ext = 0.0
+ dnats = 0
+ fv_sg_adj = 450
+ d2_bg = 0.
+ nord = 2
+ dddmp = 0.1
+ d4_bg = 0.12
+ vtdm4 = 0.02
+ delt_max = 0.002
+ ke_bg = 0.
+ do_vort_damp = .true.
+ external_ic = @[EXTERNAL_IC]
+ external_eta = .T.
+ gfs_phil = .false.
+ nggps_ic = @[NGGPS_IC]
+ mountain = @[MOUNTAIN]
+ ncep_ic = .F.
+ d_con = 1.
+ hord_mt = 6
+ hord_vt = 6
+ hord_tm = 6
+ hord_dp = 6
+ hord_tr = 8
+ adjust_dry_mass = .F.
+ consv_te = 1.
+ consv_am = .F.
+ fill = .T.
+ dwind_2d = .F.
+ print_freq = 6
+ warm_start = @[WARM_START]
+ no_dycore = .false.
+ z_tracer = .T.
+ read_increment = @[READ_INCREMENT]
+ res_latlon_dynamics = "fv3_increment.nc"
+/
+
+#&coupler_nml
+# months = 0
+# days = 1
+# hours = 0
+# dt_atmos = 225
+# dt_ocean = 225
+# current_date = 2016,10,03,00,0,0
+# calendar = 'julian'
+# memuse_verbose = .false.
+# atmos_nthreads = 1
+# use_hyper_thread = .false.
+# ncores_per_node = 24
+#/
+
+&external_ic_nml
+ filtered_terrain = .true.
+ levp = @[NPZP]
+ gfs_dwinds = .true.
+ checker_tr = .F.
+ nt_checker = 0
+/
+
+&gfs_physics_nml
+ fhzero = 6.
+ ldiag3d = .false.
+ fhcyc = 24.
+ nst_anl = .true.
+ use_ufo = .true.
+ pre_rad = .false.
+ ncld = 1
+ imp_physics = 99
+ 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 = .true.
+ redrag = .true.
+ dspheat = .true.
+ hybedmf = @[HYBEDMF]
+ satmedmf = @[SATMEDMF]
+ lheatstrg = @[LHEATSTRG]
+ random_clds = .true.
+ trans_trac = .false.
+ cnvcld = .true.
+ imfshalcnv = 2
+ imfdeepcnv = 2
+ isatmedmf = @[ISATMEDMF]
+ cdmbgwd = 3.5,0.25
+ prslrd0 = 0.
+ ivegsrc = 1
+ isot = 1
+ lsm = 1
+ iopt_dveg = 2
+ iopt_crs = 1
+ iopt_btr = 1
+ iopt_run = 1
+ iopt_sfc = 1
+ iopt_frz = 1
+ iopt_inf = 1
+ iopt_rad = 1
+ iopt_alb = 2
+ iopt_snf = 4
+ iopt_tbot = 2
+ iopt_stc = 1
+ debug = .false.
+ oz_phys = @[OZ_PHYS_OLD]
+ oz_phys_2015 = @[OZ_PHYS_NEW]
+ h2o_phys = @[H2O_PHYS]
+ nstf_name = @[NSTF_NAME]
+ cplflx = @[CPLFLX]
+ cplwav = @[CPLWAV]
+ iau_delthrs = 6
+ iaufhrs = 30
+ iau_inc_files = @[IAU_INC_FILES]
+ do_sppt = @[DO_SPPT]
+ do_shum = @[DO_SHUM]
+ do_skeb = @[DO_SKEB]
+ do_sfcperts = @[DO_SFCPERTS]
+ do_ca = @[DO_CA]
+ ca_sgs = @[CA_SGS]
+ ca_global = @[CA_GLOBAL]
+ nca = 1
+ ncells = 5
+ nlives = 24
+ nseed = 100
+ nfracseed = 0.5
+ nthresh = 8
+ ca_trigger = .True.
+ ca_entr = .False.
+ ca_closure = .False.
+ nca_g = 1
+ ncells_g = 1
+ nlives_g = 100
+ nseed_g = 100
+ ca_smooth = .True.
+ nspinup = 100
+ iseed_ca = 12345
+ nsmooth = 100
+ ca_amplitude = 0.35
+/
+
+&interpolator_nml
+ interp_method = 'conserve_great_circle'
+/
+
+&namsfc
+ FNGLAC = "global_glacier.2x2.grb",
+ FNMXIC = "global_maxice.2x2.grb",
+ FNTSFC = "RTGSST.1982.2012.monthly.clim.grb",
+ FNSNOC = "global_snoclim.1.875.grb",
+ FNZORC = "igbp",
+ FNALBC = @[FNALBC]
+ FNALBC2 = "global_albedo4.1x1.grb",
+ FNAISC = "CFSR.SEAICE.1982.2012.monthly.clim.grb",
+ FNTG3C = "global_tg3clim.2.6x1.5.grb",
+ FNVEGC = "global_vegfrac.0.144.decpercent.grb",
+ FNVETC = @[FNVETC]
+ FNSOTC = @[FNSOTC]
+ FNSMCC = @[FNSMCC]
+ FNMSKH = "seaice_newland.grb",
+ FNTSFA = "",
+ FNACNA = "",
+ FNSNOA = "",
+ FNVMNC = "global_shdmin.0.144x0.144.grb",
+ FNVMXC = "global_shdmax.0.144x0.144.grb",
+ FNSLPC = "global_slope.1x1.grb",
+ FNABSC = @[FNABSC]
+ LDEBUG =.false.,
+ FSMCL(2) = 99999
+ FSMCL(3) = 99999
+ FSMCL(4) = 99999
+ FTSFS = 90
+ FAISS = 99999
+ FSNOL = 99999
+ FSICL = 99999
+ FTSFL = 99999,
+ FAISL = 99999,
+ FVETL = 99999,
+ FSOTL = 99999,
+ FvmnL = 99999,
+ FvmxL = 99999,
+ FSLPL = 99999,
+ FABSL = 99999,
+ FSNOS = 99999,
+ FSICS = 99999,
+/
+&nam_stochy
+ lon_s=768,
+ lat_s=384,
+ ntrunc=382,
+ SKEBNORM=1,
+ SKEB_NPASS=30,
+ SKEB_VDOF=5,
+ SKEB=@[SKEB],
+ SKEB_TAU=2.16E4,
+ SKEB_LSCALE=1000.E3,
+ SHUM=@[SHUM],
+ SHUM_TAU=21600,
+ SHUM_LSCALE=500000,
+ SPPT=@[SPPT],
+ SPPT_TAU=21600,
+ SPPT_LSCALE=500000,
+ SPPT_LOGIT=.TRUE.,
+ SPPT_SFCLIMIT=.TRUE.,
+ ISEED_SHUM=1,
+ ISEED_SKEB=2,
+ ISEED_SPPT=3,
+/
+&nam_sfcperts
+ NSFCPERT=6,
+ PERTZ0=-999.,
+ PERTSHC=-999.,
+ PERTZT=-999.,
+ PERTLAI=-999.,
+ PERTVEGF=-999.,
+ PERTALB=-999.,
+ SFC_TAU=21600,
+ SFC_LSCALE=500000,
+ ISEED_SFC=0,
+ SPPT_LAND=.FALSE.,
+/
+
+&cires_ugwp_nml
+ knob_ugwp_solver = 2
+ knob_ugwp_source = 1,1,0,0
+ knob_ugwp_wvspec = 1,25,25,25
+ knob_ugwp_azdir = 2,4,4,4
+ knob_ugwp_stoch = 0,0,0,0
+ knob_ugwp_effac = 1,1,1,1
+ knob_ugwp_doaxyz = 1
+ knob_ugwp_doheat = 1
+ knob_ugwp_dokdis = 1
+ knob_ugwp_ndx4lh = 1
+ knob_ugwp_version = 0
+ launch_level = 25
+/
diff --git a/parm/ccpp_control.nml.IN b/parm/ccpp_control.nml.IN
@@ -193,19 +193,27 @@
  do_shum = @[DO_SHUM]
  do_skeb = @[DO_SKEB]
  do_sfcperts = @[DO_SFCPERTS]
+ do_ca = .false.
+ ca_sgs = .false.
  nca = 1
  ncells = 5
- nfracseed = 0.5
  nlives = 30
- nthresh = 0.5
- nseed = 100000
- do_ca = .false.
+ nseed = 1000000
+ nfracseed = 0.5
+ nthresh = 0.
+ ca_trigger = .true. 
+ ca_entr = .false.
+ ca_closure = .false.
  ca_global = .false.
- ca_sgs = .true.
+ nca_g = 1
+ ncells_g = 1
+ nlives_g = 100
+ nseed_g = 100
  ca_smooth = .false.
- nspinup = 1000
- iseed_ca = 0
- isppt_deep = .false.
+ nspinup = 100
+ iseed_ca = 1
+ nsmooth = 1
+ ca_amplitude = 0.5
  do_RRTMGP = @[DO_RRTMGP]
  active_gases = 'h2o_co2_o3_n2o_ch4_o2'
  ngases = 6