Merge pull request #1 from climbfuji/first_n_radiation_calls_on_physi…

…cs_timestep_for_spinup

Add capability to call radiation on physics timestep for first N time steps

gfsphysics/GFS_layer/GFS_typedefs.F90

@@ -590,6 +590,9 @@ module GFS_typedefs
     real(kind=kind_phys) :: fhlwr           !< frequency for longwave radiation (secs)
     integer              :: nsswr           !< integer trigger for shortwave radiation
     integer              :: nslwr           !< integer trigger for longwave  radiation
+#ifdef CCPP
+    integer              :: nhfrad          !< number of timesteps for which to call radiation on physics timestep (coldstarts)
+#endif
     integer              :: levr            !< number of vertical levels for radiation calculations
 #ifdef CCPP
     integer              :: levrp1          !< number of vertical levels for radiation calculations plus one
@@ -2703,6 +2706,9 @@ subroutine control_initialize (Model, nlunit, fn_nml, me, master, &
 !--- radiation parameters
     real(kind=kind_phys) :: fhswr          = 3600.           !< frequency for shortwave radiation (secs)
     real(kind=kind_phys) :: fhlwr          = 3600.           !< frequency for longwave radiation (secs)
+#ifdef CCPP
+    integer              :: nhfrad         = 0               !< number of timesteps for which to call radiation on physics timestep (coldstarts)
+#endif
     integer              :: levr           = -99             !< number of vertical levels for radiation calculations
     integer              :: nfxr           = 39+6            !< second dimension of input/output array fluxr   
     logical              :: aero_in        = .false.         !< flag for initializing aero data 
@@ -3072,6 +3078,9 @@ subroutine control_initialize (Model, nlunit, fn_nml, me, master, &
                                fhswr, fhlwr, levr, nfxr, aero_in, iflip, isol, ico2, ialb,  &
                                isot, iems, iaer, icliq_sw, iovr_sw, iovr_lw, ictm, isubc_sw,&
                                isubc_lw, crick_proof, ccnorm, lwhtr, swhtr,                 &
+#ifdef CCPP
+                               nhfrad,                                                      &
+#endif
                           ! IN CCN forcing
                                iccn,                                                        &
                           !--- microphysical parameterizations
@@ -3277,6 +3286,17 @@ subroutine control_initialize (Model, nlunit, fn_nml, me, master, &
     Model%fhlwr            = fhlwr
     Model%nsswr            = nint(fhswr/Model%dtp)
     Model%nslwr            = nint(fhlwr/Model%dtp)
+#ifdef CCPP
+    if (restart) then
+      Model%nhfrad         = 0
+      if (Model%me == Model%master .and. nhfrad>0) &
+        write(*,'(a)') 'Disable high-frequency radiation calls for restart run'
+    else
+      Model%nhfrad         = nhfrad
+      if (Model%me == Model%master .and. nhfrad>0) &
+        write(*,'(a,i0)') 'Number of high-frequency radiation calls for coldstart run: ', nhfrad
+    endif
+#endif
     if (levr < 0) then
       Model%levr           = levs
     else
@@ -4384,6 +4404,9 @@ subroutine control_print(Model)
       print *, ' fhlwr             : ', Model%fhlwr
       print *, ' nsswr             : ', Model%nsswr
       print *, ' nslwr             : ', Model%nslwr
+#ifdef CCPP
+      print *, ' nhfrad            : ', Model%nhfrad
+#endif
       print *, ' levr              : ', Model%levr
       print *, ' nfxr              : ', Model%nfxr
       print *, ' aero_in           : ', Model%aero_in

gfsphysics/GFS_layer/GFS_typedefs.meta

@@ -2025,6 +2025,12 @@
   dimensions = ()
   type = real
   kind = kind_phys
+[nhfrad]
+  standard_name = number_of_timesteps_for_radiation_calls_on_physics_timestep
+  long_name = number of timesteps for radiation calls on physics timestep (coldstarts only)
+  units = count
+  dimensions = ()
+  type = integer
 [levr]
   standard_name = number_of_vertical_layers_for_radiation_calculations
   long_name = number of vertical levels for radiation calculations