SCALE-RM
Functions/Subroutines | Variables
scale_atmos_refstate Module Reference

module atmosphere / reference state More...

Functions/Subroutines

subroutine, public atmos_refstate_setup (KA, KS, KE, IA, IS, IE, JA, JS, JE, CZ, FZ, REAL_CZ, REAL_FZ, REAL_PHI)
 Setup. More...
 
subroutine, public atmos_refstate_read (KA, KS, KE, IA, IS, IE, JA, JS, JE, CZ, FZ, REAL_CZ, REAL_FZ, REAL_PHI)
 Read reference state profile. More...
 
subroutine, public atmos_refstate_write
 Write reference state profile. More...
 
subroutine, public atmos_refstate_update (KA, KS, KE, IA, IS, IE, ISB, IEB, JA, JS, JE, JSB, JEB, DENS, POTT, TEMP, PRES, QV, CZ, FZ, FDZ, RCDZ, REAL_CZ, REAL_FZ, REAL_PHI, AREA, nowsec, force)
 Update reference state profile (Horizontal average) More...
 
subroutine atmos_refstate_calc3d (KA, KS, KE, IA, IS, IE, JA, JS, JE, CZ, FZ, REAL_CZ, REAL_FZ, REAL_PHI)
 apply 1D reference to 3D (terrain-following) with re-calc hydrostatic balance More...
 
subroutine atmos_refstate_smoothing (KA, KS, KE, FDZ, RCDZ, phi)
 

Variables

logical, public atmos_refstate_update_flag = .false.
 
real(rp), dimension(:,:,:), allocatable, public atmos_refstate_pres
 refernce pressure [Pa] More...
 
real(rp), dimension(:,:,:), allocatable, public atmos_refstate_temp
 refernce temperature [K] More...
 
real(rp), dimension(:,:,:), allocatable, public atmos_refstate_dens
 refernce density [kg/m3] More...
 
real(rp), dimension(:,:,:), allocatable, public atmos_refstate_pott
 refernce potential temperature [K] More...
 
real(rp), dimension(:,:,:), allocatable, public atmos_refstate_qv
 refernce vapor [kg/kg] More...
 

Detailed Description

module atmosphere / reference state

Description
Reference state of Atmosphere
Author
Team SCALE
NAMELIST
  • PARAM_ATMOS_REFSTATE
    nametypedefault valuecomment
    ATMOS_REFSTATE_IN_BASENAME character(len=H_LONG) '' basename of the input file
    ATMOS_REFSTATE_OUT_BASENAME character(len=H_LONG) '' basename of the output file
    ATMOS_REFSTATE_OUT_TITLE character(len=H_MID) 'SCALE-RM RefState' title of the output file
    ATMOS_REFSTATE_OUT_DTYPE character(len=H_SHORT) 'DEFAULT' REAL4 or REAL8
    ATMOS_REFSTATE_TYPE character(len=H_SHORT) 'UNIFORM' profile type
    ATMOS_REFSTATE_TEMP_SFC real(RP) 300.0_RP surface temperature [K]
    ATMOS_REFSTATE_RH real(RP) 0.0_RP surface & environment RH [%]
    ATMOS_REFSTATE_POTT_UNIFORM real(RP) 300.0_RP uniform potential temperature [K]
    ATMOS_REFSTATE_UPDATE_DT real(DP) -1.0_DP

History Output
No history output

Function/Subroutine Documentation

◆ atmos_refstate_setup()

subroutine, public scale_atmos_refstate::atmos_refstate_setup ( integer, intent(in)  KA,
integer, intent(in)  KS,
integer, intent(in)  KE,
integer, intent(in)  IA,
integer, intent(in)  IS,
integer, intent(in)  IE,
integer, intent(in)  JA,
integer, intent(in)  JS,
integer, intent(in)  JE,
real(rp), dimension ( ka), intent(in)  CZ,
real(rp), dimension (0:ka), intent(in)  FZ,
real(rp), dimension ( ka,ia,ja), intent(in)  REAL_CZ,
real(rp), dimension (0:ka,ia,ja), intent(in)  REAL_FZ,
real(rp), dimension( ka,ia,ja), intent(in)  REAL_PHI 
)

Setup.

Definition at line 83 of file scale_atmos_refstate.F90.

83  use scale_const, only: &
84  undef => const_undef
85  use scale_prc, only: &
86  prc_abort
87  implicit none
88  integer, intent(in) :: KA, KS, KE
89  integer, intent(in) :: IA, IS, IE
90  integer, intent(in) :: JA, JS, JE
91 
92  real(RP), intent(in) :: CZ ( KA)
93  real(RP), intent(in) :: FZ (0:KA)
94  real(RP), intent(in) :: REAL_CZ ( KA,IA,JA)
95  real(RP), intent(in) :: REAL_FZ (0:KA,IA,JA)
96  real(RP), intent(in) :: REAL_PHI( KA,IA,JA)
97 
98  namelist / param_atmos_refstate / &
99  atmos_refstate_in_basename, &
100  atmos_refstate_out_basename, &
101  atmos_refstate_out_title, &
102  atmos_refstate_out_dtype, &
103  atmos_refstate_type, &
104  atmos_refstate_temp_sfc, &
105  atmos_refstate_rh, &
106  atmos_refstate_pott_uniform, &
107  atmos_refstate_update_dt
108 
109  integer :: ierr
110  !---------------------------------------------------------------------------
111 
112  log_newline
113  log_info("ATMOS_REFSTATE_setup",*) 'Setup'
114 
115  !--- read namelist
116  rewind(io_fid_conf)
117  read(io_fid_conf,nml=param_atmos_refstate,iostat=ierr)
118  if( ierr < 0 ) then !--- missing
119  log_info("ATMOS_REFSTATE_setup",*) 'Not found namelist. Default used.'
120  elseif( ierr > 0 ) then !--- fatal error
121  log_error("ATMOS_REFSTATE_setup",*) 'Not appropriate names in namelist PARAM_ATMOS_REFSTATE. Check!'
122  call prc_abort
123  endif
124  log_nml(param_atmos_refstate)
125 
126  allocate( atmos_refstate_pres(ka,ia,ja) )
127  allocate( atmos_refstate_temp(ka,ia,ja) )
128  allocate( atmos_refstate_dens(ka,ia,ja) )
129  allocate( atmos_refstate_pott(ka,ia,ja) )
130  allocate( atmos_refstate_qv(ka,ia,ja) )
131  atmos_refstate_pres(:,:,:) = undef
132  atmos_refstate_temp(:,:,:) = undef
133  atmos_refstate_dens(:,:,:) = undef
134  atmos_refstate_pott(:,:,:) = undef
135  atmos_refstate_qv(:,:,:) = undef
136 
137  allocate( atmos_refstate1d_pres(ka) )
138  allocate( atmos_refstate1d_temp(ka) )
139  allocate( atmos_refstate1d_dens(ka) )
140  allocate( atmos_refstate1d_pott(ka) )
141  allocate( atmos_refstate1d_qv(ka) )
142  atmos_refstate1d_pres(:) = undef
143  atmos_refstate1d_temp(:) = undef
144  atmos_refstate1d_dens(:) = undef
145  atmos_refstate1d_pott(:) = undef
146  atmos_refstate1d_qv(:) = undef
147 
148  log_newline
149  log_info("ATMOS_REFSTATE_setup",*) 'Reference state settings'
150  if ( atmos_refstate_in_basename /= '' ) then
151  log_info_cont(*) 'Input file of reference state : ', trim(atmos_refstate_in_basename)
152  else
153  log_info_cont(*) 'Input file of reference state : Nothing, generate internally'
154  endif
155 
156  if ( atmos_refstate_out_basename /= '' ) then
157  log_info_cont(*) 'Output file of reference state : ', trim(atmos_refstate_out_basename)
158  else
159  log_info_cont(*) 'Output file of reference state : No output'
160  endif
161 
162  ! input or generate reference profile
163  if ( atmos_refstate_in_basename /= '' ) then
164 
165  call atmos_refstate_read( ka, ks, ke, ia, is, ie, ja, js, je, & ! [IN]
166  cz(:), fz(:), & ! [IN]
167  real_cz(:,:,:), real_fz(:,:,:), & ! [IN]
168  real_phi(:,:,:) ) ! [IN]
169 
170  else
171  if ( atmos_refstate_type == 'ISA' ) then
172 
173  log_info_cont(*) 'Reference type : ISA'
174  log_info_cont(*) 'Surface temperature [K] : ', atmos_refstate_temp_sfc
175  log_info_cont(*) 'Surface & environment RH [%] : ', atmos_refstate_rh
176 
177  call atmos_refstate_generate_isa( ka, ks, ke, ia, is, ie, ja, js, je, & ! [IN]
178  cz(:), fz(:), & ! [IN]
179  real_cz(:,:,:), real_fz(:,:,:), & ! [IN]
180  real_phi(:,:,:) ) ! [IN]
181 
182  elseif ( atmos_refstate_type == 'UNIFORM' ) then
183 
184  log_info_cont(*) 'Reference type : UNIFORM POTT'
185  log_info_cont(*) 'Potential temperature : ', atmos_refstate_pott_uniform
186 
187  call atmos_refstate_generate_uniform( ka, ks, ke, ia, is, ie, ja, js, je, & ! [IN]
188  cz(:), fz(:), & ! [IN]
189  real_cz(:,:,:), real_fz(:,:,:), & ! [IN]
190  real_phi(:,:,:) ) ! [IN]
191 
192  elseif ( atmos_refstate_type == 'ZERO' ) then
193 
194  log_info_cont(*) 'Reference type : ZERO'
195 
196  call atmos_refstate_generate_zero( ka, ia, ja )
197 
198  elseif ( atmos_refstate_type == 'INIT' ) then
199 
200  if ( atmos_refstate_update_dt >= 0.0_rp ) then
201  atmos_refstate_update_flag = .true.
202  endif
203 
204  log_info_cont(*) 'Reference type : Generate from initial data'
205  log_info_cont(*) 'Update state? : ', atmos_refstate_update_flag
206  log_info_cont(*) 'Update interval [sec] : ', atmos_refstate_update_dt
207 
208  else
209  log_error("ATMOS_REFSTATE_setup",*) 'ATMOS_REFSTATE_TYPE must be "ISA", "UNIFORM", "ZERO", or "INIT". Check! : ', trim(atmos_refstate_type)
210  call prc_abort
211  endif
212 
213  endif
214 
215  return

References atmos_refstate_dens, atmos_refstate_pott, atmos_refstate_pres, atmos_refstate_qv, atmos_refstate_read(), atmos_refstate_temp, atmos_refstate_update_flag, scale_const::const_undef, scale_io::io_fid_conf, and scale_prc::prc_abort().

Referenced by mod_atmos_driver::atmos_driver_setup().

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

◆ atmos_refstate_read()

subroutine, public scale_atmos_refstate::atmos_refstate_read ( integer, intent(in)  KA,
integer, intent(in)  KS,
integer, intent(in)  KE,
integer, intent(in)  IA,
integer, intent(in)  IS,
integer, intent(in)  IE,
integer, intent(in)  JA,
integer, intent(in)  JS,
integer, intent(in)  JE,
real(rp), dimension ( ka), intent(in)  CZ,
real(rp), dimension (0:ka), intent(in)  FZ,
real(rp), dimension ( ka,ia,ja), intent(in)  REAL_CZ,
real(rp), dimension (0:ka,ia,ja), intent(in)  REAL_FZ,
real(rp), dimension( ka,ia,ja), intent(in)  REAL_PHI 
)

Read reference state profile.

Definition at line 223 of file scale_atmos_refstate.F90.

223  use scale_file_cartesc, only: &
224  file_cartesc_open, &
225  file_cartesc_check_coordinates, &
226  file_cartesc_read, &
227  file_cartesc_close
228  use scale_prc, only: &
229  prc_abort
230  implicit none
231 
232  integer, intent(in) :: KA, KS, KE
233  integer, intent(in) :: IA, IS, IE
234  integer, intent(in) :: JA, JS, JE
235  real(RP), intent(in) :: CZ ( KA)
236  real(RP), intent(in) :: FZ (0:KA)
237  real(RP), intent(in) :: REAL_CZ ( KA,IA,JA)
238  real(RP), intent(in) :: REAL_FZ (0:KA,IA,JA)
239  real(RP), intent(in) :: REAL_PHI( KA,IA,JA)
240 
241  integer :: fid
242  !---------------------------------------------------------------------------
243 
244  log_newline
245  log_info("ATMOS_REFSTATE_read",*) 'Input reference state profile '
246 
247  if ( atmos_refstate_in_basename /= '' ) then
248 
249  call file_cartesc_open( atmos_refstate_in_basename, fid )
250 
251  if ( atmos_refstate_in_check_coordinates ) then
252  call file_cartesc_check_coordinates( fid, atmos=.true. )
253  end if
254 
255  ! 1D
256  call file_cartesc_read( fid, 'PRES_ref', 'Z', atmos_refstate1d_pres(:) )
257  call file_cartesc_read( fid, 'TEMP_ref', 'Z', atmos_refstate1d_temp(:) )
258  call file_cartesc_read( fid, 'DENS_ref', 'Z', atmos_refstate1d_dens(:) )
259  call file_cartesc_read( fid, 'POTT_ref', 'Z', atmos_refstate1d_pott(:) )
260  call file_cartesc_read( fid, 'QV_ref', 'Z', atmos_refstate1d_qv(:) )
261 
262  ! 3D
263  call file_cartesc_read( fid, 'PRES_ref3D', 'ZXY', atmos_refstate_pres(:,:,:) )
264  call file_cartesc_read( fid, 'TEMP_ref3D', 'ZXY', atmos_refstate_temp(:,:,:) )
265  call file_cartesc_read( fid, 'DENS_ref3D', 'ZXY', atmos_refstate_dens(:,:,:) )
266  call file_cartesc_read( fid, 'POTT_ref3D', 'ZXY', atmos_refstate_pott(:,:,:) )
267  call file_cartesc_read( fid, 'QV_ref3D', 'ZXY', atmos_refstate_qv(:,:,:) )
268 
269  else
270  log_error("ATMOS_REFSTATE_read",*) 'refstate file is not specified.'
271  call prc_abort
272  endif
273 
274  call atmos_refstate_calc3d( ka, ks, ke, ia, is, ie, ja, js, je, &
275  cz(:), fz(:), real_cz(:,:,:), real_fz(:,:,:), real_phi(:,:,:) )
276 
277  return

References atmos_refstate_calc3d(), atmos_refstate_dens, atmos_refstate_pott, atmos_refstate_pres, atmos_refstate_qv, atmos_refstate_temp, scale_file_cartesc::file_cartesc_close(), scale_file_cartesc::file_cartesc_open(), and scale_prc::prc_abort().

Referenced by atmos_refstate_setup().

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

◆ atmos_refstate_write()

subroutine, public scale_atmos_refstate::atmos_refstate_write

Write reference state profile.

Definition at line 283 of file scale_atmos_refstate.F90.

283  use scale_file_cartesc, only: &
284  file_cartesc_write
285  implicit none
286 
287  logical, save :: first = .true.
288  !---------------------------------------------------------------------------
289 
290  if ( .not. first ) return
291  first = .false.
292 
293  if ( atmos_refstate_out_basename /= '' ) then
294 
295  log_newline
296  log_info("ATMOS_REFSTATE_write",*) 'Output reference state profile '
297 
298  ! 1D
299  call file_cartesc_write( atmos_refstate1d_pres(:), atmos_refstate_out_basename, atmos_refstate_out_title, & ! [IN]
300  'PRES_ref', 'Reference profile of pres.', 'Pa', 'Z', atmos_refstate_out_dtype ) ! [IN]
301  call file_cartesc_write( atmos_refstate1d_temp(:), atmos_refstate_out_basename, atmos_refstate_out_title, & ! [IN]
302  'TEMP_ref', 'Reference profile of temp.', 'K', 'Z', atmos_refstate_out_dtype ) ! [IN]
303  call file_cartesc_write( atmos_refstate1d_dens(:), atmos_refstate_out_basename, atmos_refstate_out_title, & ! [IN]
304  'DENS_ref', 'Reference profile of rho', 'kg/m3', 'Z', atmos_refstate_out_dtype ) ! [IN]
305  call file_cartesc_write( atmos_refstate1d_pott(:), atmos_refstate_out_basename, atmos_refstate_out_title, & ! [IN]
306  'POTT_ref', 'Reference profile of theta', 'K', 'Z', atmos_refstate_out_dtype ) ! [IN]
307  call file_cartesc_write( atmos_refstate1d_qv(:), atmos_refstate_out_basename, atmos_refstate_out_title, & ! [IN]
308  'QV_ref', 'Reference profile of qv', 'kg/kg', 'Z', atmos_refstate_out_dtype ) ! [IN]
309 
310  ! 3D
311  call file_cartesc_write( atmos_refstate_pres(:,:,:), atmos_refstate_out_basename, atmos_refstate_out_title, & ! [IN]
312  'PRES_ref3D', 'Reference profile of pres.', 'Pa', 'ZXY', atmos_refstate_out_dtype ) ! [IN]
313  call file_cartesc_write( atmos_refstate_temp(:,:,:), atmos_refstate_out_basename, atmos_refstate_out_title, & ! [IN]
314  'TEMP_ref3D', 'Reference profile of temp.', 'K', 'ZXY', atmos_refstate_out_dtype ) ! [IN]
315  call file_cartesc_write( atmos_refstate_dens(:,:,:), atmos_refstate_out_basename, atmos_refstate_out_title, & ! [IN]
316  'DENS_ref3D', 'Reference profile of rho', 'kg/m3', 'ZXY', atmos_refstate_out_dtype ) ! [IN]
317  call file_cartesc_write( atmos_refstate_pott(:,:,:), atmos_refstate_out_basename, atmos_refstate_out_title, & ! [IN]
318  'POTT_ref3D', 'Reference profile of theta', 'K', 'ZXY', atmos_refstate_out_dtype ) ! [IN]
319  call file_cartesc_write( atmos_refstate_qv(:,:,:), atmos_refstate_out_basename, atmos_refstate_out_title, & ! [IN]
320  'QV_ref3D', 'Reference profile of qv', 'kg/kg', 'ZXY', atmos_refstate_out_dtype ) ! [IN]
321 
322  endif
323 
324  return

References atmos_refstate_calc3d(), atmos_refstate_dens, atmos_refstate_pott, atmos_refstate_pres, atmos_refstate_qv, atmos_refstate_temp, scale_const::const_epsvap, and scale_const::const_pstd.

Referenced by atmos_refstate_update().

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

◆ atmos_refstate_update()

subroutine, public scale_atmos_refstate::atmos_refstate_update ( integer, intent(in)  KA,
integer, intent(in)  KS,
integer, intent(in)  KE,
integer, intent(in)  IA,
integer, intent(in)  IS,
integer, intent(in)  IE,
integer, intent(in)  ISB,
integer, intent(in)  IEB,
integer, intent(in)  JA,
integer, intent(in)  JS,
integer, intent(in)  JE,
integer, intent(in)  JSB,
integer, intent(in)  JEB,
real(rp), dimension (ka,ia,ja), intent(in)  DENS,
real(rp), dimension (ka,ia,ja), intent(in)  POTT,
real(rp), dimension (ka,ia,ja), intent(in)  TEMP,
real(rp), dimension (ka,ia,ja), intent(in)  PRES,
real(rp), dimension (ka,ia,ja), intent(in)  QV,
real(rp), dimension ( ka), intent(in)  CZ,
real(rp), dimension (0:ka), intent(in)  FZ,
real(rp), dimension ( ka-1), intent(in)  FDZ,
real(rp), dimension ( ka), intent(in)  RCDZ,
real(rp), dimension ( ka,ia,ja), intent(in)  REAL_CZ,
real(rp), dimension (0:ka,ia,ja), intent(in)  REAL_FZ,
real(rp), dimension( ka,ia,ja), intent(in)  REAL_PHI,
real(rp), dimension ( ia,ja), intent(in)  AREA,
real(dp), intent(in)  nowsec,
logical, intent(in), optional  force 
)

Update reference state profile (Horizontal average)

Definition at line 559 of file scale_atmos_refstate.F90.

559  use scale_statistics, only: &
560  statistics_horizontal_mean
561  use scale_interp_vert, only: &
562  interp_vert_xi2z
563  implicit none
564 
565  integer, intent(in) :: KA, KS, KE
566  integer, intent(in) :: IA, IS, IE, ISB, IEB
567  integer, intent(in) :: JA, JS, JE, JSB, JEB
568  real(RP), intent(in) :: DENS (KA,IA,JA)
569  real(RP), intent(in) :: POTT (KA,IA,JA)
570  real(RP), intent(in) :: TEMP (KA,IA,JA)
571  real(RP), intent(in) :: PRES (KA,IA,JA)
572  real(RP), intent(in) :: QV (KA,IA,JA)
573  real(RP), intent(in) :: CZ ( KA)
574  real(RP), intent(in) :: FZ (0:KA)
575  real(RP), intent(in) :: FDZ ( KA-1)
576  real(RP), intent(in) :: RCDZ ( KA)
577  real(RP), intent(in) :: REAL_CZ ( KA,IA,JA)
578  real(RP), intent(in) :: REAL_FZ (0:KA,IA,JA)
579  real(RP), intent(in) :: REAL_PHI( KA,IA,JA)
580  real(RP), intent(in) :: AREA ( IA,JA)
581  real(DP), intent(in) :: nowsec
582  logical, intent(in), optional :: force
583 
584  real(RP) :: work(KA,IA,JA)
585  logical :: force_
586 
587  integer :: k
588  !---------------------------------------------------------------------------
589 
590  if ( present(force) ) then
591  force_ = force
592  else
593  force_ = .false.
594  end if
595 
596  if ( force_ .or. ( nowsec - last_updated >= atmos_refstate_update_dt ) ) then
597 
598  log_info("ATMOS_REFSTATE_update",*) 'update reference state'
599 
600  call interp_vert_xi2z( ka, ks, ke, ia, is, ie, ja, js, je, &
601  cz(:), real_cz(:,:,:), temp(:,:,:), & ! [IN]
602  work(:,:,:) ) ! [OUT]
603  call statistics_horizontal_mean( ka, ks, ke, ia, is, ie, ja, js, je, &
604  work(:,:,:), area(:,:), & ! [IN]
605  atmos_refstate1d_temp(:) ) ! [IN]
606 
607  call interp_vert_xi2z( ka, ks, ke, ia, is, ie, ja, js, je, &
608  cz(:), real_cz(:,:,:), pres(:,:,:), & ! [IN]
609  work(:,:,:) ) ! [OUT]
610  call statistics_horizontal_mean( ka, ks, ke, ia, is, ie, ja, js, je, &
611  work(:,:,:), area(:,:), & ! [IN]
612  atmos_refstate1d_pres(:) ) ! [OUT]
613 
614  call interp_vert_xi2z( ka, ks, ke, ia, is, ie, ja, js, je, &
615  cz(:), real_cz(:,:,:), dens(:,:,:), & ! [IN]
616  work(:,:,:) ) ! [OUT]
617  call statistics_horizontal_mean( ka, ks, ke, ia, is, ie, ja, js, je, &
618  work(:,:,:), area(:,:), & ! [IN]
619  atmos_refstate1d_dens(:) ) ! [OUT]
620 
621  call interp_vert_xi2z( ka, ks, ke, ia, is, ie, ja, js, je, &
622  cz(:), real_cz(:,:,:), pott(:,:,:), & ! [IN]
623  work(:,:,:) ) ! [OUT]
624  call statistics_horizontal_mean( ka, ks, ke, ia, is, ie, ja, js, je, &
625  work(:,:,:), area(:,:), & ! [IN]
626  atmos_refstate1d_pott(:) ) ! [OUT]
627 
628  call interp_vert_xi2z( ka, ks, ke, ia, is, ie, ja, js, je, &
629  cz(:), real_cz(:,:,:), qv(:,:,:), & ! [IN]
630  work(:,:,:) ) ! [OUT]
631  call statistics_horizontal_mean( ka, ks, ke, ia, is, ie, ja, js, je, &
632  work(:,:,:), area(:,:), & ! [IN]
633  atmos_refstate1d_qv(:) ) ! [OUT]
634 
635  do k = ke-1, ks, -1 ! fill undefined value
636  if( atmos_refstate1d_dens(k) <= 0.0_rp ) atmos_refstate1d_dens(k) = atmos_refstate1d_dens(k+1)
637  if( atmos_refstate1d_temp(k) <= 0.0_rp ) atmos_refstate1d_temp(k) = atmos_refstate1d_temp(k+1)
638  if( atmos_refstate1d_pres(k) <= 0.0_rp ) atmos_refstate1d_pres(k) = atmos_refstate1d_pres(k+1)
639  if( atmos_refstate1d_pott(k) <= 0.0_rp ) atmos_refstate1d_pott(k) = atmos_refstate1d_pott(k+1)
640  if( atmos_refstate1d_qv(k) <= 0.0_rp ) atmos_refstate1d_qv(k) = atmos_refstate1d_qv(k+1)
641  enddo
642 
643  call atmos_refstate_smoothing( ka, ks, ke, fdz(:), rcdz(:), atmos_refstate1d_pott(:) )
644  call atmos_refstate_smoothing( ka, ks, ke, fdz(:), rcdz(:), atmos_refstate1d_qv(:) )
645 
646  call atmos_refstate_calc3d( ka, ks, ke, ia, isb, ieb, ja, jsb, jeb, &
647  cz(:), fz(:), real_cz(:,:,:), real_fz(:,:,:), real_phi(:,:,:) )
648 
649  last_updated = nowsec
650 
651 
652  log_newline
653  log_info("ATMOS_REFSTATE_update",*) 'Generated reference state of atmosphere:'
654  log_info_cont(*) '============================================================================='
655  log_info_cont(*) ' z*-coord.: pressure: temperature: density: pot.temp.: water vapor'
656  do k = ks, ke
657  log_info_cont('(6F13.5)') cz(k), &
658  atmos_refstate1d_pres(k), &
659  atmos_refstate1d_temp(k), &
660  atmos_refstate1d_dens(k), &
661  atmos_refstate1d_pott(k), &
662  atmos_refstate1d_qv(k)
663  enddo
664  log_info_cont(*) '============================================================================='
665 
666  ! output reference profile
667  call atmos_refstate_write
668 
669  endif
670 
671  return

References atmos_refstate_calc3d(), atmos_refstate_smoothing(), atmos_refstate_write(), and scale_interp_vert::interp_vert_xi2z().

Referenced by mod_atmos_driver::atmos_driver_update(), and mod_rm_driver::restart_read().

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

◆ atmos_refstate_calc3d()

subroutine scale_atmos_refstate::atmos_refstate_calc3d ( integer, intent(in)  KA,
integer, intent(in)  KS,
integer, intent(in)  KE,
integer, intent(in)  IA,
integer, intent(in)  IS,
integer, intent(in)  IE,
integer, intent(in)  JA,
integer, intent(in)  JS,
integer, intent(in)  JE,
real(rp), dimension ( ka), intent(in)  CZ,
real(rp), dimension (0:ka), intent(in)  FZ,
real(rp), dimension ( ka,ia,ja), intent(in)  REAL_CZ,
real(rp), dimension (0:ka,ia,ja), intent(in)  REAL_FZ,
real(rp), dimension( ka,ia,ja), intent(in)  REAL_PHI 
)

apply 1D reference to 3D (terrain-following) with re-calc hydrostatic balance

Definition at line 679 of file scale_atmos_refstate.F90.

679  use scale_const, only: &
680  undef => const_undef, &
681  rdry => const_rdry, &
682  cpdry => const_cpdry, &
683  p00 => const_pre00
684  use scale_comm_cartesc, only: &
685  comm_vars8, &
686  comm_wait
687  use scale_interp_vert, only: &
688  interp_vert_z2xi
689  use scale_atmos_hydrostatic, only: &
690  hydrostatic_buildrho_atmos => atmos_hydrostatic_buildrho_atmos, &
691  hydrostatic_buildrho_atmos_rev => atmos_hydrostatic_buildrho_atmos_rev
692  implicit none
693  integer, intent(in) :: KA, KS, KE
694  integer, intent(in) :: IA, IS, IE
695  integer, intent(in) :: JA, JS, JE
696 
697  real(RP), intent(in) :: CZ ( KA)
698  real(RP), intent(in) :: FZ (0:KA)
699  real(RP), intent(in) :: REAL_CZ ( KA,IA,JA)
700  real(RP), intent(in) :: REAL_FZ (0:KA,IA,JA)
701  real(RP), intent(in) :: REAL_PHI( KA,IA,JA)
702 
703  real(RP) :: dens(KA,IA,JA)
704  real(RP) :: temp(KA,IA,JA)
705  real(RP) :: pres(KA,IA,JA)
706  real(RP) :: pott(KA,IA,JA)
707  real(RP) :: qv (KA,IA,JA)
708  real(RP) :: qc (KA,IA,JA)
709  real(RP) :: dz (KA,IA,JA)
710 
711  real(RP) :: dens_toa_1D
712  real(RP) :: temp_toa_1D
713  real(RP) :: pres_toa_1D
714  real(RP) :: qc_1D
715  real(RP) :: dz_1D
716 
717  real(RP) :: work(KA,IA,JA)
718  real(RP) :: RovCP
719  integer :: k, i, j
720  !---------------------------------------------------------------------------
721 
722  rovcp = rdry / cpdry
723 
724  !--- potential temperature
725  do j = js, je
726  do i = is, ie
727  work(:,i,j) = atmos_refstate1d_pott(:)
728  enddo
729  enddo
730 
731  call interp_vert_z2xi( ka, ks, ke, ia, is, ie, ja, js, je, & ! [IN]
732  real_cz(:,:,:), cz(:), & ! [IN]
733  work(:,:,:), & ! [IN]
734  pott(:,:,:) ) ! [OUT]
735 
736  !--- water vapor
737  do j = js, je
738  do i = is, ie
739  work(:,i,j) = atmos_refstate1d_qv(:)
740  enddo
741  enddo
742 
743  call interp_vert_z2xi( ka, ks, ke, ia, is, ie, ja, js, je, & ! [IN]
744  real_cz(:,:,:), cz(:), & ! [IN]
745  work(:,:,:), & ! [IN]
746  qv(:,:,:) ) ! [OUT]
747 
748 
749  !--- build up density to TOA (1D)
750  qc_1d = 0.0_rp
751  dz_1d = fz(ke) - cz(ke)
752 
753  call hydrostatic_buildrho_atmos( atmos_refstate1d_pott(ke), & ! [IN]
754  atmos_refstate1d_qv(ke), & ! [IN]
755  qc_1d, & ! [IN]
756  atmos_refstate1d_dens(ke), & ! [IN]
757  atmos_refstate1d_pott(ke), & ! [IN]
758  atmos_refstate1d_qv(ke), & ! [IN]
759  qc_1d, & ! [IN]
760  dz_1d, & ! [IN]
761  ke+1, & ! [IN]
762  dens_toa_1d, & ! [OUT]
763  temp_toa_1d, & ! [OUT]
764  pres_toa_1d ) ! [OUT]
765 
766  ! build down density from TOA (3D)
767  do j = js, je
768  do i = is, ie
769  dz(ks-1,i,j) = real_cz(ks,i,j) - real_fz(ks-1,i,j) ! distance from surface to cell center
770  do k = ks, ke-1
771  dz(k,i,j) = real_cz(k+1,i,j) - real_cz(k,i,j) ! distance from cell center to cell center
772  enddo
773  dz(ke,i,j) = real_fz(ke,i,j) - real_cz(ke,i,j) ! distance from cell center to TOA
774  enddo
775  enddo
776 
777  do j = js, je
778  do i = is, ie
779  dens(ke+1,i,j) = dens_toa_1d
780  temp(ke+1,i,j) = temp_toa_1d
781  pres(ke+1,i,j) = pres_toa_1d
782  pott(ke+1,i,j) = pott(ke,i,j)
783  qv(ke+1,i,j) = qv(ke,i,j)
784  enddo
785  enddo
786 
787  do j = js, je
788  do i = is, ie
789  pott(ks-1,i,j) = pott(ks,i,j)
790  qv(ks-1,i,j) = qv(ks,i,j)
791  enddo
792  enddo
793 
794  qc(:,:,:) = 0.0_rp
795 
796  call hydrostatic_buildrho_atmos_rev( ia, is, ie, ja, js, je, &
797  pott(ke,:,:), qv(ke,:,:), qc(ke,:,:), & ! [IN]
798  dens(ke+1,:,:), pott(ke+1,:,:), qv(ke+1,:,:), qc(ke+1,:,:), & ! [IN]
799  dz(ke,:,:), ke+1, & ! [IN]
800  dens(ke ,:,:), temp(ke ,:,:), pres(ke ,:,:) ) ! [OUT]
801 
802  call hydrostatic_buildrho_atmos_rev( ka, ks, ke, ia, is, ie, ja, js, je, &
803  pott(:,:,:), qv(:,:,:), qc(:,:,:), & ! [IN]
804  dz(:,:,:), & ! [IN]
805  dens(:,:,:), & ! [INOUT]
806  temp(:,:,:), pres(:,:,:) ) ! [OUT]
807 
808  do j = js, je
809  do i = is, ie
810  do k = ks, ke
811  atmos_refstate_dens(k,i,j) = dens(k,i,j)
812  atmos_refstate_temp(k,i,j) = temp(k,i,j)
813  atmos_refstate_pres(k,i,j) = pres(k,i,j)
814  atmos_refstate_pott(k,i,j) = pott(k,i,j)
815  atmos_refstate_qv(k,i,j) = qv(k,i,j)
816  enddo
817  enddo
818  enddo
819 
820  ! boundary condition
821  do j = js, je
822  do i = is, ie
823 
824  atmos_refstate_temp(1:ks-1, i,j) = temp(ks,i,j)
825  atmos_refstate_temp(ke+1:ka,i,j) = temp_toa_1d
826 
827  atmos_refstate_qv(1:ks-1, i,j) = qv(ks,i,j)
828  atmos_refstate_qv(ke+1:ka,i,j) = qv(ke,i,j)
829 
830  atmos_refstate_pres(1:ks-2, i,j) = undef
831  atmos_refstate_pres(ks-1, i,j) = atmos_refstate_pres(ks+1,i,j) &
832  - atmos_refstate_dens(ks ,i,j) * ( real_phi(ks-1,i,j) - real_phi(ks+1,i,j) )
833  atmos_refstate_pres(ke+1, i,j) = atmos_refstate_pres(ke-1,i,j) &
834  - atmos_refstate_dens(ke ,i,j) * ( real_phi(ke+1,i,j) - real_phi(ke-1,i,j) )
835  atmos_refstate_pres(ke+2:ka,i,j) = undef
836 
837  atmos_refstate_dens(1:ks-2, i,j) = undef
838  atmos_refstate_dens(ks-1, i,j) = atmos_refstate_pres(ks-1,i,j) / ( atmos_refstate_temp(ks-1,i,j) * rdry )
839  atmos_refstate_dens(ke+1, i,j) = atmos_refstate_pres(ke+1,i,j) / ( atmos_refstate_temp(ke+1,i,j) * rdry )
840  atmos_refstate_dens(ke+2:ka,i,j) = undef
841 
842  atmos_refstate_pott(1:ks-2, i,j) = undef
843  atmos_refstate_pott(ks-1, i,j) = atmos_refstate_temp(ks-1,i,j) * ( p00 / atmos_refstate_pres(ks-1,i,j) )**rovcp
844  atmos_refstate_pott(ke+1, i,j) = atmos_refstate_temp(ke+1,i,j) * ( p00 / atmos_refstate_pres(ke+1,i,j) )**rovcp
845  atmos_refstate_pott(ke+2:ka,i,j) = undef
846  enddo
847  enddo
848 
849  call comm_vars8( atmos_refstate_dens(:,:,:), 1 )
850  call comm_vars8( atmos_refstate_temp(:,:,:), 2 )
851  call comm_vars8( atmos_refstate_pres(:,:,:), 3 )
852  call comm_vars8( atmos_refstate_pott(:,:,:), 4 )
853  call comm_vars8( atmos_refstate_qv(:,:,:), 5 )
854  call comm_wait ( atmos_refstate_dens(:,:,:), 1, .false. )
855  call comm_wait ( atmos_refstate_temp(:,:,:), 2, .false. )
856  call comm_wait ( atmos_refstate_pres(:,:,:), 3, .false. )
857  call comm_wait ( atmos_refstate_pott(:,:,:), 4, .false. )
858  call comm_wait ( atmos_refstate_qv(:,:,:), 5, .false. )
859 
860  return

References atmos_refstate_dens, atmos_refstate_pott, atmos_refstate_pres, atmos_refstate_qv, atmos_refstate_temp, scale_const::const_cpdry, scale_const::const_pre00, scale_const::const_rdry, scale_const::const_undef, and scale_interp_vert::interp_vert_z2xi().

Referenced by atmos_refstate_read(), atmos_refstate_update(), and atmos_refstate_write().

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

◆ atmos_refstate_smoothing()

subroutine scale_atmos_refstate::atmos_refstate_smoothing ( integer, intent(in)  KA,
integer, intent(in)  KS,
integer, intent(in)  KE,
real(rp), dimension (ka-1), intent(in)  FDZ,
real(rp), dimension(ka), intent(in)  RCDZ,
real(rp), dimension(ka), intent(inout)  phi 
)

Definition at line 869 of file scale_atmos_refstate.F90.

869  use scale_const, only: &
870  eps => const_eps
871  implicit none
872  integer, intent(in) :: KA, KS, KE
873 
874  real(RP), intent(in) :: FDZ (KA-1)
875  real(RP), intent(in) :: RCDZ(KA)
876 
877  real(RP), intent(inout) :: phi(KA)
878 
879  real(RP) :: dev (KA)
880  real(RP) :: correction(KA)
881  real(RP) :: fact(KA)
882 
883  integer, parameter :: iter_max = 100
884  real(RP) :: sig0, sig1, zerosw
885  logical :: updated
886  integer :: k, iter
887  !---------------------------------------------------------------------------
888 
889  dev(ks) = 0.0_rp
890  dev(ke) = 0.0_rp
891 
892  correction(ks-1:ks+1) = 0.0_rp
893  correction(ke-1:ke+1) = 0.0_rp
894 
895  fact(ks-1:ks+1) = 0.0_rp
896  fact(ke-1:ke+1) = 0.0_rp
897 
898  do iter = 1, iter_max
899  updated = .false.
900 
901  do k = ks+1, ke-1
902  dev(k) = phi(k) - ( fdz(k-1)*phi(k+1) + fdz(k)*phi(k-1) ) / ( fdz(k) + fdz(k-1) )
903  enddo
904 
905  do k = ks+2, ke-2
906  sig0 = dev(k) * dev(k-1)
907  sig1 = dev(k) * dev(k+1)
908  if ( sig0 < -eps .and. sig1 < -eps ) then
909  correction(k) = dev(k) &
910  / ( 2.0_rp*rcdz(k) + ( fdz(k-1)*rcdz(k+1) + fdz(k)*rcdz(k-1) ) / ( fdz(k) + fdz(k-1) ) )
911  updated = .true.
912  else
913  correction(k) = 0.0_rp
914  end if
915  enddo
916 
917  sig1 = dev(ks+1) * dev(ks+2)
918  if ( sig1 < -eps ) then
919  correction(ks+1) = dev(ks+1) &
920  / ( 2.0_rp*rcdz(ks+1) + (fdz(ks)*rcdz(ks+2)+fdz(ks+1)*rcdz(ks))/(fdz(ks+1)+fdz(ks)) )
921  updated = .true.
922  else
923  correction(ks+1) = 0.0_rp
924  end if
925 
926  sig0 = dev(ke-1) * dev(ke-2)
927  if ( sig0 < -eps ) then
928  correction(ke-1) = dev(ke-1) &
929  / ( 2.0_rp*rcdz(ke-1) + (fdz(ke-2)*rcdz(ke)+fdz(ke-1)*rcdz(ke-2))/(fdz(ke-1)+fdz(ke-2)) )
930  updated = .true.
931  else
932  correction(ke-1) = 0.0_rp
933  end if
934 
935  if ( .NOT. updated ) exit
936 
937  do k = ks+1, ke-1
938  zerosw = 0.5_rp - sign( 0.5_rp, abs(correction(k))-eps ) ! if correction(k) == 0 then fact(k) = 0.0
939  fact(k) = correction(k) / ( correction(k) - correction(k+1) - correction(k-1) + zerosw )
940  enddo
941 
942  do k = ks, ke
943  phi(k) = phi(k) + ( correction(k+1) * fact(k+1) &
944  - correction(k ) * fact(k ) * 2.0_rp &
945  + correction(k-1) * fact(k-1) ) * rcdz(k)
946  enddo
947 
948  if ( iter == iter_max ) then
949  log_info("ATMOS_REFSTATE_smoothing",*) "iteration not converged!", phi
950  endif
951  enddo
952 
953  return

References scale_const::const_eps.

Referenced by atmos_refstate_update().

Here is the caller graph for this function:

Variable Documentation

◆ atmos_refstate_update_flag

logical, public scale_atmos_refstate::atmos_refstate_update_flag = .false.

Definition at line 36 of file scale_atmos_refstate.F90.

36  logical, public :: ATMOS_REFSTATE_UPDATE_FLAG = .false.

Referenced by mod_atmos_driver::atmos_driver_update(), and atmos_refstate_setup().

◆ atmos_refstate_pres

real(rp), dimension(:,:,:), allocatable, public scale_atmos_refstate::atmos_refstate_pres

refernce pressure [Pa]

Definition at line 38 of file scale_atmos_refstate.F90.

38  real(RP), public, allocatable :: ATMOS_REFSTATE_pres(:,:,:)

Referenced by mod_atmos_dyn_driver::atmos_dyn_driver(), atmos_refstate_calc3d(), atmos_refstate_read(), atmos_refstate_setup(), and atmos_refstate_write().

◆ atmos_refstate_temp

real(rp), dimension(:,:,:), allocatable, public scale_atmos_refstate::atmos_refstate_temp

refernce temperature [K]

Definition at line 39 of file scale_atmos_refstate.F90.

39  real(RP), public, allocatable :: ATMOS_REFSTATE_temp(:,:,:)

Referenced by atmos_refstate_calc3d(), atmos_refstate_read(), atmos_refstate_setup(), and atmos_refstate_write().

◆ atmos_refstate_dens

real(rp), dimension(:,:,:), allocatable, public scale_atmos_refstate::atmos_refstate_dens

refernce density [kg/m3]

Definition at line 40 of file scale_atmos_refstate.F90.

40  real(RP), public, allocatable :: ATMOS_REFSTATE_dens(:,:,:)

Referenced by mod_atmos_bnd_driver::atmos_boundary_driver_set(), mod_atmos_dyn_driver::atmos_dyn_driver(), mod_atmos_phy_mp_driver::atmos_phy_mp_driver_calc_tendency(), atmos_refstate_calc3d(), atmos_refstate_read(), atmos_refstate_setup(), atmos_refstate_write(), and mod_atmos_bnd_driver::calc_mass().

◆ atmos_refstate_pott

real(rp), dimension(:,:,:), allocatable, public scale_atmos_refstate::atmos_refstate_pott

refernce potential temperature [K]

Definition at line 41 of file scale_atmos_refstate.F90.

41  real(RP), public, allocatable :: ATMOS_REFSTATE_pott(:,:,:)

Referenced by mod_atmos_dyn_driver::atmos_dyn_driver(), atmos_refstate_calc3d(), atmos_refstate_read(), atmos_refstate_setup(), and atmos_refstate_write().

◆ atmos_refstate_qv

real(rp), dimension (:,:,:), allocatable, public scale_atmos_refstate::atmos_refstate_qv

refernce vapor [kg/kg]

Definition at line 42 of file scale_atmos_refstate.F90.

42  real(RP), public, allocatable :: ATMOS_REFSTATE_qv (:,:,:)

Referenced by mod_atmos_dyn_driver::atmos_dyn_driver(), atmos_refstate_calc3d(), atmos_refstate_read(), atmos_refstate_setup(), and atmos_refstate_write().

scale_statistics
module Statistics
Definition: scale_statistics.F90:11
scale_prc::prc_abort
subroutine, public prc_abort
Abort Process.
Definition: scale_prc.F90:342
scale_const::const_eps
real(rp), public const_eps
small number
Definition: scale_const.F90:33
scale_interp_vert::interp_vert_z2xi
subroutine, public interp_vert_z2xi(KA, KS, KE, IA, IS, IE, JA, JS, JE, Z, Xi, var, var_Xi)
Definition: scale_interp_vert.F90:227
scale_prc
module PROCESS
Definition: scale_prc.F90:11
scale_const
module CONSTANT
Definition: scale_const.F90:11
scale_file_cartesc::file_cartesc_close
subroutine, public file_cartesc_close(fid)
Close a netCDF file.
Definition: scale_file_cartesC.F90:1023
scale_const::const_cpdry
real(rp), public const_cpdry
specific heat (dry air,constant pressure) [J/kg/K]
Definition: scale_const.F90:56
scale_interp_vert::interp_vert_xi2z
subroutine, public interp_vert_xi2z(KA, KS, KE, IA, IS, IE, JA, JS, JE, Xi, Z, var, var_Z)
Definition: scale_interp_vert.F90:186
scale_atmos_hydrostatic
module atmosphere / hydrostatic barance
Definition: scale_atmos_hydrostatic.F90:12
scale_file_cartesc::file_cartesc_open
subroutine, public file_cartesc_open(basename, fid, aggregate)
open a netCDF file for read
Definition: scale_file_cartesC.F90:746
scale_comm_cartesc
module COMMUNICATION
Definition: scale_comm_cartesC.F90:11
scale_const::const_rdry
real(rp), public const_rdry
specific gas constant (dry air) [J/kg/K]
Definition: scale_const.F90:55
scale_const::const_undef
real(rp), public const_undef
Definition: scale_const.F90:41
scale_interp_vert
module INTERPOLATION vertical
Definition: scale_interp_vert.F90:11
scale_const::const_pre00
real(rp), public const_pre00
pressure reference [Pa]
Definition: scale_const.F90:88
scale_file_cartesc
module file / cartesianC
Definition: scale_file_cartesC.F90:11
Generated by   doxygen 1.8.17