(Remove interstitial variables for land and ice emissivity and update the land and ice emissivity in the routine setemis) Sm sept21 pr

physics/GFS_PBL_generic.F90

@@ -336,7 +336,7 @@ subroutine GFS_PBL_generic_post_run (im, levs, nvdiff, ntrac,
         dqsfc_diag, dusfci_diag, dvsfci_diag, dtsfci_diag, dqsfci_diag,                                                        &
         rd, cp, fvirt, hvap, t1, q1, prsl, hflx, ushfsfci, oceanfrac, kdt, dusfc_cice, dvsfc_cice,                             &
         dtsfc_cice, dqsfc_cice, wet, dry, icy, wind, stress_wat, hflx_wat, evap_wat, ugrs1, vgrs1, hffac, &
-        ugrs, vgrs, tgrs, qgrs, save_u, save_v, save_t, save_q, errmsg, errflg)
+        ugrs, vgrs, tgrs, qgrs, save_u, save_v, save_t, save_q, huge, errmsg, errflg)
 
       use machine,                only : kind_phys
       use GFS_PBL_generic_common, only : set_aerosol_tracer_index
@@ -357,7 +357,7 @@ subroutine GFS_PBL_generic_post_run (im, levs, nvdiff, ntrac,
       real(kind=kind_phys), dimension(:,:, :), intent(in) :: save_q
 
       real(kind=kind_phys), intent(in) :: dtf
-      real(kind=kind_phys), intent(in) :: rd, cp, fvirt, hvap
+      real(kind=kind_phys), intent(in) :: rd, cp, fvirt, hvap, huge
       real(kind=kind_phys), dimension(:), intent(in) :: t1, q1, hflx, oceanfrac
       real(kind=kind_phys), dimension(:,:), intent(in) :: prsl
       real(kind=kind_phys), dimension(:), intent(in) :: dusfc_cice, dvsfc_cice, dtsfc_cice, dqsfc_cice, &
@@ -392,7 +392,7 @@ subroutine GFS_PBL_generic_post_run (im, levs, nvdiff, ntrac,
       integer, intent(out) :: errflg
 
       real(kind=kind_phys), parameter :: zero  = 0.0_kp, one = 1.0_kp
-      real(kind=kind_phys), parameter :: huge  = 9.9692099683868690E36 ! NetCDF float FillValue, same as in GFS_typedefs.F90
+!     real(kind=kind_phys), parameter :: huge  = 9.9692099683868690E36 ! NetCDF float FillValue, same as in GFS_typedefs.F90
       real(kind=kind_phys), parameter :: qmin  = 1.0e-8_kp
       integer :: i, k, kk, k1, n
       real(kind=kind_phys) :: tem, rho

physics/GFS_PBL_generic.meta

@@ -1383,6 +1383,15 @@
   kind = kind_phys
   intent = in
   optional = F
+[huge]
+  standard_name = netcdf_float_fillvalue
+  long_name = definition of NetCDF float FillValue
+  units = none
+  dimensions = ()
+  type = real
+  kind = kind_phys
+  intent = in
+  optional = F
 [errmsg]
   standard_name = ccpp_error_message
   long_name = error message for error handling in CCPP

physics/GFS_debug.F90

@@ -1307,7 +1307,7 @@ subroutine GFS_interstitialtoscreen_run (Model, Statein, Stateout, Sfcprop, Coup
                      call print_var(mpirank, omprank, blkno, Grid%xlat_d, Grid%xlon_d, 'Interstitial%scmpsw%visdf        ', Interstitial%scmpsw%visdf            )
                      call print_var(mpirank, omprank, blkno, Grid%xlat_d, Grid%xlon_d, 'Interstitial%semis_ice           ', Interstitial%semis_ice               )
                      call print_var(mpirank, omprank, blkno, Grid%xlat_d, Grid%xlon_d, 'Interstitial%semis_land          ', Interstitial%semis_land              )
-                     call print_var(mpirank, omprank, blkno, Grid%xlat_d, Grid%xlon_d, 'Interstitial%semis_water         ', Interstitial%semis_water             )
+!                    call print_var(mpirank, omprank, blkno, Grid%xlat_d, Grid%xlon_d, 'Interstitial%semis_water         ', Interstitial%semis_water             )
                      call print_var(mpirank, omprank, blkno, Grid%xlat_d, Grid%xlon_d, 'Interstitial%sfcalb              ', Interstitial%sfcalb                  )
                      call print_var(mpirank, omprank, blkno, Grid%xlat_d, Grid%xlon_d, 'Interstitial%sigma               ', Interstitial%sigma                   )
                      call print_var(mpirank, omprank, blkno, Grid%xlat_d, Grid%xlon_d, 'Interstitial%sigmaf              ', Interstitial%sigmaf                  )

physics/GFS_phys_time_vary.fv3.meta

@@ -1744,8 +1744,8 @@
   intent = inout
   optional = F
 [tisfc]
-  standard_name = sea_ice_temperature
-  long_name = sea ice surface skin temperature
+  standard_name = surface_skin_temperature_over_ice
+  long_name = surface skin temperature over ice
   units = K
   dimensions = (horizontal_dimension)
   type = real

physics/GFS_radiation_surface.F90

@@ -58,11 +58,11 @@ end subroutine GFS_radiation_surface_init
       subroutine GFS_radiation_surface_run (                            &
         im, frac_grid, lslwr, lsswr, lsm, lsm_noahmp, lsm_ruc,          &
         xlat, xlon, slmsk, lndp_type, n_var_lndp, sfc_alb_pert,         &
-        lndp_var_list, lndp_prt_list, landfrac, snowd, sncovr,          &
+        lndp_var_list, lndp_prt_list, landfrac, snodl, snodi, sncovr,   &
         sncovr_ice, fice, zorl, hprime, tsfg, tsfa, tisfc, coszen,      &
-        min_seaice, min_lakeice, lakefrac,                              &
+        cplice, min_seaice, min_lakeice, lakefrac, use_flake,           &
         alvsf, alnsf, alvwf, alnwf, facsf, facwf,                       &
-        semis_lnd, semis_ice, snoalb, use_cice_alb,                     &
+        semis_lnd, semis_ice, semis_wat, snoalb, use_cice_alb,          &
         albdvis_lnd, albdnir_lnd, albivis_lnd, albinir_lnd,             &
         albdvis_ice, albdnir_ice, albivis_ice, albinir_ice,             &
         semisbase, semis, sfcalb, sfc_alb_dif, errmsg, errflg)
@@ -73,28 +73,31 @@ subroutine GFS_radiation_surface_run (                            &
 
       implicit none
 
-      integer,              intent(in) :: im
-      logical,              intent(in) :: frac_grid, lslwr, lsswr, use_cice_alb
-      integer,              intent(in) :: lsm, lsm_noahmp, lsm_ruc, lndp_type, n_var_lndp
-      real(kind=kind_phys), intent(in) :: min_seaice, min_lakeice
+      integer,               intent(in) :: im
+      logical,               intent(in) :: frac_grid, lslwr, lsswr, use_cice_alb, cplice
+      integer,               intent(in) :: lsm, lsm_noahmp, lsm_ruc, lndp_type, n_var_lndp
+      real(kind=kind_phys),  intent(in) :: min_seaice, min_lakeice
+      logical, dimension(:), intent(in) :: use_flake
 
       real(kind=kind_phys), dimension(:),   intent(in)  :: xlat, xlon, slmsk,           &
                                                            sfc_alb_pert, lndp_prt_list, &
                                                            landfrac, lakefrac,          &
-                                                           snowd, sncovr,               &
+                                                           snodl, snodi, sncovr,        &
                                                            sncovr_ice, fice, zorl,      &
                                                            hprime, tsfg, tsfa, tisfc,   &
                                                            coszen, alvsf, alnsf, alvwf, &
-                                                           alnwf, facsf, facwf,         &
-                                                           semis_lnd, semis_ice, snoalb
+                                                           alnwf, facsf, facwf, snoalb
       character(len=3)    , dimension(:),   intent(in)  :: lndp_var_list
-      real(kind=kind_phys), dimension(:),   intent(inout)  :: albdvis_lnd, albdnir_lnd, &
-                                                           albivis_lnd, albinir_lnd
       real(kind=kind_phys), dimension(:),   intent(in)  :: albdvis_ice, albdnir_ice,    &
                                                            albivis_ice, albinir_ice
+
+      real(kind=kind_phys), dimension(:),   intent(inout) :: albdvis_lnd, albdnir_lnd,  &
+                                                             albivis_lnd, albinir_lnd,  &
+                                                             semis_lnd,   semis_ice, semis_wat
       real(kind=kind_phys), dimension(:),   intent(inout) :: semisbase, semis
       real(kind=kind_phys), dimension(:,:), intent(inout) :: sfcalb
       real(kind=kind_phys), dimension(:),   intent(inout) :: sfc_alb_dif
+
       character(len=*),                     intent(out) :: errmsg
       integer,                              intent(out) :: errflg
 
@@ -158,13 +161,13 @@ subroutine GFS_radiation_surface_run (                            &
       if (lslwr) then
 !>  - Call module_radiation_surface::setemis(),to set up surface
 !! emissivity for LW radiation.
-        call setemis (lsm, lsm_noahmp, lsm_ruc,                    &
-                      frac_grid,             xlon, xlat, slmsk,    &
-!                     frac_grid, min_seaice, xlon, xlat, slmsk,    &
-                      snowd, sncovr, sncovr_ice, zorl, tsfg, tsfa, &
-                      hprime, semis_lnd, semis_ice, im,            &
-                      fracl, fraco, fraci, icy,                    & !  ---  inputs
-                      semisbase, semis)                              !  ---  outputs
+        call setemis (lsm, lsm_noahmp, lsm_ruc, frac_grid, cplice,  &
+                      use_flake, lakefrac, xlon, xlat, slmsk,       &
+!                     frac_grid, min_seaice, xlon, xlat, slmsk,     &
+                      snodl, snodi, sncovr, sncovr_ice, zorl, tsfg, &
+                      tsfa, hprime, semis_lnd, semis_ice, semis_wat,&
+                      im, fracl, fraco, fraci, icy,                 & !  ---  inputs
+                      semisbase, semis)                               !  ---  outputs
       endif
 
       if (lsswr) then
@@ -181,7 +184,7 @@ subroutine GFS_radiation_surface_run (                            &
 !>  - Call module_radiation_surface::setalb(),to set up surface
 !! albedor for SW radiation.
 
-        call setalb (slmsk, lsm, lsm_noahmp, lsm_ruc, use_cice_alb, snowd, sncovr, sncovr_ice, &
+        call setalb (slmsk, lsm, lsm_noahmp, lsm_ruc, use_cice_alb, snodi, sncovr, sncovr_ice, &
                      snoalb, zorl, coszen, tsfg, tsfa, hprime, frac_grid, lakefrac,            &
 !                    snoalb, zorl, coszen, tsfg, tsfa, hprime, frac_grid, min_seaice,          &
                      alvsf, alnsf, alvwf, alnwf, facsf, facwf, fice, tisfc,                    &

physics/GFS_radiation_surface.meta

@@ -197,9 +197,18 @@
   kind = kind_phys
   intent = in
   optional = F
-[snowd]
-  standard_name = lwe_surface_snow
-  long_name = water equivalent snow depth
+[snodl]
+  standard_name = surface_snow_thickness_water_equivalent_over_land
+  long_name = water equivalent snow depth over land
+  units = mm
+  dimensions = (horizontal_loop_extent)
+  type = real
+  kind = kind_phys
+  intent = in
+  optional = F
+[snodi]
+  standard_name = surface_snow_thickness_water_equivalent_over_ice
+  long_name = water equivalent snow depth over ice
   units = mm
   dimensions = (horizontal_loop_extent)
   type = real
@@ -270,8 +279,8 @@
   intent = in
   optional = F
 [tisfc]
-  standard_name = sea_ice_temperature
-  long_name = sea ice surface skin temperature
+  standard_name = surface_skin_temperature_over_ice
+  long_name = surface_skin_temperature_over_ice
   units = K
   dimensions = (horizontal_loop_extent)
   type = real
@@ -287,6 +296,14 @@
   kind = kind_phys
   intent = in
   optional = F
+[cplice]
+  standard_name = flag_for_sea_ice_coupling
+  long_name = flag controlling cplice collection (default on)
+  units = flag
+  dimensions = ()
+  type = logical
+  intent = in
+  optional = F
 [min_seaice]
   standard_name = min_sea_ice_area_fraction
   long_name = minimum sea ice value
@@ -314,6 +331,14 @@
   kind = kind_phys
   intent = in
   optional = F
+[use_flake]
+  standard_name = flag_for_using_flake
+  long_name = flag indicating lake points using flake model
+  units = flag
+  dimensions = (horizontal_loop_extent)
+  type = logical
+  intent = inout
+  optional = F
 [alvsf]
   standard_name = vis_albedo_strong_cosz
   long_name = mean vis albedo with strong cosz dependency
@@ -386,6 +411,15 @@
   kind = kind_phys
   intent = in
   optional = F
+[semis_wat]
+  standard_name = surface_longwave_emissivity_over_water
+  long_name = surface lw emissivity in fraction over water
+  units = frac
+  dimensions = (horizontal_loop_extent)
+  type = real
+  kind = kind_phys
+  intent = inout
+  optional = F
 [snoalb]
   standard_name = upper_bound_of_max_albedo_assuming_deep_snow
   long_name = maximum snow albedo