Skip to content
New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

Sign up for GitHub

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

Jump to bottom

UFS-dev PR#126 #1054

Merged
merged 10 commits into from
Mar 8, 2024
Merged

UFS-dev PR#126 #1054

Show file tree
Hide file tree
Changes from all commits
Commits
Show all changes
10 commits
Select commit Hold shift + click to select a range
9d59cca
add kind_phys where missing
DeniseWorthen Nov 6, 2023
737fced
Merge branch 'ufs-community:ufs/dev' into feature/fix_nstkind
DeniseWorthen Nov 10, 2023
6397387
clean up type mix-matches
DeniseWorthen Nov 10, 2023
3714e76
add parameter zero and clean up nst_parameters
DeniseWorthen Nov 13, 2023
ade9c10
Merge pull request #131 from DeniseWorthen/bugfix/grvfunction
dustinswales Nov 20, 2023
b24beb1
Merge branch 'ufs/dev' into feature/fix_nstkind
DeniseWorthen Nov 25, 2023
f8e1601
rename and update files
DeniseWorthen Nov 26, 2023
9b8848e
Merge branch 'feature/fix_nstkind' into ufs-dev-PR126
grantfirl Feb 6, 2024
b4a3436
Merge branch 'ufs-dev-PR120' into ufs-dev-PR126
grantfirl Feb 9, 2024
e136fb2
Merge branch 'main' into ufs-dev-PR126
grantfirl Feb 23, 2024
File filter

Filter by extension

Filter by extension
Conversations
Failed to load comments.
Loading
Jump to
Jump to file
Failed to load files.
Loading
Diff view
Diff view
1,819 changes: 912 additions & 907 deletions physics/module_nst_model.f90
View file
Open in desktop

Large diffs are not rendered by default.

133 changes: 73 additions & 60 deletions physics/module_nst_parameters.f90
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -9,74 +9,87 @@
!! history:
!! 20210305: X.Li, reduce z_w_max from 30 m to 20 m
module module_nst_parameters

use machine, only : kind_phys
!
! air constants and coefficients from the atmospehric model
use physcons, only: &
eps => con_eps &
,cp_a => con_cp & !< spec heat air @p (j/kg/k)
, epsm1 => con_epsm1 &
, hvap => con_hvap & !< lat heat h2o cond (j/kg)
,sigma_r => con_sbc & !< stefan-boltzmann (w/m2/k4)
,grav => con_g & !< acceleration due to gravity (kg/m/s^2)
,omega => con_omega & !< ang vel of earth (1/s)
,rvrdm1 => con_fvirt &
,rd => con_rd &
,rocp => con_rocp & !< r/cp
,pi => con_pi
use physcons, only: &
eps => con_eps & !< con_rd/con_rv (nd)
,cp_a => con_cp & !< spec heat air @p (j/kg/k)
,epsm1 => con_epsm1 & !< eps - 1 (nd)
,hvap => con_hvap & !< lat heat h2o cond (j/kg)
,sigma_r => con_sbc & !< stefan-boltzmann (w/m2/k4)
,grav => con_g & !< acceleration due to gravity (kg/m/s^2)
,omega => con_omega & !< ang vel of earth (1/s)
,rvrdm1 => con_fvirt & !< con_rv/con_rd-1. (nd)
,rd => con_rd & !< gas constant air (j/kg/k)
,rocp => con_rocp & !< r/cp
,pi => con_pi

implicit none

private

public :: sigma_r
public :: zero, one, half
public :: niter_conv, niter_z_w, niter_sfs
public :: z_w_max, z_w_min, z_w_ini, z_c_max, z_c_ini, eps_z_w, eps_conv, eps_sfs
public :: ri_c, ri_g, omg_m, omg_sh, tc_w, visw, cp_w, t0k, ustar_a_min, delz, exp_const
public :: rad2deg, const_rot, tw_max, sst_max, solar_time_6am, tau_min, wd_max

real(kind_phys), parameter :: zero = 0.0_kind_phys, one = 1.0_kind_phys, half = 0.5_kind_phys
!
! note: take timestep from here later
public
integer :: &
niter_conv = 5, &
niter_z_w = 5, &
niter_sfs = 5
real (kind=kind_phys), parameter :: &
!
! general constants
sec_in_day=86400. &
,sec_in_hour=3600. &
,solar_time_6am=21600.0 &
,const_rot=0.000073 & !< constant to calculate corioli force
,ri_c=0.65 &
,ri_g=0.25 &
,eps_z_w=0.01 & !< criteria to finish iterations for z_w
,eps_conv=0.01 & !< criteria to finish iterations for d_conv
,eps_sfs=0.01 & !< criteria to finish iterations for d_sfs
,z_w_max=20.0 & !< max warm layer thickness
,z_w_min=0.2 & !< min warm layer thickness
,z_w_ini=0.2 & !< initial warm layer thickness in dtl_onset
,z_c_max=0.01 & !< maximum of sub-layer thickness (m)
,z_c_ini=0.001 & !< initial value of z_c
,ustar_a_min=0.031 & !< minimum of friction wind speed (m/s): 0.031 ~ 1m/s at 10 m hight
,tau_min=0.005 & !< minimum of wind stress for dtm
,exp_const=9.5 & !< coefficient in exponet profile
,delz=0.1 & !< vertical increment for integral calculation (m)
,von=0.4 & !< von karman's "constant"
,t0k=273.16 & !< celsius to kelvin
,gray=0.97 &
,sst_max=308.16 &
,tw_max=5.0 &
,wd_max=2.0 &
,omg_m =1.0 & !< trace factor to apply salinity effect
,omg_rot = 1.0 & !< trace factor to apply rotation effect
,omg_sh = 1.0 & !< trace factor to apply sensible heat due to rainfall effect
,visw=1.e-6 & !< m2/s kinematic viscosity water
,novalue=0 &
,smallnumber=1.e-6 &
,timestep_oc=sec_in_day/8. & !< time step in the ocean model (3 hours)
,radian=2.*pi/180. &
,rad2deg=180./pi &
,cp_w=4000. & !< specific heat water (j/kg/k )
,rho0_w=1022.0 & !< density water (kg/m3 ) (or 1024.438)
,vis_w=1.e-6 & !< kinematic viscosity water (m2/s )
,tc_w=0.6 & !< thermal conductivity water (w/m/k )
,capa_w =3950.0 & !< heat capacity of sea water !
,thref =1.0e-3 !< reference value of specific volume (m**3/kg)
!
! general constants
real (kind=kind_phys), parameter :: &
sec_in_day = 86400. &
,sec_in_hour = 3600. &
,solar_time_6am = 21600.0 &
,const_rot = 0.000073 & !< constant to calculate corioli force
,ri_c = 0.65 &
,ri_g = 0.25 &
,eps_z_w = 0.01 & !< criteria to finish iterations for z_w
,eps_conv = 0.01 & !< criteria to finish iterations for d_conv
,eps_sfs = 0.01 & !< criteria to finish iterations for d_sfs
,z_w_max = 20.0 & !< max warm layer thickness
,z_w_min = 0.2 & !< min warm layer thickness
,z_w_ini = 0.2 & !< initial warm layer thickness in dtl_onset
,z_c_max = 0.01 & !< maximum of sub-layer thickness (m)
,z_c_ini = 0.001 & !< initial value of z_c
,ustar_a_min = 0.031 & !< minimum of friction wind speed (m/s): 0.031 ~ 1m/s at 10 m hight
,tau_min = 0.005 & !< minimum of wind stress for dtm
,exp_const = 9.5 & !< coefficient in exponet profile
,delz = 0.1 & !< vertical increment for integral calculation (m)
,von = 0.4 & !< von karman's "constant"
,t0k = 273.16 & !< celsius to kelvin
,gray = 0.97 &
,sst_max = 308.16 &
,tw_max = 5.0 &
,wd_max = 2.0 &
,omg_m = 1.0 & !< trace factor to apply salinity effect
,omg_rot = 1.0 & !< trace factor to apply rotation effect
,omg_sh = 1.0 & !< trace factor to apply sensible heat due to rainfall effect
,visw = 1.e-6 & !< m2/s kinematic viscosity water
,novalue = 0 &
,smallnumber = 1.e-6 &
,timestep_oc = sec_in_day/8. & !< time step in the ocean model (3 hours)
,radian = 2.*pi/180. &
,rad2deg = 180./pi &
,cp_w = 4000. & !< specific heat water (j/kg/k )
,rho0_w = 1022.0 & !< density water (kg/m3 ) (or 1024.438)
,vis_w = 1.e-6 & !< kinematic viscosity water (m2/s )
,tc_w = 0.6 & !< thermal conductivity water (w/m/k )
,capa_w = 3950.0 & !< heat capacity of sea water !
,thref = 1.0e-3 !< reference value of specific volume (m**3/kg)

!!$!============================================
!!$
!!$ ,lvapor=2.453e6 & ! latent heat of vaporization note: make it function of t ????? note the same as hvap
!!$ ,lvapor=2.453e6 & ! latent heat of vaporization note: make it function of t ????? note the same as hvap
!!$ ,alpha=1 ! thermal expansion coefficient
!!$ ,beta ! saline contraction coefficient
!!$ ,cp=1 !=1 specific heat of sea water
Expand All @@ -95,7 +108,7 @@ module module_nst_parameters
!!$ fdg=1.00 !based on results from flux workshop august 1995
!!$ tok=273.16 ! celsius to kelvin
!!$ twopi=3.14159*2.
!!$
!!$
!!$c air constants and coefficients
!!$ rgas=287.1 !j/kg/k gas const. dry air
!!$ xlv=(2.501-0.00237*ts)*1e+6 !j/kg latent heat of vaporization at ts
Expand All @@ -104,7 +117,7 @@ module module_nst_parameters
!!$ rhoa=p*100./(rgas*(t+tok)*(1.+.61*q)) !kg/m3 moist air density ( " )
!!$ visa=1.326e-5*(1+6.542e-3*t+8.301e-6*t*t-4.84e-9*t*t*t) !m2/s
!!$ !kinematic viscosity of dry air - andreas (1989) crrel rep. 89-11
!!$c
!!$c
!!$c cool skin constants
!!$ al=2.1e-5*(ts+3.2)**0.79 !water thermal expansion coefft.
!!$ be=0.026 !salinity expansion coefft.
Expand All @@ -126,11 +139,11 @@ module module_nst_parameters
!!$ real, parameter :: rhoref = 1024.438 !sea water reference density, kg/m^3
!!$ real , parameter :: hslab=50.0 !slab ocean depth
!!$ real , parameter :: bad=-1.0e+10
!!$ real , parameter :: tmin=2.68e+02
!!$ real , parameter :: tmin=2.68e+02
!!$ real , parameter :: tmax=3.11e+02
!!$
!!$ real, parameter :: grav =9.81 !gravity, kg/m/s^2
!!$ real, parameter :: capa =3950.0 !heat capacity of sea water
!!$ real, parameter :: capa =3950.0 !heat capacity of sea water
!!$ real, parameter :: rhoref = 1024.438 !sea water reference density, kg/m^3
!!$ real, parameter :: tmin=2.68e+02 !normal minimal temp
!!$ real, parameter :: tmax=3.11e+02 !normal max temp
Expand Down
Loading
Loading