SCALE-RM
Functions/Subroutines | Variables
mod_mkinit Module Reference

module INITIAL More...

Functions/Subroutines

subroutine, public mkinit_setup
 Setup. More...
 
subroutine, public mkinit (output)
 Driver. More...
 
subroutine rect_setup
 Bubble. More...
 
real(rp) function faero (f0, r0, x, alpha, rhoa)
 
subroutine flux_setup
 flux setup More...
 
subroutine land_setup
 Land setup. More...
 
subroutine ocean_setup
 Ocean setup. More...
 
subroutine urban_setup
 Urban setup. More...
 
subroutine tke_setup
 TKE setup. More...
 
subroutine read_sounding (DENS, VELX, VELY, POTT, QV)
 Read sounding data from file. More...
 

Variables

integer, public mkinit_type = -1
 
integer, parameter, public i_ignore = 0
 
integer, parameter, public i_planestate = 1
 
integer, parameter, public i_tracerbubble = 2
 
integer, parameter, public i_coldbubble = 3
 
integer, parameter, public i_lambwave = 4
 
integer, parameter, public i_gravitywave = 5
 
integer, parameter, public i_khwave = 6
 
integer, parameter, public i_turbulence = 7
 
integer, parameter, public i_mountainwave = 8
 
integer, parameter, public i_warmbubble = 9
 
integer, parameter, public i_supercell = 10
 
integer, parameter, public i_squallline = 11
 
integer, parameter, public i_wk1982 = 12
 
integer, parameter, public i_dycoms2_rf01 = 13
 
integer, parameter, public i_dycoms2_rf02 = 14
 
integer, parameter, public i_rico = 15
 
integer, parameter, public i_interporation = 16
 
integer, parameter, public i_landcouple = 17
 
integer, parameter, public i_oceancouple = 18
 
integer, parameter, public i_urbancouple = 19
 
integer, parameter, public i_triplecouple = 20
 
integer, parameter, public i_bubblecouple = 21
 
integer, parameter, public i_seabreeze = 22
 
integer, parameter, public i_heatisland = 23
 
integer, parameter, public i_dycoms2_rf02_dns = 24
 
integer, parameter, public i_real = 25
 
integer, parameter, public i_grayzone = 26
 
integer, parameter, public i_boxaero = 27
 
integer, parameter, public i_warmbubbleaero = 28
 
integer, parameter, public i_cavityflow = 29
 
integer, parameter, public i_barocwave = 30
 
integer, parameter, public i_bomex = 31
 

Detailed Description

module INITIAL

Description
subroutines for preparing initial data
Author
Team SCALE
NAMELIST
  • PARAM_MKINIT
    nametypedefault valuecomment
    MKINIT_INITNAME character(len=H_SHORT) 'NONE'

  • PARAM_BUBBLE
    nametypedefault valuecomment
    BBL_EACHNODE logical .false. Arrange bubble at each node? [kg/kg]
    BBL_CZ real(RP) 2.E3_RP center location [m]: z
    BBL_CX real(RP) 2.E3_RP center location [m]: x
    BBL_CY real(RP) 2.E3_RP center location [m]: y
    BBL_RZ real(RP) 0.0_RP bubble radius [m]: z
    BBL_RX real(RP) 0.0_RP bubble radius [m]: x
    BBL_RY real(RP) 0.0_RP bubble radius [m]: y
    BBL_FUNCTYPE character(len=H_SHORT) 'COSBELL' COSBELL or GAUSSIAN

  • PARAM_RECT
    nametypedefault valuecomment
    RCT_EACHNODE logical .false. Arrange rectangle at each node? [kg/kg]
    RCT_CZ real(RP) 2.E3_RP center location [m]: z
    RCT_CX real(RP) 2.E3_RP center location [m]: x
    RCT_CY real(RP) 2.E3_RP center location [m]: y
    RCT_RZ real(RP) 2.E3_RP rectangle z width [m]: z
    RCT_RX real(RP) 2.E3_RP rectangle x width [m]: x
    RCT_RY real(RP) 2.E3_RP rectangle y width [m]: y

  • PARAM_AERO
    nametypedefault valuecomment
    CCN_INIT real(RP) 50.E+6_RP initial cloud condensation nucrei [#/m3]
    M0_INIT real(RP) 0.0_RP initial total num. conc. of modes (Atk,Acm,Cor) [#/m3]
    DG_INIT real(RP) 80.e-9_RP initial number equivalen diameters of modes [m]
    SG_INIT real(RP) 1.6_RP initial standard deviation [-]
    D_MIN_INP real(RP), dimension(3) d_min_def
    D_MAX_INP real(RP), dimension(3) d_max_def
    K_MIN_INP real(RP), dimension(3) k_min_def
    K_MAX_INP real(RP), dimension(3) k_max_def
    N_KAP_INP integer, dimension(3) n_kap_def

  • PARAM_MKINIT_FLUX
    nametypedefault valuecomment
    FLX_RAIN real(RP) 0.0_RP surface rain flux [kg/m2/s]
    FLX_SNOW real(RP) 0.0_RP surface snow flux [kg/m2/s]
    FLX_IR_DN real(RP) 0.0_RP surface downwad radiation flux [J/m2/s]
    FLX_NIR_DN real(RP) 0.0_RP surface downwad radiation flux [J/m2/s]
    FLX_VIS_DN real(RP) 0.0_RP surface downwad radiation flux [J/m2/s]

  • PARAM_MKINIT_LAND
    nametypedefault valuecomment
    LND_TEMP real(RP) land soil temperature [K]
    LND_WATER real(RP) 0.15_RP land soil moisture [m3/m3]
    LND_ICE real(RP) 0.00_RP land soil ice [m3/m3]
    SFC_TEMP real(RP) ocean skin temperature [K]
    SFC_ALBEDO_LW real(RP) 0.04_RP ocean surface albedo for LW (0-1)
    SFC_ALBEDO_SW real(RP) 0.05_RP ocean surface albedo for SW (0-1)

  • PARAM_MKINIT_OCEAN
    nametypedefault valuecomment
    OCN_TEMP real(RP) ocean temperature [K]
    OCN_SALT real(RP) 0.0_RP ocean salinity [psu]
    OCN_UVEL real(RP) 0.0_RP ocean u-velocity [m/s]
    OCN_VVEL real(RP) 0.0_RP ocean v-velocity [m/s]
    ICE_TEMP real(RP) ocean temperature [K]
    ICE_MASS real(RP) 0.0_RP ocean temperature [K]
    SFC_TEMP real(RP) ocean skin temperature [K]
    SFC_ALBEDO_LW real(RP) 0.04_RP ocean surface albedo for LW (0-1)
    SFC_ALBEDO_SW real(RP) 0.05_RP ocean surface albedo for SW (0-1)
    SFC_Z0M real(RP) 1.E-4_RP ocean surface roughness length (momentum) [m]
    SFC_Z0H real(RP) 1.E-4_RP ocean surface roughness length (heat) [m]
    SFC_Z0E real(RP) 1.E-4_RP ocean surface roughness length (vapor) [m]

  • PARAM_MKINIT_URBAN
    nametypedefault valuecomment
    URB_ROOF_TEMP real(RP) Surface temperature of roof [K]
    URB_BLDG_TEMP real(RP) Surface temperature of building [K]
    URB_GRND_TEMP real(RP) Surface temperature of ground [K]
    URB_CNPY_TEMP real(RP) Diagnostic canopy air temperature [K]
    URB_CNPY_HMDT real(RP) 0.0_RP Diagnostic canopy humidity [kg/kg]
    URB_CNPY_WIND real(RP) 0.0_RP Diagnostic canopy wind [m/s]
    URB_ROOF_LAYER_TEMP real(RP) temperature in layer of roof [K]
    URB_BLDG_LAYER_TEMP real(RP) temperature in layer of building [K]
    URB_GRND_LAYER_TEMP real(RP) temperature in layer of ground [K]
    URB_ROOF_RAIN real(RP) 0.0_RP temperature in layer of roof [kg/m2]
    URB_BLDG_RAIN real(RP) 0.0_RP temperature in layer of building [kg/m2]
    URB_GRND_RAIN real(RP) 0.0_RP temperature in layer of ground [kg/m2]
    URB_SFC_TEMP real(RP) Grid average of surface temperature [K]
    URB_ALB_LW real(RP) 0.10_RP Grid average of surface albedo for LW (0-1)
    URB_ALB_SW real(RP) 0.20_RP Grid average of surface albedo for SW (0-1)

  • PARAM_MKINIT_TKE
    nametypedefault valuecomment
    TKE_CONST real(RP)
    ZI_CONST real(RP)

  • PARAM_MKINIT_SOUNDING
    nametypedefault valuecomment
    ENV_IN_SOUNDING_FILE character(len=H_LONG) ''

  • PARAM_MKINIT_PLANESTATE
    nametypedefault valuecomment
    SFC_THETA real(RP) surface potential temperature [K]
    SFC_PRES real(RP) surface pressure [Pa]
    SFC_RH real(RP) 80.0_RP surface relative humidity [%]
    ENV_THETA real(RP) potential temperature of environment [K]
    ENV_TLAPS real(RP) 4.E-3_RP Lapse rate of THETA [K/m]
    ENV_U real(RP) 0.0_RP velocity u of environment [m/s]
    ENV_V real(RP) 0.0_RP velocity v of environment [m/s]
    ENV_RH real(RP) 80.0_RP Relative Humidity of environment [%]
    RANDOM_THETA real(RP) 0.01_RP
    RANDOM_U real(RP) 0.0_RP amplitude of random disturbance u
    RANDOM_V real(RP) 0.0_RP amplitude of random disturbance v
    RANDOM_RH real(RP) 0.0_RP amplitude of random disturbance RH

  • PARAM_MKINIT_TRACERBUBBLE
    nametypedefault valuecomment
    SFC_THETA real(RP) surface potential temperature [K]
    SFC_PRES real(RP) surface pressure [Pa]
    ENV_THETA real(RP) potential temperature of environment [K]
    ENV_U real(RP) 0.0_RP velocity u of environment [m/s]
    ENV_V real(RP) 0.0_RP velocity v of environment [m/s]
    SHAPE_PTRACER character(len=H_SHORT) 'BUBBLE' BUBBLE or RECT
    BBL_PTRACER real(RP) 0.0_RP extremum of passive tracer in bubble [kg/kg]

  • PARAM_MKINIT_COLDBUBBLE
    nametypedefault valuecomment
    SFC_THETA real(RP) surface potential temperature [K]
    SFC_PRES real(RP) surface pressure [Pa]
    ENV_THETA real(RP) potential temperature of environment [K]
    BBL_TEMP real(RP) -15.0_RP extremum of temperature in bubble [K]

  • PARAM_MKINIT_LAMBWAVE
    nametypedefault valuecomment
    SFC_PRES real(RP) surface pressure [Pa]
    ENV_U real(RP) 0.0_RP velocity u of environment [m/s]
    ENV_V real(RP) 0.0_RP velocity v of environment [m/s]
    ENV_TEMP real(RP) 300.0_RP temperature of environment [K]
    BBL_PRES real(RP) 100._RP extremum of pressure in bubble [Pa]

  • PARAM_MKINIT_GRAVITYWAVE
    nametypedefault valuecomment
    SFC_THETA real(RP) surface potential temperature [K]
    SFC_PRES real(RP) surface pressure [Pa]
    ENV_U real(RP) 0.0_RP velocity u of environment [m/s]
    ENV_V real(RP) 0.0_RP velocity v of environment [m/s]
    ENV_BVF real(RP) 0.01_RP Brunt Vaisala frequencies of environment [1/s]
    BBL_THETA real(RP) 1.0_RP extremum of temperature in bubble [K]

  • PARAM_MKINIT_KHWAVE
    nametypedefault valuecomment
    SFC_THETA real(RP) surface potential temperature [K]
    SFC_PRES real(RP) surface pressure [Pa]
    ENV_L1_ZTOP real(RP) 1.E3_RP top height of the layer1 (constant THETA) [m]
    ENV_L3_ZBOTTOM real(RP) 2100.0_RP bottom height of the layer3 (high THETA) [m]
    ENV_L1_THETA real(RP) 300.0_RP THETA in the layer1 (low THETA) [K]
    ENV_L3_THETA real(RP) 301.0_RP THETA in the layer3 (high THETA) [K]
    ENV_L1_U real(RP) 0.0_RP velocity u in the layer1 (low THETA) [K]
    ENV_L3_U real(RP) 20.0_RP velocity u in the layer3 (high THETA) [K]
    RANDOM_U real(RP) 0.0_RP amplitude of random disturbance u

  • PARAM_MKINIT_TURBULENCE
    nametypedefault valuecomment
    SFC_THETA real(RP) surface potential temperature [K]
    SFC_PRES real(RP) surface pressure [Pa]
    SFC_RH real(RP) 80.0_RP surface relative humidity [%]
    ENV_THETA real(RP) potential temperature of environment [K]
    ENV_TLAPS real(RP) 4.E-3_RP Lapse rate of THETA [K/m]
    ENV_U real(RP) 0.0_RP velocity u of environment [m/s]
    ENV_V real(RP) 0.0_RP velocity v of environment [m/s]
    ENV_RH real(RP) 80.0_RP Relative Humidity of environment [%]
    RANDOM_THETA real(RP) 0.01_RP
    RANDOM_U real(RP) 0.0_RP amplitude of random disturbance u
    RANDOM_V real(RP) 0.0_RP amplitude of random disturbance v
    RANDOM_RH real(RP) 0.0_RP amplitude of random disturbance RH

  • PARAM_MKINIT_CAVITYFLOW
    nametypedefault valuecomment
    ULID real(RP) 1.D01
    PRES0 real(RP) 1.D05
    REYNOLDS_NUM real(RP) 1.D03
    MACH_NUM real(RP) 3.D-2

  • PARAM_MKINIT_MOUNTAINWAVE
    nametypedefault valuecomment
    SFC_THETA real(RP) surface potential temperature [K]
    SFC_PRES real(RP) surface pressure [Pa]
    ENV_U real(RP) 0.0_RP velocity u of environment [m/s]
    ENV_V real(RP) 0.0_RP velocity v of environment [m/s]
    SCORER real(RP) 2.E-3_RP Scorer parameter (~=N/U) [1/m]
    BBL_PTRACER real(RP) 0.0_RP extremum of passive tracer in bubble [kg/kg]

  • PARAM_MKINIT_BAROCWAVE
    nametypedefault valuecomment
    REF_TEMP real(RP) 288.E0_RP The reference temperature [K]
    REF_PRES real(RP) 1.E5_RP The reference pressure [Pa]
    LAPSE_RATE real(RP) 5.E-3_RP The lapse rate [K/m]
    PHI0DEG real(RP) 45.E0_RP The central latitude [degree_north]
    U0 real(RP) 35.E0_RP The parameter associated with zonal jet maximum amplitude [m/s]
    B real(RP) 2.E0_RP The vertical half-width [1]
    UP real(RP) 1.E0_RP The maximum amplitude of zonal wind perturbation [m/s]
    LP real(RP) 600.E3_RP The width of Gaussian profile
    XC real(RP) 2000.E3_RP The center point (x) of inital perturbation
    YC real(RP) 2500.E3_RP The center point (y) of inital perturbation

  • PARAM_MKINIT_WARMBUBBLE
    nametypedefault valuecomment
    SFC_THETA real(RP) surface potential temperature [K]
    SFC_PRES real(RP) surface pressure [Pa]
    ENV_U real(RP) 0.0_RP velocity u of environment [m/s]
    ENV_V real(RP) 0.0_RP velocity v of environment [m/s]
    ENV_RH real(RP) 80.0_RP Relative Humidity of environment [%]
    ENV_L1_ZTOP real(RP) 1.E3_RP top height of the layer1 (constant THETA) [m]
    ENV_L2_ZTOP real(RP) 14.E3_RP top height of the layer2 (small THETA gradient) [m]
    ENV_L2_TLAPS real(RP) 4.E-3_RP Lapse rate of THETA in the layer2 (small THETA gradient) [K/m]
    ENV_L3_TLAPS real(RP) 3.E-2_RP Lapse rate of THETA in the layer3 (large THETA gradient) [K/m]
    BBL_THETA real(RP) 1.0_RP extremum of temperature in bubble [K]

  • PARAM_MKINIT_SUPERCELL
    nametypedefault valuecomment
    BBL_THETA real(RP) 1.0_RP extremum of temperature in bubble [K]

  • PARAM_MKINIT_SQUALLLINE
    nametypedefault valuecomment
    RANDOM_THETA real(RP) 0.01_RP
    OFFSET_VELX real(RP) 12.0_RP
    OFFSET_VELY real(RP) -2.0_RP

  • PARAM_MKINIT_WK1982
    nametypedefault valuecomment
    SFC_THETA real(RP) surface potential temperature [K]
    SFC_PRES real(RP) surface pressure [Pa]
    TR_Z real(RP) 12000.0_RP height of tropopause [m]
    TR_THETA real(RP) 343.0_RP pot. temperature at tropopause [K]
    TR_TEMP real(RP) 213.0_RP temperature at tropopause [K]
    SHEAR_Z real(RP) 3000.0_RP center height of shear layer [m]
    SHEAR_U real(RP) 15.0_RP velocity u over the shear layer [m/s]
    QV0 real(RP) 14.0_RP maximum vapor mixing ration [g/kg]
    BBL_THETA real(RP) 1.0_RP extremum of temperature in bubble [K]

  • PARAM_MKINIT_RF01
    nametypedefault valuecomment
    PERTURB_AMP real(RP) 0.0_RP
    RANDOM_LIMIT integer 0
    RANDOM_FLAG integer 0 0 -> no perturbation
    USE_LWSET logical .false. use liq. water. static energy temp.?

  • PARAM_MKINIT_RF02
    nametypedefault valuecomment
    PERTURB_AMP real(RP) 0.0_RP
    RANDOM_LIMIT integer 0
    RANDOM_FLAG integer 0 0 -> no perturbation

  • PARAM_MKINIT_RF02_DNS
    nametypedefault valuecomment
    ZB real(RP) 750.0_RP domain bottom
    CONST_U real(RP) 0.0_RP
    CONST_V real(RP) 0.0_RP
    PRES_ZB real(RP) 93060.0_RP
    PERTURB_AMP real(RP) 0.0_RP
    RANDOM_LIMIT integer 0
    RANDOM_FLAG integer 0 0 -> no perturbation

  • PARAM_MKINIT_RICO
    nametypedefault valuecomment
    PERTURB_AMP_PT real(RP) 0.1_RP
    PERTURB_AMP_QV real(RP) 2.5E-5_RP

  • PARAM_MKINIT_BOMEX
    nametypedefault valuecomment
    PERTURB_AMP_PT real(RP) 0.1_RP
    PERTURB_AMP_QV real(RP) 2.5E-5_RP

  • PARAM_MKINIT_SEABREEZE
    nametypedefault valuecomment
    LAND_SIZE real(RP)

  • PARAM_MKINIT_GRAYZONE
    nametypedefault valuecomment
    PERTURB_AMP real(RP) 0.0_RP
    RANDOM_LIMIT integer 0
    RANDOM_FLAG integer 0 0 -> no perturbation

  • PARAM_MKINIT_BOXAERO
    nametypedefault valuecomment
    INIT_DENS real(RP) 1.12_RP [kg/m3]
    INIT_TEMP real(RP) 298.18_RP [K]
    INIT_PRES real(RP) 1.E+5_RP [Pa]
    INIT_SSLIQ real(RP) 0.01_RP [%]

History Output
No history output

Function/Subroutine Documentation

◆ mkinit_setup()

subroutine, public mod_mkinit::mkinit_setup

Setup.

Definition at line 211 of file mod_mkinit.F90.

211  implicit none
212 
213  character(len=H_SHORT) :: MKINIT_initname = 'NONE'
214 
215  namelist / param_mkinit / &
216  mkinit_initname
217 
218  integer :: ierr
219  !---------------------------------------------------------------------------
220 
221  log_newline
222  log_info("MKINIT_setup",*) 'Setup'
223 
224  !--- read namelist
225  rewind(io_fid_conf)
226  read(io_fid_conf,nml=param_mkinit,iostat=ierr)
227  if( ierr < 0 ) then !--- missing
228  log_info("MKINIT_setup",*) 'Not found namelist. Default used.'
229  elseif( ierr > 0 ) then !--- fatal error
230  log_error("MKINIT_setup",*) 'Not appropriate names in namelist PARAM_MKINIT. Check!'
231  call prc_abort
232  endif
233  log_nml(param_mkinit)
234 
235  allocate( pres(ka,ia,ja) )
236  allocate( temp(ka,ia,ja) )
237  allocate( pott(ka,ia,ja) )
238  allocate( qdry(ka,ia,ja) )
239  allocate( qsat(ka,ia,ja) )
240  allocate( qv(ka,ia,ja) )
241  allocate( qc(ka,ia,ja) )
242  allocate( nc(ka,ia,ja) )
243  allocate( velx(ka,ia,ja) )
244  allocate( vely(ka,ia,ja) )
245  allocate( ptrc(ka,ia,ja) )
246 
247  allocate( pres_sfc(ia,ja) )
248  allocate( temp_sfc(ia,ja) )
249  allocate( pott_sfc(ia,ja) )
250  allocate( psat_sfc(ia,ja) )
251  allocate( qsat_sfc(ia,ja) )
252  allocate( qv_sfc(ia,ja) )
253  allocate( qc_sfc(ia,ja) )
254 
255  allocate( rndm(ka,ia,ja) )
256  allocate( bubble(ka,ia,ja) )
257  allocate( rect(ka,ia,ja) )
258 
259 
260  select case(trim(mkinit_initname))
261  case('NONE')
262  mkinit_type = i_ignore
263  case('PLANESTATE')
264  mkinit_type = i_planestate
265  case('TRACERBUBBLE')
266  mkinit_type = i_tracerbubble
267  case('COLDBUBBLE')
268  mkinit_type = i_coldbubble
269  call bubble_setup
270  case('LAMBWAVE')
271  mkinit_type = i_lambwave
272  call bubble_setup
273  case('GRAVITYWAVE')
274  mkinit_type = i_gravitywave
275  call bubble_setup
276  case('KHWAVE')
277  mkinit_type = i_khwave
278  case('TURBULENCE')
279  mkinit_type = i_turbulence
280  case('MOUNTAINWAVE')
281  mkinit_type = i_mountainwave
282  call bubble_setup
283  case('WARMBUBBLE')
284  mkinit_type = i_warmbubble
285  call bubble_setup
286  case('SUPERCELL')
287  mkinit_type = i_supercell
288  call bubble_setup
289  case('SQUALLLINE')
290  mkinit_type = i_squallline
291  case('WK1982')
292  mkinit_type = i_wk1982
293  call bubble_setup
294  case('DYCOMS2_RF01')
295  mkinit_type = i_dycoms2_rf01
296  case('DYCOMS2_RF02')
297  mkinit_type = i_dycoms2_rf02
298  case('RICO')
299  mkinit_type = i_rico
300  case('BOMEX')
301  mkinit_type = i_bomex
302  case('INTERPORATION')
303  mkinit_type = i_interporation
304  case('LANDCOUPLE')
305  mkinit_type = i_landcouple
306  case('OCEANCOUPLE')
307  mkinit_type = i_oceancouple
308  case('URBANCOUPLE')
309  mkinit_type = i_urbancouple
310  case('TRIPLECOUPLE')
311  mkinit_type = i_triplecouple
312  case('BUBBLECOUPLE')
313  mkinit_type = i_bubblecouple
314  call bubble_setup
315  case('SEABREEZE')
316  mkinit_type = i_seabreeze
317  case('HEATISLAND')
318  mkinit_type = i_heatisland
319  case('DYCOMS2_RF02_DNS')
320  mkinit_type = i_dycoms2_rf02_dns
321  case('REAL')
322  mkinit_type = i_real
323  case('GRAYZONE')
324  mkinit_type = i_grayzone
325  case('BOXAERO')
326  mkinit_type = i_boxaero
327  case('WARMBUBBLEAERO')
328  mkinit_type = i_warmbubbleaero
329  call bubble_setup
330  case('CAVITYFLOW')
331  mkinit_type = i_cavityflow
332  case('BAROCWAVE')
333  mkinit_type = i_barocwave
334  case default
335  log_error("MKINIT_setup",*) 'Unsupported TYPE:', trim(mkinit_initname)
336  call prc_abort
337  endselect
338 
339  return

References i_barocwave, i_bomex, i_boxaero, i_bubblecouple, i_cavityflow, i_coldbubble, i_dycoms2_rf01, i_dycoms2_rf02, i_dycoms2_rf02_dns, i_gravitywave, i_grayzone, i_heatisland, i_ignore, i_interporation, i_khwave, i_lambwave, i_landcouple, i_mountainwave, i_oceancouple, i_planestate, i_real, i_rico, i_seabreeze, i_squallline, i_supercell, i_tracerbubble, i_triplecouple, i_turbulence, i_urbancouple, i_warmbubble, i_warmbubbleaero, i_wk1982, scale_atmos_grid_cartesc_index::ia, scale_io::io_fid_conf, scale_atmos_grid_cartesc_index::ja, scale_atmos_grid_cartesc_index::ka, mkinit_type, and scale_prc::prc_abort().

Referenced by mod_rm_prep::rm_prep().

Here is the call graph for this function:
Here is the caller graph for this function:

◆ mkinit()

subroutine, public mod_mkinit::mkinit ( logical, intent(out)  output)

Driver.

Definition at line 345 of file mod_mkinit.F90.

345  use scale_const, only: &
346  const_undef8
347  use scale_atmos_hydrometeor, only: &
348  atmos_hydrometeor_dry, &
349  n_hyd, &
350  i_hc
351  use mod_atmos_phy_mp_vars, only: &
352  qs_mp, &
353  qe_mp
354  use mod_atmos_admin, only: &
355  atmos_phy_mp_type
356  use mod_atmos_phy_mp_driver, only: &
357  atmos_phy_mp_driver_qhyd2qtrc
358  implicit none
359  logical, intent(out) :: output
360 
361  real(RP) :: QHYD(KA,IA,JA,N_HYD)
362  real(RP) :: QNUM(KA,IA,JA,N_HYD)
363 
364  logical :: convert_qtrc
365  integer :: iq
366  !---------------------------------------------------------------------------
367 
368  if ( mkinit_type == i_ignore ) then
369  log_newline
370  log_progress(*) 'skip making initial data'
371  output = .false.
372  else
373  log_newline
374  log_progress(*) 'start making initial data'
375 
376  !--- Initialize variables
377  pres(:,:,:) = const_undef8
378  temp(:,:,:) = const_undef8
379  pott(:,:,:) = const_undef8
380  qsat(:,:,:) = const_undef8
381  velx(:,:,:) = const_undef8
382  vely(:,:,:) = const_undef8
383 
384  rndm(:,:,:) = const_undef8
385 
386  pres_sfc(:,:) = const_undef8
387  temp_sfc(:,:) = const_undef8
388  pott_sfc(:,:) = const_undef8
389  psat_sfc(:,:) = const_undef8
390  qsat_sfc(:,:) = const_undef8
391 
392  qv(:,:,:) = 0.0_rp
393  qc(:,:,:) = 0.0_rp
394  nc(:,:,:) = 0.0_rp
395  qv_sfc(:,:) = 0.0_rp
396  qc_sfc(:,:) = 0.0_rp
397 
398  ptrc(:,:,:) = const_undef8
399 
400 !OCL XFILL
401  qtrc(:,:,:,:) = 0.0_rp
402 !OCL XFILL
403  qhyd(:,:,:,:) = 0.0_rp
404 !OCL XFILL
405  qnum(:,:,:,:) = 0.0_rp
406 
407  call prof_rapstart('_MkInit_main',3)
408 
409  convert_qtrc = .true.
410 
411  select case(mkinit_type)
412  case(i_planestate)
413  call mkinit_planestate
414  case(i_tracerbubble)
415  call mkinit_tracerbubble
416  case(i_coldbubble)
417  call mkinit_coldbubble
418  case(i_lambwave)
419  call mkinit_lambwave
420  case(i_gravitywave)
421  call mkinit_gravitywave
422  case(i_khwave)
423  call mkinit_khwave
424  case(i_turbulence)
425  call mkinit_turbulence
426  case(i_mountainwave)
427  call mkinit_mountainwave
428  case(i_warmbubble)
429  call mkinit_warmbubble
430  case(i_supercell)
431  call mkinit_supercell
432  case(i_squallline)
433  call mkinit_squallline
434  case(i_wk1982)
435  call mkinit_wk1982
436  case(i_dycoms2_rf01)
437  call mkinit_dycoms2_rf01
438  case(i_dycoms2_rf02)
439  call mkinit_dycoms2_rf02
440  case(i_rico)
441  call mkinit_rico
442  case(i_bomex)
443  call mkinit_bomex
444  case(i_oceancouple)
445  call mkinit_planestate
446  call mkinit_oceancouple
447  case(i_landcouple)
448  call mkinit_planestate
449  call mkinit_landcouple
450  case(i_urbancouple)
451  call mkinit_planestate
452  call mkinit_urbancouple
453  case(i_triplecouple)
454  call mkinit_planestate
455  call mkinit_oceancouple
456  call mkinit_landcouple
457  call mkinit_urbancouple
458  case(i_bubblecouple)
459  call mkinit_planestate
460  call mkinit_warmbubble
461  call mkinit_oceancouple
462  call mkinit_landcouple
463  call mkinit_urbancouple
464  case(i_seabreeze)
465  call mkinit_planestate
466  call mkinit_seabreeze
467  case(i_heatisland)
468  call mkinit_planestate
469  call mkinit_heatisland
470  case(i_dycoms2_rf02_dns)
471  call mkinit_dycoms2_rf02_dns
472  case(i_real)
473  call mkinit_real
474  convert_qtrc = .false.
475  case(i_grayzone)
476  call mkinit_grayzone
477  case(i_boxaero)
478  call mkinit_boxaero
479  case(i_warmbubbleaero)
480  call mkinit_warmbubbleaero
481  case(i_cavityflow)
482  call mkinit_cavityflow
483  case(i_barocwave)
484  call mkinit_barocwave
485  case default
486  log_error("MKINIT",*) 'Unsupported TYPE:', mkinit_type
487  call prc_abort
488  endselect
489 
490  call tke_setup
491 
492  call aerosol_setup
493 
494  call sbmaero_setup( convert_qtrc ) ! [INOUT]
495 
496  ! water content
497  if ( ( .not. atmos_hydrometeor_dry ) .AND. convert_qtrc ) then
498 !OCL XFILL
499  qhyd(:,:,:,i_hc) = qc(:,:,:)
500 !OCL XFILL
501  qnum(:,:,:,i_hc) = nc(:,:,:)
502  call atmos_phy_mp_driver_qhyd2qtrc( ka, ks, ke, ia, is, ie, ja, js, je, &
503  qv(:,:,:), qhyd(:,:,:,:), & ! [IN]
504  qtrc(:,:,:,qs_mp:qe_mp), & ! [OUT]
505  qnum=qnum(:,:,:,:) ) ! [IN]
506  end if
507 
508  ! passive tracer
509  call tracer_inq_id( "PTracer", iq)
510  if ( iq > 0 ) qtrc(:,:,:,iq) = ptrc(:,:,:)
511 
512  call prof_rapend ('_MkInit_main',3)
513 
514  log_progress(*) 'end making initial data'
515 
516  output = .true.
517 
518  endif
519 
520  return

References scale_atmos_hydrometeor::atmos_hydrometeor_dry, mod_atmos_phy_mp_driver::atmos_phy_mp_driver_qhyd2qtrc(), mod_atmos_admin::atmos_phy_mp_type, scale_const::const_undef8, i_barocwave, i_bomex, i_boxaero, i_bubblecouple, i_cavityflow, i_coldbubble, i_dycoms2_rf01, i_dycoms2_rf02, i_dycoms2_rf02_dns, i_gravitywave, i_grayzone, scale_atmos_hydrometeor::i_hc, i_heatisland, i_ignore, i_khwave, i_lambwave, i_landcouple, i_mountainwave, i_oceancouple, i_planestate, i_real, i_rico, i_seabreeze, i_squallline, i_supercell, i_tracerbubble, i_triplecouple, i_turbulence, i_urbancouple, i_warmbubble, i_warmbubbleaero, i_wk1982, scale_atmos_grid_cartesc_index::ia, scale_atmos_grid_cartesc_index::iag, scale_atmos_grid_cartesc_index::ie, scale_atmos_grid_cartesc_index::ihalo, scale_io::io_fid_conf, scale_atmos_grid_cartesc_index::is, scale_atmos_grid_cartesc_index::ja, scale_atmos_grid_cartesc_index::jag, scale_atmos_grid_cartesc_index::je, scale_atmos_grid_cartesc_index::jhalo, scale_atmos_grid_cartesc_index::js, scale_tracer::k, scale_atmos_grid_cartesc_index::ka, scale_atmos_grid_cartesc_index::ke, scale_atmos_grid_cartesc_index::ks, mkinit_type, scale_atmos_hydrometeor::n_hyd, scale_prc::prc_abort(), scale_prof::prof_rapend(), scale_prof::prof_rapstart(), mod_atmos_phy_mp_vars::qe_mp, mod_atmos_phy_mp_vars::qs_mp, mod_atmos_vars::qtrc, tke_setup(), and scale_tracer::tracer_inq_id().

Referenced by mod_rm_prep::rm_prep().

Here is the call graph for this function:
Here is the caller graph for this function:

◆ rect_setup()

subroutine mod_mkinit::rect_setup

Bubble.

Definition at line 631 of file mod_mkinit.F90.

631  use scale_const, only: &
632  const_undef8
633  implicit none
634 
635  ! Bubble
636  logical :: RCT_eachnode = .false. ! Arrange rectangle at each node? [kg/kg]
637  real(RP) :: RCT_CZ = 2.e3_rp ! center location [m]: z
638  real(RP) :: RCT_CX = 2.e3_rp ! center location [m]: x
639  real(RP) :: RCT_CY = 2.e3_rp ! center location [m]: y
640  real(RP) :: RCT_RZ = 2.e3_rp ! rectangle z width [m]: z
641  real(RP) :: RCT_RX = 2.e3_rp ! rectangle x width [m]: x
642  real(RP) :: RCT_RY = 2.e3_rp ! rectangle y width [m]: y
643 
644  namelist / param_rect / &
645  rct_eachnode, &
646  rct_cz, &
647  rct_cx, &
648  rct_cy, &
649  rct_rz, &
650  rct_rx, &
651  rct_ry
652 
653  real(RP) :: CZ_offset
654  real(RP) :: CX_offset
655  real(RP) :: CY_offset
656  real(RP) :: dist
657 
658  integer :: ierr
659  integer :: k, i, j
660  !---------------------------------------------------------------------------
661 
662  log_newline
663  log_info("RECT_setup",*) 'Setup'
664 
665  !--- read namelist
666  rewind(io_fid_conf)
667  read(io_fid_conf,nml=param_rect,iostat=ierr)
668  if( ierr < 0 ) then !--- missing
669  log_error("RECT_setup",*) 'Not found namelist. Check!'
670  call prc_abort
671  elseif( ierr > 0 ) then !--- fatal error
672  log_error("RECT_setup",*) 'Not appropriate names in namelist PARAM_RECT. Check!'
673  call prc_abort
674  endif
675  log_nml(param_rect)
676 
677  rect(:,:,:) = const_undef8
678 
679  if ( rct_eachnode ) then
680  cz_offset = cz(ks)
681  cx_offset = cx(is)
682  cy_offset = cy(js)
683  else
684  cz_offset = 0.0_rp
685  cx_offset = 0.0_rp
686  cy_offset = 0.0_rp
687  endif
688 
689  do j = 1, ja
690  do i = 1, ia
691  do k = ks, ke
692 
693  ! make tracer rectangle
694  dist = 2.0_rp * max( &
695  abs(cz(k) - cz_offset - rct_cz)/rct_rz, &
696  abs(cx(i) - cx_offset - rct_cx)/rct_rx, &
697  abs(cy(j) - cy_offset - rct_cy)/rct_ry &
698  & )
699  if ( dist <= 1.0_rp ) then
700  rect(k,i,j) = 1.0_rp
701  else
702  rect(k,i,j) = 0.0_rp
703  end if
704  enddo
705  enddo
706  enddo
707 
708  return

References scale_atmos_phy_ae_kajino13::atmos_phy_ae_kajino13_mkinit(), scale_atmos_phy_ae_offline::atmos_phy_ae_offline_mkinit(), mod_atmos_admin::atmos_phy_ae_type, mod_atmos_admin::atmos_phy_mp_type, scale_const::const_undef8, mod_atmos_vars::dens, scale_atmos_hydrometeor::i_qv, scale_atmos_grid_cartesc_index::ia, scale_atmos_grid_cartesc_index::ie, scale_atmos_grid_cartesc_index::ieb, scale_io::io_fid_conf, scale_atmos_grid_cartesc_index::is, scale_atmos_grid_cartesc_index::isb, scale_atmos_grid_cartesc_index::ja, scale_atmos_grid_cartesc_index::je, scale_atmos_grid_cartesc_index::jeb, scale_atmos_grid_cartesc_index::js, scale_atmos_grid_cartesc_index::jsb, scale_atmos_grid_cartesc_index::ka, scale_atmos_grid_cartesc_index::ke, scale_atmos_grid_cartesc_index::ks, scale_atmos_phy_mp_suzuki10::nbin, scale_atmos_phy_mp_suzuki10::nccn, scale_prc::prc_abort(), mod_atmos_phy_ae_vars::qa_ae, mod_atmos_phy_ae_vars::qe_ae, scale_atmos_hydrometeor::qhe, scale_atmos_hydrometeor::qhs, mod_atmos_phy_ae_vars::qs_ae, and mod_atmos_vars::qtrc.

Referenced by read_sounding().

Here is the call graph for this function:
Here is the caller graph for this function:

◆ faero()

real(rp) function mod_mkinit::faero ( real(rp), intent(in)  f0,
real(rp), intent(in)  r0,
real(rp), intent(in)  x,
real(rp), intent(in)  alpha,
real(rp), intent(in)  rhoa 
)

Definition at line 861 of file mod_mkinit.F90.

861  use scale_const, only: &
862  pi => const_pi
863  implicit none
864 
865  real(RP), intent(in) :: x, f0, r0, alpha, rhoa
866  real(RP) :: faero
867  real(RP) :: rad
868  !---------------------------------------------------------------------------
869 
870  rad = ( exp(x) * 3.0_rp / 4.0_rp / pi / rhoa )**(1.0_rp/3.0_rp)
871 
872  faero = f0 * (rad/r0)**(-alpha)
873 
874  return

References scale_const::const_pi.

◆ flux_setup()

subroutine mod_mkinit::flux_setup

flux setup

Definition at line 880 of file mod_mkinit.F90.

880  use mod_atmos_phy_mp_vars, only: &
881  sflx_rain => atmos_phy_mp_sflx_rain, &
882  sflx_snow => atmos_phy_mp_sflx_snow
883  use mod_atmos_phy_rd_vars, only: &
884  sflx_lw_up => atmos_phy_rd_sflx_lw_up, &
885  sflx_lw_dn => atmos_phy_rd_sflx_lw_dn, &
886  sflx_sw_up => atmos_phy_rd_sflx_sw_up, &
887  sflx_sw_dn => atmos_phy_rd_sflx_sw_dn
888  implicit none
889 
890  ! Flux from Atmosphere
891  real(RP) :: FLX_rain = 0.0_rp ! surface rain flux [kg/m2/s]
892  real(RP) :: FLX_snow = 0.0_rp ! surface snow flux [kg/m2/s]
893  real(RP) :: FLX_IR_dn = 0.0_rp ! surface downwad radiation flux [J/m2/s]
894  real(RP) :: FLX_NIR_dn = 0.0_rp ! surface downwad radiation flux [J/m2/s]
895  real(RP) :: FLX_VIS_dn = 0.0_rp ! surface downwad radiation flux [J/m2/s]
896 
897  namelist / param_mkinit_flux / &
898  flx_rain, &
899  flx_snow, &
900  flx_ir_dn, &
901  flx_nir_dn, &
902  flx_vis_dn
903 
904  integer :: i, j
905  integer :: ierr
906  !---------------------------------------------------------------------------
907 
908  !--- read namelist
909  rewind(io_fid_conf)
910  read(io_fid_conf,nml=param_mkinit_flux,iostat=ierr)
911  if( ierr < 0 ) then !--- missing
912  log_info("flux_setup",*) 'Not found namelist. Default used.'
913  elseif( ierr > 0 ) then !--- fatal error
914  log_error("flux_setup",*) 'Not appropriate names in namelist PARAM_MKINIT_FLUX. Check!'
915  call prc_abort
916  endif
917  log_nml(param_mkinit_flux)
918 
919  do j = jsb, jeb
920  do i = isb, ieb
921  sflx_rain(i,j) = flx_rain
922  sflx_snow(i,j) = flx_snow
923 
924  sflx_lw_up(i,j) = 0.0_rp
925  sflx_lw_dn(i,j) = flx_ir_dn
926  sflx_sw_up(i,j) = 0.0_rp
927  sflx_sw_dn(i,j) = flx_nir_dn + flx_vis_dn
928  enddo
929  enddo
930 
931  return

References mod_atmos_phy_mp_vars::atmos_phy_mp_sflx_rain, mod_atmos_phy_mp_vars::atmos_phy_mp_sflx_snow, mod_atmos_phy_rd_vars::atmos_phy_rd_sflx_lw_dn, mod_atmos_phy_rd_vars::atmos_phy_rd_sflx_lw_up, mod_atmos_phy_rd_vars::atmos_phy_rd_sflx_sw_dn, mod_atmos_phy_rd_vars::atmos_phy_rd_sflx_sw_up, scale_atmos_grid_cartesc_index::ieb, scale_io::io_fid_conf, scale_atmos_grid_cartesc_index::isb, scale_atmos_grid_cartesc_index::jeb, scale_atmos_grid_cartesc_index::jsb, and scale_prc::prc_abort().

Referenced by read_sounding().

Here is the call graph for this function:
Here is the caller graph for this function:

◆ land_setup()

subroutine mod_mkinit::land_setup

Land setup.

Definition at line 937 of file mod_mkinit.F90.

937  use mod_land_vars, only: &
938  land_temp, &
939  land_water, &
940  land_ice, &
941  land_sfc_temp, &
942  land_sfc_albedo, &
943  snow_flag, &
944  snow_sfc_temp, &
945  snow_swe, &
946  snow_depth, &
947  snow_dzero, &
948  snow_nosnowsec
949  implicit none
950 
951  real(RP) :: LND_TEMP ! land soil temperature [K]
952  real(RP) :: LND_WATER = 0.15_rp ! land soil moisture [m3/m3]
953  real(RP) :: LND_ICE = 0.00_rp ! land soil ice [m3/m3]
954  real(RP) :: SFC_TEMP ! land skin temperature [K]
955  real(RP) :: SFC_albedo_LW = 0.01_rp ! land surface albedo for LW (0-1)
956  real(RP) :: SFC_albedo_SW = 0.20_rp ! land surface albedo for SW (0-1)
957 
958  namelist / param_mkinit_land / &
959  lnd_temp, &
960  lnd_water, &
961  lnd_ice, &
962  sfc_temp, &
963  sfc_albedo_lw, &
964  sfc_albedo_sw
965 
966  integer :: ierr
967  !---------------------------------------------------------------------------
968 
969  lnd_temp = thetastd
970  sfc_temp = thetastd
971 
972  !--- read namelist
973  rewind(io_fid_conf)
974  read(io_fid_conf,nml=param_mkinit_land,iostat=ierr)
975  if( ierr < 0 ) then !--- missing
976  log_info("land_setup",*) 'Not found namelist. Default used.'
977  elseif( ierr > 0 ) then !--- fatal error
978  log_error("land_setup",*) 'Not appropriate names in namelist PARAM_MKINIT_LAND. Check!'
979  call prc_abort
980  endif
981  log_nml(param_mkinit_land)
982 
983  land_temp(:,:,:) = lnd_temp
984  land_water(:,:,:) = lnd_water
985  land_ice(:,:,:) = lnd_ice
986 
987  land_sfc_temp(:,:) = sfc_temp
988  land_sfc_albedo(:,:,:,i_r_ir) = sfc_albedo_lw
989  land_sfc_albedo(:,:,:,i_r_nir) = sfc_albedo_sw
990  land_sfc_albedo(:,:,:,i_r_vis) = sfc_albedo_sw
991 
992  if ( snow_flag ) then
993  !!!!! Tentative for snow model !!!!!
994  snow_sfc_temp(:,:) = 273.15_rp
995  snow_swe(:,:) = 0.0_rp
996  snow_depth(:,:) = 0.0_rp
997  snow_dzero(:,:) = 0.0_rp
998  snow_nosnowsec(:,:) = 0.0_rp
999  end if
1000 
1001  return

References scale_cpl_sfc_index::i_r_ir, scale_cpl_sfc_index::i_r_nir, scale_cpl_sfc_index::i_r_vis, scale_io::io_fid_conf, mod_land_vars::land_ice, mod_land_vars::land_sfc_albedo, mod_land_vars::land_sfc_temp, mod_land_vars::land_temp, mod_land_vars::land_water, scale_prc::prc_abort(), mod_land_vars::snow_depth, mod_land_vars::snow_dzero, mod_land_vars::snow_flag, mod_land_vars::snow_nosnowsec, mod_land_vars::snow_sfc_temp, and mod_land_vars::snow_swe.

Referenced by read_sounding().

Here is the call graph for this function:
Here is the caller graph for this function:

◆ ocean_setup()

subroutine mod_mkinit::ocean_setup

Ocean setup.

Definition at line 1007 of file mod_mkinit.F90.

1007  use mod_ocean_vars, only: &
1008  ice_flag, &
1009  ocean_temp, &
1010  ocean_salt, &
1011  ocean_uvel, &
1012  ocean_vvel, &
1013  ocean_ocn_z0m, &
1014  ocean_ice_temp, &
1015  ocean_ice_mass, &
1016  ocean_sfc_temp, &
1017  ocean_sfc_albedo, &
1018  ocean_sfc_z0m, &
1019  ocean_sfc_z0h, &
1020  ocean_sfc_z0e
1021  implicit none
1022 
1023  real(RP) :: OCN_TEMP ! ocean temperature [K]
1024  real(RP) :: OCN_SALT = 0.0_rp ! ocean salinity [psu]
1025  real(RP) :: OCN_UVEL = 0.0_rp ! ocean u-velocity [m/s]
1026  real(RP) :: OCN_VVEL = 0.0_rp ! ocean v-velocity [m/s]
1027  real(RP) :: ICE_TEMP ! ocean temperature [K]
1028  real(RP) :: ICE_MASS = 0.0_rp ! ocean temperature [K]
1029  real(RP) :: SFC_TEMP ! ocean skin temperature [K]
1030  real(RP) :: SFC_albedo_LW = 0.04_rp ! ocean surface albedo for LW (0-1)
1031  real(RP) :: SFC_albedo_SW = 0.05_rp ! ocean surface albedo for SW (0-1)
1032  real(RP) :: SFC_Z0M = 1.e-4_rp ! ocean surface roughness length (momentum) [m]
1033  real(RP) :: SFC_Z0H = 1.e-4_rp ! ocean surface roughness length (heat) [m]
1034  real(RP) :: SFC_Z0E = 1.e-4_rp ! ocean surface roughness length (vapor) [m]
1035 
1036  namelist / param_mkinit_ocean / &
1037  ocn_temp, &
1038  ocn_salt, &
1039  ocn_uvel, &
1040  ocn_vvel, &
1041  ice_temp, &
1042  ice_mass, &
1043  sfc_temp, &
1044  sfc_albedo_lw, &
1045  sfc_albedo_sw, &
1046  sfc_z0m, &
1047  sfc_z0h, &
1048  sfc_z0e
1049 
1050  integer :: ierr
1051  !---------------------------------------------------------------------------
1052 
1053  ocn_temp = thetastd
1054  ice_temp = 271.35_rp ! freezing point of the ocean
1055  sfc_temp = thetastd
1056 
1057  !--- read namelist
1058  rewind(io_fid_conf)
1059  read(io_fid_conf,nml=param_mkinit_ocean,iostat=ierr)
1060  if( ierr < 0 ) then !--- missing
1061  log_info("ocean_setup",*) 'Not found namelist. Default used.'
1062  elseif( ierr > 0 ) then !--- fatal error
1063  log_error("ocean_setup",*) 'Not appropriate names in namelist PARAM_MKINIT_OCEAN. Check!'
1064  call prc_abort
1065  endif
1066  log_nml(param_mkinit_ocean)
1067 
1068  ocean_temp(:,:,:) = ocn_temp
1069  ocean_salt(:,:,:) = ocn_salt
1070  ocean_uvel(:,:,:) = ocn_uvel
1071  ocean_vvel(:,:,:) = ocn_vvel
1072  ocean_ocn_z0m(:,:) = sfc_z0m
1073  ocean_sfc_temp(:,:) = sfc_temp
1074  ocean_sfc_albedo(:,:,:,i_r_ir) = sfc_albedo_lw
1075  ocean_sfc_albedo(:,:,:,i_r_nir) = sfc_albedo_sw
1076  ocean_sfc_albedo(:,:,:,i_r_vis) = sfc_albedo_sw
1077  ocean_sfc_z0m(:,:) = sfc_z0m
1078  ocean_sfc_z0h(:,:) = sfc_z0h
1079  ocean_sfc_z0e(:,:) = sfc_z0e
1080 
1081  if ( ice_flag ) then
1082  ocean_ice_temp(:,:) = ice_temp
1083  ocean_ice_mass(:,:) = ice_mass
1084  end if
1085 
1086  return

References scale_cpl_sfc_index::i_r_ir, scale_cpl_sfc_index::i_r_nir, scale_cpl_sfc_index::i_r_vis, mod_ocean_vars::ice_flag, scale_io::io_fid_conf, mod_ocean_vars::ocean_ice_mass, mod_ocean_vars::ocean_ice_temp, mod_ocean_vars::ocean_ocn_z0m, mod_ocean_vars::ocean_salt, mod_ocean_vars::ocean_sfc_albedo, mod_ocean_vars::ocean_sfc_temp, mod_ocean_vars::ocean_sfc_z0e, mod_ocean_vars::ocean_sfc_z0h, mod_ocean_vars::ocean_sfc_z0m, mod_ocean_vars::ocean_temp, mod_ocean_vars::ocean_uvel, mod_ocean_vars::ocean_vvel, and scale_prc::prc_abort().

Referenced by read_sounding().

Here is the call graph for this function:
Here is the caller graph for this function:

◆ urban_setup()

subroutine mod_mkinit::urban_setup

Urban setup.

Definition at line 1092 of file mod_mkinit.F90.

1092  use mod_urban_vars, only: &
1093  urban_tr, &
1094  urban_tb, &
1095  urban_tg, &
1096  urban_tc, &
1097  urban_qc, &
1098  urban_uc, &
1099  urban_trl, &
1100  urban_tbl, &
1101  urban_tgl, &
1102  urban_rainr, &
1103  urban_rainb, &
1104  urban_raing, &
1105  urban_roff, &
1106  urban_sfc_temp, &
1107  urban_sfc_albedo
1108  implicit none
1109 
1110  real(RP) :: URB_ROOF_TEMP ! Surface temperature of roof [K]
1111  real(RP) :: URB_BLDG_TEMP ! Surface temperature of building [K]
1112  real(RP) :: URB_GRND_TEMP ! Surface temperature of ground [K]
1113  real(RP) :: URB_CNPY_TEMP ! Diagnostic canopy air temperature [K]
1114  real(RP) :: URB_CNPY_HMDT = 0.0_rp ! Diagnostic canopy humidity [kg/kg]
1115  real(RP) :: URB_CNPY_WIND = 0.0_rp ! Diagnostic canopy wind [m/s]
1116  real(RP) :: URB_ROOF_LAYER_TEMP ! temperature in layer of roof [K]
1117  real(RP) :: URB_BLDG_LAYER_TEMP ! temperature in layer of building [K]
1118  real(RP) :: URB_GRND_LAYER_TEMP ! temperature in layer of ground [K]
1119  real(RP) :: URB_ROOF_RAIN = 0.0_rp ! temperature in layer of roof [kg/m2]
1120  real(RP) :: URB_BLDG_RAIN = 0.0_rp ! temperature in layer of building [kg/m2]
1121  real(RP) :: URB_GRND_RAIN = 0.0_rp ! temperature in layer of ground [kg/m2]
1122  real(RP) :: URB_SFC_TEMP ! Grid average of surface temperature [K]
1123  real(RP) :: URB_ALB_LW = 0.10_rp ! Grid average of surface albedo for LW (0-1)
1124  real(RP) :: URB_ALB_SW = 0.20_rp ! Grid average of surface albedo for SW (0-1)
1125 
1126  namelist / param_mkinit_urban / &
1127  urb_roof_temp, &
1128  urb_bldg_temp, &
1129  urb_grnd_temp, &
1130  urb_cnpy_temp, &
1131  urb_cnpy_hmdt, &
1132  urb_cnpy_wind, &
1133  urb_roof_layer_temp, &
1134  urb_bldg_layer_temp, &
1135  urb_grnd_layer_temp, &
1136  urb_roof_rain, &
1137  urb_bldg_rain, &
1138  urb_grnd_rain, &
1139  urb_sfc_temp, &
1140  urb_alb_lw, &
1141  urb_alb_sw
1142 
1143  integer :: ierr
1144  !---------------------------------------------------------------------------
1145 
1146  urb_roof_temp = thetastd
1147  urb_bldg_temp = thetastd
1148  urb_grnd_temp = thetastd
1149  urb_cnpy_temp = thetastd
1150  urb_roof_layer_temp = thetastd
1151  urb_bldg_layer_temp = thetastd
1152  urb_grnd_layer_temp = thetastd
1153  urb_sfc_temp = thetastd
1154 
1155  !--- read namelist
1156  rewind(io_fid_conf)
1157  read(io_fid_conf,nml=param_mkinit_urban,iostat=ierr)
1158  if( ierr < 0 ) then !--- missing
1159  log_info("urban_setup",*) 'Not found namelist. Default used.'
1160  elseif( ierr > 0 ) then !--- fatal error
1161  log_error("urban_setup",*) 'Not appropriate names in namelist PARAM_MKINIT_URBAN. Check!'
1162  call prc_abort
1163  endif
1164  log_nml(param_mkinit_urban)
1165 
1166  urban_trl(:,:,:) = urb_roof_layer_temp
1167  urban_tbl(:,:,:) = urb_bldg_layer_temp
1168  urban_tgl(:,:,:) = urb_grnd_layer_temp
1169 
1170  urban_tr(:,:) = urb_roof_temp
1171  urban_tb(:,:) = urb_bldg_temp
1172  urban_tg(:,:) = urb_grnd_temp
1173  urban_tc(:,:) = urb_cnpy_temp
1174  urban_qc(:,:) = urb_cnpy_hmdt
1175  urban_uc(:,:) = urb_cnpy_wind
1176  urban_rainr(:,:) = urb_roof_rain
1177  urban_rainb(:,:) = urb_bldg_rain
1178  urban_raing(:,:) = urb_grnd_rain
1179  urban_sfc_temp(:,:) = urb_sfc_temp
1180  urban_sfc_albedo(:,:,:,i_r_ir) = urb_alb_lw
1181  urban_sfc_albedo(:,:,:,i_r_nir) = urb_alb_sw
1182  urban_sfc_albedo(:,:,:,i_r_vis) = urb_alb_sw
1183 
1184  return

References scale_cpl_sfc_index::i_r_ir, scale_cpl_sfc_index::i_r_nir, scale_cpl_sfc_index::i_r_vis, scale_io::io_fid_conf, scale_prc::prc_abort(), mod_urban_vars::urban_qc, mod_urban_vars::urban_rainb, mod_urban_vars::urban_raing, mod_urban_vars::urban_rainr, mod_urban_vars::urban_roff, mod_urban_vars::urban_sfc_albedo, mod_urban_vars::urban_sfc_temp, mod_urban_vars::urban_tb, mod_urban_vars::urban_tbl, mod_urban_vars::urban_tc, mod_urban_vars::urban_tg, mod_urban_vars::urban_tgl, mod_urban_vars::urban_tr, mod_urban_vars::urban_trl, and mod_urban_vars::urban_uc.

Referenced by read_sounding().

Here is the call graph for this function:
Here is the caller graph for this function:

◆ tke_setup()

subroutine mod_mkinit::tke_setup

TKE setup.

Definition at line 1190 of file mod_mkinit.F90.

1190  use scale_const, only: &
1191  eps => const_eps
1192  use mod_atmos_phy_tb_vars, only: &
1193  i_tke
1194  use mod_atmos_phy_bl_vars, only: &
1195  qs_bl => qs, &
1196  qe_bl => qe, &
1197  zi => atmos_phy_bl_zi
1198  implicit none
1199 
1200  real(RP) :: TKE_CONST
1201  real(RP) :: Zi_CONST
1202 
1203  namelist / param_mkinit_tke / &
1204  tke_const, &
1205  zi_const
1206 
1207  integer :: k, i, j
1208  integer :: ierr
1209  !---------------------------------------------------------------------------
1210 
1211  tke_const = eps
1212  zi_const = 100.0_rp
1213 
1214  !--- read namelist
1215  rewind(io_fid_conf)
1216  read(io_fid_conf,nml=param_mkinit_tke,iostat=ierr)
1217  if( ierr < 0 ) then !--- missing
1218  log_info("tke_setup",*) 'Not found namelist. Default used.'
1219  elseif( ierr > 0 ) then !--- fatal error
1220  log_error("tke_setup",*) 'Not appropriate names in namelist PARAM_MKINIT_TKE. Check!'
1221  call prc_abort
1222  endif
1223  log_nml(param_mkinit_tke)
1224 
1225  if ( i_tke > 0 ) then
1226  do j = 1, ja
1227  do i = 1, ia
1228  do k = 1, ka
1229  qtrc(k,i,j,i_tke) = tke_const
1230  enddo
1231  enddo
1232  enddo
1233  end if
1234  if ( qs_bl > 0 ) then
1235  do j = 1, ja
1236  do i = 1, ia
1237  do k = 1, ka
1238  qtrc(k,i,j,qs_bl) = tke_const
1239  qtrc(k,i,j,qs_bl+1:qe_bl) = 0.0_rp
1240  enddo
1241  enddo
1242  enddo
1243  end if
1244 
1245  do j = 1, ja
1246  do i = 1, ia
1247  zi(i,j) = zi_const
1248  end do
1249  end do
1250 
1251  return

References mod_atmos_phy_bl_vars::atmos_phy_bl_zi, scale_const::const_eps, mod_atmos_phy_tb_vars::i_tke, scale_atmos_grid_cartesc_index::ia, scale_io::io_fid_conf, scale_atmos_grid_cartesc_index::ja, scale_atmos_grid_cartesc_index::ka, scale_prc::prc_abort(), mod_atmos_phy_bl_vars::qe, mod_atmos_phy_bl_vars::qs, and mod_atmos_vars::qtrc.

Referenced by mkinit().

Here is the call graph for this function:
Here is the caller graph for this function:

◆ read_sounding()

subroutine mod_mkinit::read_sounding ( real(rp), dimension(ka), intent(out)  DENS,
real(rp), dimension(ka), intent(out)  VELX,
real(rp), dimension(ka), intent(out)  VELY,
real(rp), dimension(ka), intent(out)  POTT,
real(rp), dimension (ka), intent(out)  QV 
)

Read sounding data from file.

Definition at line 1258 of file mod_mkinit.F90.

1258  use scale_atmos_hydrometeor, only: &
1259  atmos_hydrometeor_dry
1260  implicit none
1261 
1262  real(RP), intent(out) :: DENS(KA)
1263  real(RP), intent(out) :: VELX(KA)
1264  real(RP), intent(out) :: VELY(KA)
1265  real(RP), intent(out) :: POTT(KA)
1266  real(RP), intent(out) :: QV (KA)
1267 
1268  real(RP) :: TEMP(KA)
1269  real(RP) :: PRES(KA)
1270  real(RP) :: QC (KA)
1271 
1272  character(len=H_LONG) :: ENV_IN_SOUNDING_file = ''
1273 
1274  integer, parameter :: EXP_klim = 100
1275  integer :: EXP_kmax
1276 
1277  real(RP) :: SFC_THETA ! surface potential temperature [K]
1278  real(RP) :: SFC_PRES ! surface pressure [hPa]
1279  real(RP) :: SFC_QV ! surface watervapor [g/kg]
1280 
1281  real(RP) :: EXP_z (EXP_klim+1) ! height [m]
1282  real(RP) :: EXP_pott(EXP_klim+1) ! potential temperature [K]
1283  real(RP) :: EXP_qv (EXP_klim+1) ! water vapor [g/kg]
1284  real(RP) :: EXP_u (EXP_klim+1) ! velocity u [m/s]
1285  real(RP) :: EXP_v (EXP_klim+1) ! velocity v [m/s]
1286 
1287  real(RP) :: fact1, fact2
1288  integer :: k, kref
1289  integer :: fid
1290  integer :: ierr
1291 
1292  namelist / param_mkinit_sounding / &
1293  env_in_sounding_file
1294 
1295  !--- read namelist
1296  rewind(io_fid_conf)
1297  read(io_fid_conf,nml=param_mkinit_sounding,iostat=ierr)
1298 
1299  if( ierr < 0 ) then !--- missing
1300  log_info("read_sounding",*) 'Not found namelist. Default used.'
1301  elseif( ierr > 0 ) then !--- fatal error
1302  log_error("read_sounding",*) 'Not appropriate names in namelist PARAM_MKINIT_SOUNDING. Check!'
1303  call prc_abort
1304  endif
1305  log_nml(param_mkinit_sounding)
1306 
1307  !--- prepare sounding profile
1308  log_info("read_sounding",*) 'Input sounding file:', trim(env_in_sounding_file)
1309  fid = io_get_available_fid()
1310  open( fid, &
1311  file = trim(env_in_sounding_file), &
1312  form = 'formatted', &
1313  status = 'old', &
1314  iostat = ierr )
1315 
1316  if ( ierr /= 0 ) then
1317  log_error("read_sounding",*) '[mod_mkinit/read_sounding] Input file not found!'
1318  endif
1319 
1320  !--- read sounding file till end
1321  read(fid,*) sfc_pres, sfc_theta, sfc_qv
1322 
1323  log_info("read_sounding",*) '+ Surface pressure [hPa]', sfc_pres
1324  log_info("read_sounding",*) '+ Surface pot. temp [K]', sfc_theta
1325  log_info("read_sounding",*) '+ Surface water vapor [g/kg]', sfc_qv
1326 
1327  do k = 2, exp_klim
1328  read(fid,*,iostat=ierr) exp_z(k), exp_pott(k), exp_qv(k), exp_u(k), exp_v(k)
1329  if ( ierr /= 0 ) exit
1330  enddo
1331 
1332  exp_kmax = k - 1
1333  close(fid)
1334 
1335  ! Boundary
1336  exp_z(1) = 0.0_rp
1337  exp_pott(1) = sfc_theta
1338  exp_qv(1) = sfc_qv
1339  exp_u(1) = exp_u(2)
1340  exp_v(1) = exp_v(2)
1341  exp_z(exp_kmax+1) = 100.e3_rp
1342  exp_pott(exp_kmax+1) = exp_pott(exp_kmax)
1343  exp_qv(exp_kmax+1) = exp_qv(exp_kmax)
1344  exp_u(exp_kmax+1) = exp_u(exp_kmax)
1345  exp_v(exp_kmax+1) = exp_v(exp_kmax)
1346 
1347  do k = 1, exp_kmax+1
1348  exp_qv(k) = exp_qv(k) * 1.e-3_rp ! [g/kg]->[kg/kg]
1349  enddo
1350 
1351  ! calc in dry condition
1352  pres_sfc(:,:) = sfc_pres * 1.e2_rp ! [hPa]->[Pa]
1353  pott_sfc(:,:) = sfc_theta
1354  if ( .not. atmos_hydrometeor_dry ) then
1355  qv_sfc(:,:) = sfc_qv * 1.e-3_rp ! [g/kg]->[kg/kg]
1356  end if
1357 
1358  !--- linear interpolate to model grid
1359  do k = ks, ke
1360  do kref = 2, exp_kmax+1
1361  if ( cz(k) > exp_z(kref-1) &
1362  .AND. cz(k) <= exp_z(kref ) ) then
1363 
1364  fact1 = ( exp_z(kref) - cz(k) ) / ( exp_z(kref)-exp_z(kref-1) )
1365  fact2 = ( cz(k) - exp_z(kref-1) ) / ( exp_z(kref)-exp_z(kref-1) )
1366 
1367  pott(k) = exp_pott(kref-1) * fact1 &
1368  + exp_pott(kref ) * fact2
1369  qv(k) = exp_qv(kref-1) * fact1 &
1370  + exp_qv(kref ) * fact2
1371  velx(k) = exp_u(kref-1) * fact1 &
1372  + exp_u(kref ) * fact2
1373  vely(k) = exp_v(kref-1) * fact1 &
1374  + exp_v(kref ) * fact2
1375  endif
1376  enddo
1377  enddo
1378  if ( atmos_hydrometeor_dry ) qv(:) = 0.0_rp
1379 
1380  qc(:) = 0.0_rp
1381 
1382  ! make density & pressure profile in moist condition
1383  call hydrostatic_buildrho( ka, ks, ke, &
1384  pott(:), qv(:), qc(:), & ! [IN]
1385  pres_sfc(1,1), pott_sfc(1,1), qv_sfc(1,1), qc_sfc(1,1), & ! [IN]
1386  cz(:), fz(:), & ! [IN]
1387  dens(:), temp(:), pres(:), temp_sfc(1,1) ) ! [OUT]
1388 
1389  return

References scale_atmos_grid_cartesc::atmos_grid_cartesc_domain_center_x, scale_atmos_grid_cartesc::atmos_grid_cartesc_domain_center_y, scale_atmos_grid_cartesc::atmos_grid_cartesc_fyg, scale_atmos_hydrometeor::atmos_hydrometeor_dry, mod_atmos_admin::atmos_phy_ae_type, scale_const::const_cl, scale_const::const_cpdry, scale_const::const_cpvap, scale_const::const_cvdry, scale_const::const_cvvap, scale_const::const_grav, scale_const::const_ohm, scale_const::const_radius, scale_const::const_rdry, scale_const::const_rvap, mod_atmos_vars::dens, flux_setup(), scale_atmos_grid_cartesc_index::ia, scale_atmos_grid_cartesc_index::ie, scale_atmos_grid_cartesc_index::ieb, scale_atmos_grid_cartesc_index::imax, scale_io::io_fid_conf, scale_io::io_get_available_fid(), scale_atmos_grid_cartesc_index::is, scale_atmos_grid_cartesc_index::isb, scale_atmos_grid_cartesc_index::ja, scale_atmos_grid_cartesc_index::jag, scale_atmos_grid_cartesc_index::jeb, scale_atmos_grid_cartesc_index::jhalo, scale_atmos_grid_cartesc_index::js, scale_atmos_grid_cartesc_index::jsb, scale_atmos_grid_cartesc_index::ka, scale_atmos_grid_cartesc_index::ke, scale_atmos_grid_cartesc_index::ks, scale_land_grid_cartesc_real::land_grid_cartesc_real_set_areavol(), land_setup(), scale_landuse::landuse_calc_fact(), scale_landuse::landuse_fillhalo(), scale_landuse::landuse_frac_land, scale_landuse::landuse_frac_urban, mod_atmos_vars::momx, mod_atmos_vars::momy, mod_atmos_vars::momz, scale_ocean_grid_cartesc_real::ocean_grid_cartesc_real_set_areavol(), ocean_setup(), scale_prc::prc_abort(), scale_prc_cartesc::prc_num_x, scale_prof::prof_rapend(), scale_prof::prof_rapstart(), mod_realinput::realinput_atmos(), mod_realinput::realinput_surface(), rect_setup(), mod_atmos_vars::rhot, and urban_setup().

Here is the call graph for this function:

Variable Documentation

◆ mkinit_type

integer, public mod_mkinit::mkinit_type = -1

Definition at line 88 of file mod_mkinit.F90.

88  integer, public :: MKINIT_TYPE = -1

Referenced by mkinit(), and mkinit_setup().

◆ i_ignore

integer, parameter, public mod_mkinit::i_ignore = 0

Definition at line 89 of file mod_mkinit.F90.

89  integer, public, parameter :: I_IGNORE = 0

Referenced by mkinit(), and mkinit_setup().

◆ i_planestate

integer, parameter, public mod_mkinit::i_planestate = 1

Definition at line 91 of file mod_mkinit.F90.

91  integer, public, parameter :: I_PLANESTATE = 1

Referenced by mkinit(), and mkinit_setup().

◆ i_tracerbubble

integer, parameter, public mod_mkinit::i_tracerbubble = 2

Definition at line 92 of file mod_mkinit.F90.

92  integer, public, parameter :: I_TRACERBUBBLE = 2

Referenced by mkinit(), and mkinit_setup().

◆ i_coldbubble

integer, parameter, public mod_mkinit::i_coldbubble = 3

Definition at line 93 of file mod_mkinit.F90.

93  integer, public, parameter :: I_COLDBUBBLE = 3

Referenced by mkinit(), and mkinit_setup().

◆ i_lambwave

integer, parameter, public mod_mkinit::i_lambwave = 4

Definition at line 95 of file mod_mkinit.F90.

95  integer, public, parameter :: I_LAMBWAVE = 4

Referenced by mkinit(), and mkinit_setup().

◆ i_gravitywave

integer, parameter, public mod_mkinit::i_gravitywave = 5

Definition at line 96 of file mod_mkinit.F90.

96  integer, public, parameter :: I_GRAVITYWAVE = 5

Referenced by mkinit(), and mkinit_setup().

◆ i_khwave

integer, parameter, public mod_mkinit::i_khwave = 6

Definition at line 97 of file mod_mkinit.F90.

97  integer, public, parameter :: I_KHWAVE = 6

Referenced by mkinit(), and mkinit_setup().

◆ i_turbulence

integer, parameter, public mod_mkinit::i_turbulence = 7

Definition at line 98 of file mod_mkinit.F90.

98  integer, public, parameter :: I_TURBULENCE = 7

Referenced by mkinit(), and mkinit_setup().

◆ i_mountainwave

integer, parameter, public mod_mkinit::i_mountainwave = 8

Definition at line 99 of file mod_mkinit.F90.

99  integer, public, parameter :: I_MOUNTAINWAVE = 8

Referenced by mkinit(), and mkinit_setup().

◆ i_warmbubble

integer, parameter, public mod_mkinit::i_warmbubble = 9

Definition at line 101 of file mod_mkinit.F90.

101  integer, public, parameter :: I_WARMBUBBLE = 9

Referenced by mkinit(), and mkinit_setup().

◆ i_supercell

integer, parameter, public mod_mkinit::i_supercell = 10

Definition at line 102 of file mod_mkinit.F90.

102  integer, public, parameter :: I_SUPERCELL = 10

Referenced by mkinit(), and mkinit_setup().

◆ i_squallline

integer, parameter, public mod_mkinit::i_squallline = 11

Definition at line 103 of file mod_mkinit.F90.

103  integer, public, parameter :: I_SQUALLLINE = 11

Referenced by mkinit(), and mkinit_setup().

◆ i_wk1982

integer, parameter, public mod_mkinit::i_wk1982 = 12

Definition at line 104 of file mod_mkinit.F90.

104  integer, public, parameter :: I_WK1982 = 12

Referenced by mkinit(), and mkinit_setup().

◆ i_dycoms2_rf01

integer, parameter, public mod_mkinit::i_dycoms2_rf01 = 13

Definition at line 105 of file mod_mkinit.F90.

105  integer, public, parameter :: I_DYCOMS2_RF01 = 13

Referenced by mkinit(), and mkinit_setup().

◆ i_dycoms2_rf02

integer, parameter, public mod_mkinit::i_dycoms2_rf02 = 14

Definition at line 106 of file mod_mkinit.F90.

106  integer, public, parameter :: I_DYCOMS2_RF02 = 14

Referenced by mkinit(), and mkinit_setup().

◆ i_rico

integer, parameter, public mod_mkinit::i_rico = 15

Definition at line 107 of file mod_mkinit.F90.

107  integer, public, parameter :: I_RICO = 15

Referenced by mkinit(), and mkinit_setup().

◆ i_interporation

integer, parameter, public mod_mkinit::i_interporation = 16

Definition at line 109 of file mod_mkinit.F90.

109  integer, public, parameter :: I_INTERPORATION = 16

Referenced by mkinit_setup().

◆ i_landcouple

integer, parameter, public mod_mkinit::i_landcouple = 17

Definition at line 111 of file mod_mkinit.F90.

111  integer, public, parameter :: I_LANDCOUPLE = 17

Referenced by mkinit(), and mkinit_setup().

◆ i_oceancouple

integer, parameter, public mod_mkinit::i_oceancouple = 18

Definition at line 112 of file mod_mkinit.F90.

112  integer, public, parameter :: I_OCEANCOUPLE = 18

Referenced by mkinit(), and mkinit_setup().

◆ i_urbancouple

integer, parameter, public mod_mkinit::i_urbancouple = 19

Definition at line 113 of file mod_mkinit.F90.

113  integer, public, parameter :: I_URBANCOUPLE = 19

Referenced by mkinit(), and mkinit_setup().

◆ i_triplecouple

integer, parameter, public mod_mkinit::i_triplecouple = 20

Definition at line 114 of file mod_mkinit.F90.

114  integer, public, parameter :: I_TRIPLECOUPLE = 20

Referenced by mkinit(), and mkinit_setup().

◆ i_bubblecouple

integer, parameter, public mod_mkinit::i_bubblecouple = 21

Definition at line 115 of file mod_mkinit.F90.

115  integer, public, parameter :: I_BUBBLECOUPLE = 21

Referenced by mkinit(), and mkinit_setup().

◆ i_seabreeze

integer, parameter, public mod_mkinit::i_seabreeze = 22

Definition at line 117 of file mod_mkinit.F90.

117  integer, public, parameter :: I_SEABREEZE = 22

Referenced by mkinit(), and mkinit_setup().

◆ i_heatisland

integer, parameter, public mod_mkinit::i_heatisland = 23

Definition at line 118 of file mod_mkinit.F90.

118  integer, public, parameter :: I_HEATISLAND = 23

Referenced by mkinit(), and mkinit_setup().

◆ i_dycoms2_rf02_dns

integer, parameter, public mod_mkinit::i_dycoms2_rf02_dns = 24

Definition at line 120 of file mod_mkinit.F90.

120  integer, public, parameter :: I_DYCOMS2_RF02_DNS = 24

Referenced by mkinit(), and mkinit_setup().

◆ i_real

integer, parameter, public mod_mkinit::i_real = 25

Definition at line 122 of file mod_mkinit.F90.

122  integer, public, parameter :: I_REAL = 25

Referenced by mkinit(), and mkinit_setup().

◆ i_grayzone

integer, parameter, public mod_mkinit::i_grayzone = 26

Definition at line 124 of file mod_mkinit.F90.

124  integer, public, parameter :: I_GRAYZONE = 26

Referenced by mkinit(), and mkinit_setup().

◆ i_boxaero

integer, parameter, public mod_mkinit::i_boxaero = 27

Definition at line 125 of file mod_mkinit.F90.

125  integer, public, parameter :: I_BOXAERO = 27

Referenced by mkinit(), and mkinit_setup().

◆ i_warmbubbleaero

integer, parameter, public mod_mkinit::i_warmbubbleaero = 28

Definition at line 126 of file mod_mkinit.F90.

126  integer, public, parameter :: I_WARMBUBBLEAERO = 28

Referenced by mkinit(), and mkinit_setup().

◆ i_cavityflow

integer, parameter, public mod_mkinit::i_cavityflow = 29

Definition at line 128 of file mod_mkinit.F90.

128  integer, public, parameter :: I_CAVITYFLOW = 29

Referenced by mkinit(), and mkinit_setup().

◆ i_barocwave

integer, parameter, public mod_mkinit::i_barocwave = 30

Definition at line 129 of file mod_mkinit.F90.

129  integer, public, parameter :: I_BAROCWAVE = 30

Referenced by mkinit(), and mkinit_setup().

◆ i_bomex

integer, parameter, public mod_mkinit::i_bomex = 31

Definition at line 130 of file mod_mkinit.F90.

130  integer, public, parameter :: I_BOMEX = 31

Referenced by mkinit(), and mkinit_setup().

mod_atmos_phy_mp_vars
module Atmosphere / Physics Cloud Microphysics
Definition: mod_atmos_phy_mp_vars.F90:12
mod_ocean_vars::ocean_ocn_z0m
real(rp), dimension(:,:), allocatable, public ocean_ocn_z0m
surface roughness length for momentum, open ocean [m]
Definition: mod_ocean_vars.F90:65
mod_land_vars::land_temp
real(rp), dimension(:,:,:), allocatable, public land_temp
temperature of each soil layer [K]
Definition: mod_land_vars.F90:61
mod_urban_vars::urban_sfc_temp
real(rp), dimension(:,:), allocatable, public urban_sfc_temp
Definition: mod_urban_vars.F90:74
scale_prc::prc_abort
subroutine, public prc_abort
Abort Process.
Definition: scale_prc.F90:342
mod_ocean_vars::ocean_ice_temp
real(rp), dimension(:,:), allocatable, public ocean_ice_temp
sea ice temperature [K]
Definition: mod_ocean_vars.F90:73
mod_land_vars::snow_sfc_temp
real(rp), dimension(:,:), allocatable, public snow_sfc_temp
snow surface temperature [K]
Definition: mod_land_vars.F90:68
mod_urban_vars::urban_tb
real(rp), dimension(:,:), allocatable, public urban_tb
Definition: mod_urban_vars.F90:64
mod_atmos_vars::pott
real(rp), dimension(:,:,:), allocatable, target, public pott
Definition: mod_atmos_vars.F90:132
mod_ocean_vars::ocean_sfc_z0e
real(rp), dimension(:,:), allocatable, public ocean_sfc_z0e
ocean surface roughness length for vapor [m]
Definition: mod_ocean_vars.F90:71
mod_land_vars::land_sfc_albedo
real(rp), dimension(:,:,:,:), allocatable, public land_sfc_albedo
land surface albedo (direct/diffuse,IR/near-IR/VIS) (0-1)
Definition: mod_land_vars.F90:65
scale_const::const_undef8
real(dp), parameter, public const_undef8
undefined value (REAL8)
Definition: scale_const.F90:40
mod_ocean_vars::ice_flag
logical, public ice_flag
Definition: mod_ocean_vars.F90:141
mod_atmos_admin
module ATMOS admin
Definition: mod_atmos_admin.F90:11
mod_urban_vars::urban_sfc_albedo
real(rp), dimension(:,:,:,:), allocatable, public urban_sfc_albedo
Definition: mod_urban_vars.F90:75
mod_atmos_phy_mp_vars::qs_mp
integer, public qs_mp
Definition: mod_atmos_phy_mp_vars.F90:78
mod_land_vars::land_water
real(rp), dimension(:,:,:), allocatable, public land_water
moisture of each soil layer [m3/m3]
Definition: mod_land_vars.F90:62
scale_const::const_eps
real(rp), public const_eps
small number
Definition: scale_const.F90:33
mod_urban_vars::urban_trl
real(rp), dimension(:,:,:), allocatable, public urban_trl
Definition: mod_urban_vars.F90:60
scale_atmos_hydrometeor
module atmosphere / hydrometeor
Definition: scale_atmos_hydrometeor.F90:12
mod_atmos_phy_rd_vars
module Atmosphere / Physics Radiation
Definition: mod_atmos_phy_rd_vars.F90:12
mod_atmos_phy_mp_vars::atmos_phy_mp_sflx_rain
real(rp), dimension(:,:), allocatable, public atmos_phy_mp_sflx_rain
Definition: mod_atmos_phy_mp_vars.F90:73
mod_urban_vars::urban_raing
real(rp), dimension(:,:), allocatable, public urban_raing
Definition: mod_urban_vars.F90:71
mod_urban_vars::urban_qc
real(rp), dimension(:,:), allocatable, public urban_qc
Definition: mod_urban_vars.F90:67
mod_ocean_vars::ocean_sfc_z0m
real(rp), dimension(:,:), allocatable, public ocean_sfc_z0m
ocean surface roughness length for momentum [m]
Definition: mod_ocean_vars.F90:69
scale_atmos_hydrometeor::atmos_hydrometeor_dry
logical, public atmos_hydrometeor_dry
Definition: scale_atmos_hydrometeor.F90:97
mod_atmos_vars::qtrc
real(rp), dimension(:,:,:,:), allocatable, target, public qtrc
Definition: mod_atmos_vars.F90:80
mod_atmos_phy_mp_driver
module atmosphere / physics / cloud microphysics
Definition: mod_atmos_phy_mp_driver.F90:12
scale_const::const_pi
real(rp), public const_pi
pi
Definition: scale_const.F90:31
mod_atmos_phy_rd_vars::atmos_phy_rd_sflx_lw_dn
real(rp), dimension(:,:), allocatable, public atmos_phy_rd_sflx_lw_dn
Definition: mod_atmos_phy_rd_vars.F90:61
mod_atmos_phy_bl_vars
module atmosphere / physics / PBL
Definition: mod_atmos_phy_bl_vars.F90:12
mod_atmos_phy_rd_vars::atmos_phy_rd_sflx_sw_dn
real(rp), dimension(:,:), allocatable, public atmos_phy_rd_sflx_sw_dn
Definition: mod_atmos_phy_rd_vars.F90:63
mod_atmos_vars::dens
real(rp), dimension(:,:,:), allocatable, target, public dens
Definition: mod_atmos_vars.F90:75
mod_ocean_vars::ocean_sfc_z0h
real(rp), dimension(:,:), allocatable, public ocean_sfc_z0h
ocean surface roughness length for heat [m]
Definition: mod_ocean_vars.F90:70
scale_const
module CONSTANT
Definition: scale_const.F90:11
mod_ocean_vars::ocean_uvel
real(rp), dimension(:,:,:), allocatable, public ocean_uvel
ocean zonal velocity [m/s]
Definition: mod_ocean_vars.F90:62
mod_atmos_phy_rd_vars::atmos_phy_rd_sflx_sw_up
real(rp), dimension(:,:), allocatable, public atmos_phy_rd_sflx_sw_up
Definition: mod_atmos_phy_rd_vars.F90:62
mod_atmos_phy_rd_vars::atmos_phy_rd_sflx_lw_up
real(rp), dimension(:,:), allocatable, public atmos_phy_rd_sflx_lw_up
Definition: mod_atmos_phy_rd_vars.F90:60
mod_urban_vars::urban_tgl
real(rp), dimension(:,:,:), allocatable, public urban_tgl
Definition: mod_urban_vars.F90:62
mod_atmos_phy_tb_vars
module Atmosphere / Physics Turbulence
Definition: mod_atmos_phy_tb_vars.F90:12
mod_land_vars
module LAND Variables
Definition: mod_land_vars.F90:11
mod_atmos_phy_bl_vars::atmos_phy_bl_zi
real(rp), dimension(:,:), allocatable, public atmos_phy_bl_zi
Definition: mod_atmos_phy_bl_vars.F90:64
mod_land_vars::snow_nosnowsec
real(rp), dimension(:,:), allocatable, public snow_nosnowsec
sec while no snow [s]
Definition: mod_land_vars.F90:72
mod_ocean_vars
module OCEAN Variables
Definition: mod_ocean_vars.F90:12
mod_land_vars::snow_swe
real(rp), dimension(:,:), allocatable, public snow_swe
snow water equivalent [kg/m2]
Definition: mod_land_vars.F90:69
scale_atmos_hydrometeor::i_hc
integer, parameter, public i_hc
liquid water cloud
Definition: scale_atmos_hydrometeor.F90:81
mod_atmos_phy_bl_vars::qe
integer, public qe
Definition: mod_atmos_phy_bl_vars.F90:45
mod_urban_vars::urban_tc
real(rp), dimension(:,:), allocatable, public urban_tc
Definition: mod_urban_vars.F90:66
mod_urban_vars::urban_tg
real(rp), dimension(:,:), allocatable, public urban_tg
Definition: mod_urban_vars.F90:65
mod_urban_vars::urban_uc
real(rp), dimension(:,:), allocatable, public urban_uc
Definition: mod_urban_vars.F90:68
mod_atmos_phy_tb_vars::i_tke
integer, public i_tke
Definition: mod_atmos_phy_tb_vars.F90:61
mod_land_vars::land_sfc_temp
real(rp), dimension(:,:), allocatable, public land_sfc_temp
land surface skin temperature [K]
Definition: mod_land_vars.F90:64
mod_ocean_vars::ocean_ice_mass
real(rp), dimension(:,:), allocatable, public ocean_ice_mass
sea ice mass [kg]
Definition: mod_ocean_vars.F90:74
mod_ocean_vars::ocean_vvel
real(rp), dimension(:,:,:), allocatable, public ocean_vvel
ocean meridional velocity [m/s]
Definition: mod_ocean_vars.F90:63
mod_land_vars::snow_depth
real(rp), dimension(:,:), allocatable, public snow_depth
snow depth [m]
Definition: mod_land_vars.F90:70
mod_urban_vars::urban_rainb
real(rp), dimension(:,:), allocatable, public urban_rainb
Definition: mod_urban_vars.F90:70
mod_land_vars::snow_flag
logical, public snow_flag
Definition: mod_land_vars.F90:132
mod_ocean_vars::ocean_sfc_albedo
real(rp), dimension(:,:,:,:), allocatable, public ocean_sfc_albedo
ocean surface albedo (direct/diffuse,IR/near-IR/VIS) (0-1)
Definition: mod_ocean_vars.F90:68
mod_atmos_phy_mp_vars::atmos_phy_mp_sflx_snow
real(rp), dimension(:,:), allocatable, public atmos_phy_mp_sflx_snow
Definition: mod_atmos_phy_mp_vars.F90:74
mod_urban_vars::urban_tr
real(rp), dimension(:,:), allocatable, public urban_tr
Definition: mod_urban_vars.F90:63
mod_ocean_vars::ocean_sfc_temp
real(rp), dimension(:,:), allocatable, public ocean_sfc_temp
ocean surface skin temperature [K]
Definition: mod_ocean_vars.F90:67
mod_atmos_admin::atmos_phy_mp_type
character(len=h_short), public atmos_phy_mp_type
Definition: mod_atmos_admin.F90:36
mod_urban_vars::urban_roff
real(rp), dimension(:,:), allocatable, public urban_roff
Definition: mod_urban_vars.F90:91
mod_atmos_phy_mp_vars::qe_mp
integer, public qe_mp
Definition: mod_atmos_phy_mp_vars.F90:79
mod_ocean_vars::ocean_salt
real(rp), dimension(:,:,:), allocatable, public ocean_salt
ocean salinity [PSU]
Definition: mod_ocean_vars.F90:61
mod_urban_vars::urban_tbl
real(rp), dimension(:,:,:), allocatable, public urban_tbl
Definition: mod_urban_vars.F90:61
mod_land_vars::land_ice
real(rp), dimension(:,:,:), allocatable, public land_ice
ice of each soil layer [m3/m3]
Definition: mod_land_vars.F90:63
mod_atmos_phy_mp_driver::atmos_phy_mp_driver_qhyd2qtrc
subroutine, public atmos_phy_mp_driver_qhyd2qtrc(KA, KS, KE, IA, IS, IE, JA, JS, JE, QV, QHYD, QTRC, QNUM)
Definition: mod_atmos_phy_mp_driver.F90:1361
mod_urban_vars::urban_rainr
real(rp), dimension(:,:), allocatable, public urban_rainr
Definition: mod_urban_vars.F90:69
scale_atmos_hydrometeor::n_hyd
integer, parameter, public n_hyd
Definition: scale_atmos_hydrometeor.F90:79
mod_urban_vars
module URBAN Variables
Definition: mod_urban_vars.F90:12
mod_land_vars::snow_dzero
real(rp), dimension(:,:), allocatable, public snow_dzero
snow depth at melting point [m]
Definition: mod_land_vars.F90:71
mod_ocean_vars::ocean_temp
real(rp), dimension(:,:,:), allocatable, public ocean_temp
ocean temperature [K]
Definition: mod_ocean_vars.F90:60
mod_atmos_phy_bl_vars::qs
integer, public qs
Definition: mod_atmos_phy_bl_vars.F90:45
Generated by   doxygen 1.8.17