clarified unit for stmp in nws=4; thx Saeed!

schism-dev · Jan 24, 2025 · aff3de3 · aff3de3
1 parent cec6f0d
commit aff3de3
Show file tree

Hide file tree

Showing 4 changed files with 7 additions and 4 deletions.
diff --git a/sample_inputs/atmos.nc b/sample_inputs/atmos.nc
diff --git a/src/Hydro/misc_subs.F90 b/src/Hydro/misc_subs.F90
@@ -314,7 +314,7 @@ subroutine other_hot_init(time)
           if(j/=NF90_NOERR) call parallel_abort('MISC: atmos.nc prmsl')
           j=nf90_get_var(ncid_atmos,mm,rwild(3,:),(/1,ninv+1/),(/np_global,1/))
           if(j/=NF90_NOERR) call parallel_abort('MISC: atmos.nc prmsl(2)')
-          j=nf90_inq_varid(ncid_atmos, "stmp",mm)
+          j=nf90_inq_varid(ncid_atmos, "stmp_in_centigrade",mm)
           if(j/=NF90_NOERR) call parallel_abort('MISC: atmos.nc stmp')
           j=nf90_get_var(ncid_atmos,mm,rwild(4,:),(/1,ninv+1/),(/np_global,1/))
           if(j/=NF90_NOERR) call parallel_abort('MISC: atmos.nc stmp(2)')
@@ -375,7 +375,7 @@ subroutine other_hot_init(time)
           j=nf90_inq_varid(ncid_atmos, "prmsl",mm)
           j=nf90_get_var(ncid_atmos,mm,rwild(3,:),(/1,ninv+2/),(/np_global,1/))
           if(j/=NF90_NOERR) call parallel_abort('MISC: atmos.nc prmsl(3)')
-          j=nf90_inq_varid(ncid_atmos, "stmp",mm)
+          j=nf90_inq_varid(ncid_atmos, "stmp_in_centigrade",mm)
           j=nf90_get_var(ncid_atmos,mm,rwild(4,:),(/1,ninv+1/),(/np_global,1/))
           if(j/=NF90_NOERR) call parallel_abort('MISC: atmos.nc stmp(3)')
           j=nf90_inq_varid(ncid_atmos, "spfh",mm)

diff --git a/src/Hydro/schism_step.F90 b/src/Hydro/schism_step.F90
@@ -599,7 +599,8 @@ subroutine schism_step(it)
             if(j/=NF90_NOERR) call parallel_abort('STEP: atmos.nc prmsl')
             j=nf90_get_var(ncid_atmos,mm,rwild6(3,:),(/1,itmp2/),(/np_global,1/))
             if(j/=NF90_NOERR) call parallel_abort('STEP: atmos.nc prmsl(2)')
-            j=nf90_inq_varid(ncid_atmos, "stmp",mm)
+            !air T in centigrade not Kelvin
+            j=nf90_inq_varid(ncid_atmos, "stmp_in_centigrade",mm)
             if(j/=NF90_NOERR) call parallel_abort('STEP: atmos.nc stmp')
             j=nf90_get_var(ncid_atmos,mm,rwild6(4,:),(/1,itmp2/),(/np_global,1/))
             if(j/=NF90_NOERR) call parallel_abort('STEP: atmos.nc stmp(2)')
@@ -658,7 +659,7 @@ subroutine schism_step(it)
           call exchange_p2d(windx2)
           call exchange_p2d(windy2)
           call exchange_p2d(pr2)
-          call exchange_p2d(airt2)
+          call exchange_p2d(airt2) !centigrade
           call exchange_p2d(shum2)
           call exchange_p2d(srad)
           call exchange_p2d(hradd)

diff --git a/src/Hydro/sflux_9c.F90 b/src/Hydro/sflux_9c.F90
@@ -67,6 +67,7 @@
 !	get_rad          (called from within surf_fluxes)
 !	get_precip_flux  (called from within surf_fluxes - IF ENABLED)
 !	surf_fluxes      (called from within SCHISM)
+!	surf_fluxes2      (alternative form called from within SCHISM)
 !
 ! In addition, there are a number of secondary routines and functions
 ! that are called by those listed above. For a complete list see below
@@ -266,6 +267,7 @@
 !
 ! List of all routines in this file:
 !   surf_fluxes
+!   surf_fluxes2
 !   turb_fluxes:  Calculate bulk aerodynamic surface fluxes over water using method of
 !                 Zeng et al or Fairall (computes the bulk parameterization of surface wind stress and
 !                 surface net heat fluxes using method of Fairall et al.