32-bit physics with FV3_RAP

CMakeLists.txt

@@ -31,6 +31,12 @@ else()
   list(APPEND _fv3atm_defs_private NO_INLINE_POST)
 endif()
 
+if(CCPP_32BIT)
+  add_definitions(-DCCPP_32BIT)
+else()
+  remove_definitions(-DCCPP_32BIT)
+endif()
+
 if(NOT PARALLEL_NETCDF)
   list(APPEND _fv3atm_defs_private NO_PARALLEL_NETCDF)
 endif()

atmos_model.F90

@@ -137,14 +137,14 @@ module atmos_model_mod
      integer                       :: mlon, mlat
      integer                       :: iau_offset         ! iau running window length
      logical                       :: pe                 ! current pe.
-     real(kind=8),             pointer, dimension(:)     :: ak, bk
+     real(kind=GFS_kind_phys), pointer, dimension(:)     :: ak, bk
      real(kind=GFS_kind_phys), pointer, dimension(:,:)   :: lon_bnd  => null() ! local longitude axis grid box corners in radians.
      real(kind=GFS_kind_phys), pointer, dimension(:,:)   :: lat_bnd  => null() ! local latitude axis grid box corners in radians.
      real(kind=GFS_kind_phys), pointer, dimension(:,:)   :: lon      => null() ! local longitude axis grid box centers in radians.
      real(kind=GFS_kind_phys), pointer, dimension(:,:)   :: lat      => null() ! local latitude axis grid box centers in radians.
      real(kind=GFS_kind_phys), pointer, dimension(:,:)   :: dx, dy
-     real(kind=8),             pointer, dimension(:,:)   :: area
-     real(kind=8),             pointer, dimension(:,:,:) :: layer_hgt, level_hgt
+     real(kind=GFS_kind_phys), pointer, dimension(:,:)   :: area
+     real(kind=GFS_kind_phys), pointer, dimension(:,:,:) :: layer_hgt, level_hgt
      type(domain2d)                :: domain             ! domain decomposition
      type(domain2d)                :: domain_for_read    ! domain decomposition
      type(time_type)               :: Time               ! current time
@@ -467,9 +467,9 @@ subroutine atmos_timestep_diagnostics(Atmos)
             psum  = psum + adiff
             if(adiff>=maxabs) then
               maxabs=adiff
-              pmaxloc(2:3) = (/ ATM_block%index(nb)%ii(i), ATM_block%index(nb)%jj(i) /)
-              pmaxloc(4:7) = (/ pdiff, GFS_data(nb)%Statein%pgr(i), &
-                   GFS_data(nb)%Grid%xlat(i), GFS_data(nb)%Grid%xlon(i) /)
+              pmaxloc(2:3) = (/ dble(ATM_block%index(nb)%ii(i)), dble(ATM_block%index(nb)%jj(i)) /)
+              pmaxloc(4:7) = (/ dble(pdiff), dble(GFS_data(nb)%Statein%pgr(i)), &
+                   dble(GFS_data(nb)%Grid%xlat(i)), dble(GFS_data(nb)%Grid%xlon(i)) /)
             endif
           enddo
           pcount = pcount+count
@@ -2883,6 +2883,7 @@ subroutine setup_exportdata(rc)
             ! Instantaneous u wind (m/s) 10 m above ground
             case ('inst_zonal_wind_height10m')
               call block_data_copy(datar82d, GFS_data(nb)%coupling%u10mi_cpl, Atm_block, nb, rc=localrc)
+              !call block_data_copy(datar82d, GFS_data(nb)%coupling%u10mi_cpl, Atm_block, nb, rc=localrc)
             ! Instantaneous v wind (m/s) 10 m above ground
             case ('inst_merid_wind_height10m')
               call block_data_copy(datar82d, GFS_data(nb)%coupling%v10mi_cpl, Atm_block, nb, rc=localrc)

ccpp/CMakeLists.txt

@@ -53,22 +53,35 @@ endif()
 #------------------------------------------------------------------------------
 # Set flag for 32bit dynamics build
 if(32BIT)
-  message(STATUS "Compile CCPP slow physics with 64-bit precision, fast physics with 32-bit precision")
+  message(STATUS "Compile CCPP fast physics with 32-bit precision")
   add_definitions(-DOVERLOAD_R4)
-  if(CMAKE_Fortran_COMPILER_ID MATCHES "Intel")
-    set(CMAKE_Fortran_FLAGS_PHYSICS  "-real-size 64")
-  elseif(CMAKE_Fortran_COMPILER_ID MATCHES "GNU")
-    set(CMAKE_Fortran_FLAGS_PHYSICS  "-fdefault-real-8 -fdefault-double-8")
-  endif()
   set(CMAKE_Fortran_FLAGS_DYNAMICS "")
 else()
-  message(STATUS "Compile CCPP physics with 64-bit precision")
+  message(STATUS "Compile CCPP fast physics with 64-bit precision")
   remove_definitions(-DOVERLOAD_R8)
   remove_definitions(-DOVERLOAD_R4)
-  set(CMAKE_Fortran_FLAGS_PHYSICS  "")
-  set(CMAKE_Fortran_FLAGS_DYNAMICS "")
 endif()
 
+if(CCPP_32BIT)
+  message(STATUS "Compile CCPP slow physics with 32-bit precision")
+  add_definitions(-DSINGLE_PREC)
+  add_definitions(-DRTE_USE_SP)
+  if(CMAKE_Fortran_COMPILER_ID MATCHES "Intel")
+    set(CMAKE_Fortran_FLAGS_PHYSICS  "-real-size 32")
+  elseif(CMAKE_Fortran_COMPILER_ID MATCHES "GNU")
+    set(CMAKE_Fortran_FLAGS_PHYSICS  "-fno-default-real-8 -fdefault-double-8")
+  endif()
+else(CCPP_32BIT)
+  message(STATUS "Compile CCPP slow physics with 64-bit precision")
+  remove_definitions(-DSINGLE_PREC)
+  remove_definitions(-DRTE_USE_SP)
+  if(CMAKE_Fortran_COMPILER_ID MATCHES "Intel")
+    set(CMAKE_Fortran_FLAGS_PHYSICS  "-real-size 64")
+  elseif(CMAKE_Fortran_COMPILER_ID MATCHES "GNU")
+    set(CMAKE_Fortran_FLAGS_PHYSICS  "-fdefault-real-8 -fdefault-double-8")
+  endif()
+endif(CCPP_32BIT)
+
 #------------------------------------------------------------------------------
 # Add model-specific flags for C/C++/Fortran preprocessor
 add_definitions(-DMOIST_CAPPA -DUSE_COND -DNEMS_GSM)

ccpp/data/GFS_typedefs.F90

@@ -1,6 +1,6 @@
 module GFS_typedefs
 
-   use machine,                  only: kind_phys, kind_dbl_prec
+   use machine,                  only: kind_phys, kind_dbl_prec, kind_sngl_prec
    use physcons,                 only: con_cp, con_fvirt, con_g,                       &
                                        con_hvap, con_hfus, con_pi, con_rd, con_rv,     &
                                        con_t0c, con_cvap, con_cliq, con_eps, con_epsq, &
@@ -2895,6 +2895,8 @@ subroutine control_initialize (Model, nlunit, fn_nml, me, master, &
     logical :: exists
     real(kind=kind_phys) :: tem
     real(kind=kind_phys) :: rinc(5)
+    real(kind=kind_sngl_prec) :: rinc4(5)
+    real(kind=kind_dbl_prec) :: rinc8(5)
     real(kind=kind_phys) :: wrk(1)
     real(kind=kind_phys), parameter :: con_hr = 3600.
 
@@ -3385,7 +3387,11 @@ subroutine control_initialize (Model, nlunit, fn_nml, me, master, &
 !  max and min lon and lat for critical relative humidity
     integer :: max_lon=5000, max_lat=2000, min_lon=192, min_lat=94
     real(kind=kind_phys) :: rhcmax = 0.9999999               !< max critical rel. hum.
+#ifdef SINGLE_PREC
+    real(kind=kind_phys) :: huge   = 9.9692099683868690E30  !  NetCDF float FillValue
+#else
     real(kind=kind_phys) :: huge   = 9.9692099683868690E36  !  NetCDF float FillValue
+#endif
 
 
 !--- stochastic physics control parameters
@@ -3429,6 +3435,7 @@ subroutine control_initialize (Model, nlunit, fn_nml, me, master, &
 
     real(kind=kind_phys) :: radar_tten_limits(2) = (/ limit_unspecified, limit_unspecified /)
     integer :: itime
+    integer :: w3kindreal,w3kindint
 
 !--- END NAMELIST VARIABLES
 
@@ -4874,7 +4881,19 @@ subroutine control_initialize (Model, nlunit, fn_nml, me, master, &
     Model%cdec             = -9999.
     Model%clstp            = -9999
     rinc(1:5)              = 0
-    call w3difdat(jdat,idat,4,rinc)
+    call w3kind(w3kindreal,w3kindint)
+    if (w3kindreal == 8) then
+       rinc8(1:5) = 0
+       call w3difdat(jdat,idat,4,rinc8)
+       rinc = rinc8
+    else if (w3kindreal == 4) then
+       rinc4(1:5) = 0
+       call w3difdat(jdat,idat,4,rinc4)
+       rinc = rinc4
+    else
+       write(0,*)' FATAL ERROR: Invalid w3kindreal'
+       call abort
+    endif
     Model%phour            = rinc(4)/con_hr
     Model%fhour            = (rinc(4) + Model%dtp)/con_hr
     Model%zhour            = mod(Model%phour,Model%fhzero)