SCALE-RM
scale_atmos_sub_hydrostatic.F90
Go to the documentation of this file.
1 !-------------------------------------------------------------------------------
10 !-------------------------------------------------------------------------------
11 module scale_atmos_hydrostatic
12  !-----------------------------------------------------------------------------
13  !
14  !++ used modules
15  !
16  use scale_precision
17  use scale_stdio
18  use scale_prof
19  use scale_grid_index
20 
21  use scale_const, only: &
22  grav => const_grav, &
23  rdry => const_rdry, &
24  rvap => const_rvap, &
25  cvdry => const_cvdry, &
26  cvvap => const_cvvap, &
27  cl => const_cl, &
28  laps => const_laps, &
29  lapsdry => const_lapsdry, &
30  p00 => const_pre00, &
31  thermodyn_type => const_thermodyn_type
32  use scale_grid, only: &
33  grid_cz, &
34  grid_fdz
35  use scale_grid_real, only: &
36  real_cz, &
37  real_fz
38  !-----------------------------------------------------------------------------
39  implicit none
40  private
41  !-----------------------------------------------------------------------------
42  !
43  !++ Public procedure
44  !
45  public :: atmos_hydrostatic_setup
46  public :: atmos_hydrostatic_buildrho
47  public :: atmos_hydrostatic_buildrho_real
48  public :: atmos_hydrostatic_buildrho_atmos
49  public :: atmos_hydrostatic_buildrho_bytemp
50  public :: atmos_hydrostatic_buildrho_bytemp_atmos
51 
52  public :: atmos_hydrostatic_buildrho_atmos_0d
53  public :: atmos_hydrostatic_buildrho_atmos_rev_2d
54  public :: atmos_hydrostatic_buildrho_atmos_rev_3d
55 
56  public :: atmos_hydrostatic_barometric_law_mslp
57  public :: atmos_hydrostatic_barometric_law_pres
58 
59  interface atmos_hydrostatic_buildrho
60  module procedure atmos_hydrostatic_buildrho_1d
61  module procedure atmos_hydrostatic_buildrho_3d
62  end interface atmos_hydrostatic_buildrho
63 
64  interface atmos_hydrostatic_buildrho_real
65 ! module procedure ATMOS_HYDROSTATIC_buildrho_1D ! buildrho_real_1D is completely equal to buildrho_1D
66  module procedure atmos_hydrostatic_buildrho_real_3d
67  end interface atmos_hydrostatic_buildrho_real
68 
69  interface atmos_hydrostatic_buildrho_atmos
70  module procedure atmos_hydrostatic_buildrho_atmos_1d
71  module procedure atmos_hydrostatic_buildrho_atmos_3d
72  end interface atmos_hydrostatic_buildrho_atmos
73 
74  interface atmos_hydrostatic_buildrho_bytemp
75  module procedure atmos_hydrostatic_buildrho_bytemp_1d
76  module procedure atmos_hydrostatic_buildrho_bytemp_3d
77  end interface atmos_hydrostatic_buildrho_bytemp
78 
79  interface atmos_hydrostatic_buildrho_bytemp_atmos
80  module procedure atmos_hydrostatic_buildrho_bytemp_atmos_1d
81  module procedure atmos_hydrostatic_buildrho_bytemp_atmos_3d
82  end interface atmos_hydrostatic_buildrho_bytemp_atmos
83 
84  interface atmos_hydrostatic_barometric_law_mslp
85  module procedure atmos_hydrostatic_barometric_law_mslp_0d
86  module procedure atmos_hydrostatic_barometric_law_mslp_2d
87  end interface atmos_hydrostatic_barometric_law_mslp
88 
89  interface atmos_hydrostatic_barometric_law_pres
90  module procedure atmos_hydrostatic_barometric_law_pres_0d
91  module procedure atmos_hydrostatic_barometric_law_pres_2d
92  end interface atmos_hydrostatic_barometric_law_pres
93 
94  !-----------------------------------------------------------------------------
95  !
96  !++ Public parameters & variables
97  !
98  !-----------------------------------------------------------------------------
99  !
100  !++ Private procedure
101  !
102  private :: atmos_hydrostatic_buildrho_atmos_2d
103 
104  !-----------------------------------------------------------------------------
105  !
106  !++ Private parameters & variables
107  !
108  integer, private, parameter :: itelim = 100
109  real(RP), private :: criteria
110  logical, private :: hydrostatic_uselapserate = .false.
111  integer, private :: hydrostatic_buildrho_real_kref = 1
112 
113 
114  real(RP), private :: cv_qv
115  real(RP), private :: cv_qc
116 
117  !-----------------------------------------------------------------------------
118 contains
119  !-----------------------------------------------------------------------------
121  subroutine atmos_hydrostatic_setup
122  use scale_process, only: &
123  prc_mpistop
124  use scale_const, only: &
125  const_eps
126  implicit none
127 
128  namelist / param_atmos_hydrostatic / &
129  hydrostatic_uselapserate, &
130  hydrostatic_buildrho_real_kref
131 
132  integer :: ierr
133  !---------------------------------------------------------------------------
134 
135  if( io_l ) write(io_fid_log,*)
136  if( io_l ) write(io_fid_log,*) '++++++ Module[HYDROSTATIC] / Categ[ATMOS SHARE] / Origin[SCALElib]'
137 
138  !--- read namelist
139  rewind(io_fid_conf)
140  read(io_fid_conf,nml=param_atmos_hydrostatic,iostat=ierr)
141  if( ierr < 0 ) then !--- missing
142  if( io_l ) write(io_fid_log,*) '*** Not found namelist. Default used.'
143  elseif( ierr > 0 ) then !--- fatal error
144  write(*,*) 'xxx Not appropriate names in namelist PARAM_ATMOS_HYDROSTATIC. Check!'
145  call prc_mpistop
146  endif
147  if( io_lnml ) write(io_fid_log,nml=param_atmos_hydrostatic)
148 
149  criteria = sqrt( const_eps )
150 
151  if( io_l ) write(io_fid_log,*)
152  if( io_l ) write(io_fid_log,*) '*** use lapse rate for estimation of surface temperature? : ', hydrostatic_uselapserate
153  if( io_l ) write(io_fid_log,*) '*** buildrho conversion criteria : ', criteria
154 
155  if ( thermodyn_type == 'EXACT' ) then
156  cv_qv = cvvap
157  cv_qc = cl
158  elseif( thermodyn_type == 'SIMPLE' ) then
159  cv_qv = cvdry
160  cv_qc = cvdry
161  endif
162 
163  return
164  end subroutine atmos_hydrostatic_setup
165 
166  !-----------------------------------------------------------------------------
168  subroutine atmos_hydrostatic_buildrho_1d( &
169  dens, &
170  temp, &
171  pres, &
172  pott, &
173  qv, &
174  qc, &
175  temp_sfc, &
176  pres_sfc, &
177  pott_sfc, &
178  qv_sfc, &
179  qc_sfc )
180  use scale_process, only: &
181  prc_mpistop
182  implicit none
183 
184  real(RP), intent(out) :: dens(ka)
185  real(RP), intent(out) :: temp(ka)
186  real(RP), intent(out) :: pres(ka)
187  real(RP), intent(in) :: pott(ka)
188  real(RP), intent(in) :: qv (ka)
189  real(RP), intent(in) :: qc (ka)
190 
191  real(RP), intent(out) :: temp_sfc
192  real(RP), intent(in) :: pres_sfc
193  real(RP), intent(in) :: pott_sfc
194  real(RP), intent(in) :: qv_sfc
195  real(RP), intent(in) :: qc_sfc
196 
197  real(RP) :: dens_sfc
198 
199  real(RP) :: Rtot_sfc
200  real(RP) :: CVtot_sfc
201  real(RP) :: CPovCV_sfc
202  real(RP) :: Rtot
203  real(RP) :: CVtot
204  real(RP) :: CPovCV
205 
206  real(RP) :: CVovCP_sfc, CPovR, CVovCP, RovCV
207  real(RP) :: dens_s, dhyd, dgrd
208  integer :: ite
209  logical :: converged
210  !---------------------------------------------------------------------------
211 
212  !--- from surface to lowermost atmosphere
213 
214  rtot_sfc = rdry * ( 1.0_rp - qv_sfc - qc_sfc ) &
215  + rvap * qv_sfc
216  cvtot_sfc = cvdry * ( 1.0_rp - qv_sfc - qc_sfc ) &
217  + cv_qv * qv_sfc &
218  + cv_qc * qc_sfc
219  cpovcv_sfc = ( cvtot_sfc + rtot_sfc ) / cvtot_sfc
220 
221  rtot = rdry * ( 1.0_rp - qv(ks) - qc(ks) ) &
222  + rvap * qv(ks)
223  cvtot = cvdry * ( 1.0_rp - qv(ks) - qc(ks) ) &
224  + cv_qv * qv(ks) &
225  + cv_qc * qc(ks)
226  cpovcv = ( cvtot + rtot ) / cvtot
227 
228  ! density at surface
229  cvovcp_sfc = 1.0_rp / cpovcv_sfc
230  dens_sfc = p00 / rtot_sfc / pott_sfc * ( pres_sfc/p00 )**cvovcp_sfc
231  temp_sfc = pres_sfc / ( dens_sfc * rtot_sfc )
232 
233  ! make density at lowermost cell center
234  if ( hydrostatic_uselapserate ) then
235 
236  cpovr = ( cvtot + rtot ) / rtot
237  cvovcp = 1.0_rp / cpovcv
238 
239  temp(ks) = pott_sfc - lapsdry * grid_cz(ks) ! use dry lapse rate
240  pres(ks) = p00 * ( temp(ks)/pott(ks) )**cpovr
241  dens(ks) = p00 / rtot / pott(ks) * ( pres(ks)/p00 )**cvovcp
242 
243  else ! use itelation
244 
245  rovcv = rtot / cvtot
246 
247  dens_s = 0.0_rp
248  dens(ks) = dens_sfc ! first guess
249 
250  converged = .false.
251  do ite = 1, itelim
252  if ( abs(dens(ks)-dens_s) <= criteria ) then
253  converged = .true.
254  exit
255  endif
256 
257  dens_s = dens(ks)
258 
259  dhyd = + ( p00 * ( dens_sfc * rtot_sfc * pott_sfc / p00 )**cpovcv_sfc &
260  - p00 * ( dens_s * rtot * pott(ks) / p00 )**cpovcv ) / grid_cz(ks) & ! dp/dz
261  - grav * 0.5_rp * ( dens_sfc + dens_s ) ! rho*g
262 
263  dgrd = - p00 * ( rtot * pott(ks) / p00 )**cpovcv / grid_cz(ks) &
264  * cpovcv * dens_s**rovcv &
265  - 0.5_rp * grav
266 
267  dens(ks) = dens_s - dhyd/dgrd
268 
269  if( dens(ks)*0.0_rp /= 0.0_rp) exit
270  enddo
271 
272  if ( .NOT. converged ) then
273  write(*,*) 'xxx [buildrho 1D sfc] iteration not converged!', &
274  dens(ks),ite,dens_s,dhyd,dgrd
275  call prc_mpistop
276  endif
277 
278  endif
279 
280  !--- from lowermost atmosphere to top of atmosphere
281  call atmos_hydrostatic_buildrho_atmos_1d( dens(:), & ! [INOUT]
282  temp(:), & ! [OUT]
283  pres(:), & ! [OUT]
284  pott(:), & ! [IN]
285  qv(:), & ! [IN]
286  qc(:) ) ! [IN]
287 
288  return
289  end subroutine atmos_hydrostatic_buildrho_1d
290 
291  !-----------------------------------------------------------------------------
293  subroutine atmos_hydrostatic_buildrho_3d( &
294  dens, &
295  temp, &
296  pres, &
297  pott, &
298  qv, &
299  qc, &
300  temp_sfc, &
301  pres_sfc, &
302  pott_sfc, &
303  qv_sfc, &
304  qc_sfc )
305  use scale_process, only: &
306  prc_mpistop
307  use scale_comm, only: &
308  comm_horizontal_mean
309  implicit none
310 
311  real(RP), intent(out) :: dens(ka,ia,ja)
312  real(RP), intent(out) :: temp(ka,ia,ja)
313  real(RP), intent(out) :: pres(ka,ia,ja)
314  real(RP), intent(in) :: pott(ka,ia,ja)
315  real(RP), intent(in) :: qv (ka,ia,ja)
316  real(RP), intent(in) :: qc (ka,ia,ja)
317 
318  real(RP), intent(out) :: temp_sfc(1,ia,ja)
319  real(RP), intent(in) :: pres_sfc(1,ia,ja)
320  real(RP), intent(in) :: pott_sfc(1,ia,ja)
321  real(RP), intent(in) :: qv_sfc (1,ia,ja)
322  real(RP), intent(in) :: qc_sfc (1,ia,ja)
323 
324  real(RP) :: dz(ka,ia,ja)
325 
326  real(RP) :: dens_sfc (1,ia,ja)
327  real(RP) :: pott_toa (ia,ja)
328  real(RP) :: qv_toa (ia,ja)
329  real(RP) :: qc_toa (ia,ja)
330  real(RP) :: dens_1D (ka)
331 
332  real(RP) :: Rtot_sfc (ia,ja)
333  real(RP) :: CVtot_sfc (ia,ja)
334  real(RP) :: CPovCV_sfc(ia,ja)
335  real(RP) :: Rtot (ia,ja)
336  real(RP) :: CVtot (ia,ja)
337  real(RP) :: CPovCV (ia,ja)
338 
339  real(RP) :: CVovCP_sfc, CPovR, CVovCP
340 
341  integer :: k, i, j
342  !---------------------------------------------------------------------------
343 
344  !--- from surface to lowermost atmosphere
345 
346  do j = jsb, jeb
347  do i = isb, ieb
348  dz(ks,i,j) = real_cz(ks,i,j) - real_fz(ks-1,i,j) ! distance from surface to cell center
349  do k = ks+1, ke
350  dz(k,i,j) = real_cz(k,i,j) - real_cz(k-1,i,j) ! distance from cell center to cell center
351  enddo
352  dz(ke+1,i,j) = real_fz(ke,i,j) - real_cz(ke,i,j) ! distance from cell center to TOA
353  enddo
354  enddo
355 
356  do j = jsb, jeb
357  do i = isb, ieb
358  pott_toa(i,j) = pott(ke,i,j)
359  qv_toa(i,j) = qv(ke,i,j)
360  qc_toa(i,j) = qc(ke,i,j)
361  enddo
362  enddo
363 
364  ! density at surface
365  do j = jsb, jeb
366  do i = isb, ieb
367  rtot_sfc(i,j) = rdry * ( 1.0_rp - qv_sfc(1,i,j) - qc_sfc(1,i,j) ) &
368  + rvap * qv_sfc(1,i,j)
369  cvtot_sfc(i,j) = cvdry * ( 1.0_rp - qv_sfc(1,i,j) - qc_sfc(1,i,j) ) &
370  + cv_qv * qv_sfc(1,i,j) &
371  + cv_qc * qc_sfc(1,i,j)
372  cpovcv_sfc(i,j) = ( cvtot_sfc(i,j) + rtot_sfc(i,j) ) / cvtot_sfc(i,j)
373  enddo
374  enddo
375 
376  do j = jsb, jeb
377  do i = isb, ieb
378  rtot(i,j) = rdry * ( 1.0_rp - qv(ks,i,j) - qc(ks,i,j) ) &
379  + rvap * qv(ks,i,j)
380  cvtot(i,j) = cvdry * ( 1.0_rp - qv(ks,i,j) - qc(ks,i,j) ) &
381  + cv_qv * qv(ks,i,j) &
382  + cv_qc * qc(ks,i,j)
383  cpovcv(i,j) = ( cvtot(i,j) + rtot(i,j) ) / cvtot(i,j)
384  enddo
385  enddo
386 
387  ! density at surface
388  do j = jsb, jeb
389  do i = isb, ieb
390  cvovcp_sfc = 1.0_rp / cpovcv_sfc(i,j)
391  dens_sfc(1,i,j) = p00 / rtot_sfc(i,j) / pott_sfc(1,i,j) * ( pres_sfc(1,i,j)/p00 )**cvovcp_sfc
392  temp_sfc(1,i,j) = pres_sfc(1,i,j) / ( dens_sfc(1,i,j) * rtot_sfc(i,j) )
393  enddo
394  enddo
395 
396  ! make density at lowermost cell center
397  if ( hydrostatic_uselapserate ) then
398 
399  do j = jsb, jeb
400  do i = isb, ieb
401  cpovr = ( cvtot(i,j) + rtot(i,j) ) / rtot(i,j)
402  cvovcp = 1.0_rp / cpovcv(i,j)
403 
404  temp(ks,i,j) = pott_sfc(1,i,j) - lapsdry * ( real_cz(ks,i,j) - real_fz(ks-1,i,j) ) ! use dry lapse rate
405  pres(ks,i,j) = p00 * ( temp(ks,i,j)/pott(ks,i,j) )**cpovr
406  dens(ks,i,j) = p00 / rtot(i,j) / pott(ks,i,j) * ( pres(ks,i,j)/p00 )**cvovcp
407  enddo
408  enddo
409 
410  else ! use itelation
411 
412  call atmos_hydrostatic_buildrho_atmos_2d( dens(ks,:,:), & ! [OUT]
413  temp(ks,:,:), & ! [OUT]->not used
414  pres(ks,:,:), & ! [OUT]->not used
415  pott(ks,:,:), & ! [IN]
416  qv(ks,:,:), & ! [IN]
417  qc(ks,:,:), & ! [IN]
418  dens_sfc(1,:,:), & ! [IN]
419  pott_sfc(1,:,:), & ! [IN]
420  qv_sfc(1,:,:), & ! [IN]
421  qc_sfc(1,:,:), & ! [IN]
422  dz(ks,:,:), & ! [IN]
423  ks ) ! [IN]
424 
425  endif
426 
427  !--- from lowermost atmosphere to top of atmosphere
428  call atmos_hydrostatic_buildrho_atmos_3d( dens(:,:,:), & ! [INOUT]
429  temp(:,:,:), & ! [OUT]
430  pres(:,:,:), & ! [OUT]
431  pott(:,:,:), & ! [IN]
432  qv(:,:,:), & ! [IN]
433  qc(:,:,:), & ! [IN]
434  dz(:,:,:) ) ! [IN]
435 
436  call atmos_hydrostatic_buildrho_atmos_2d( dens(ke+1,:,:), & ! [OUT]
437  temp(ke+1,:,:), & ! [OUT]->not used
438  pres(ke+1,:,:), & ! [OUT]->not used
439  pott_toa(:,:), & ! [IN]
440  qv_toa(:,:), & ! [IN]
441  qc_toa(:,:), & ! [IN]
442  dens(ke ,:,:), & ! [IN]
443  pott(ke ,:,:), & ! [IN]
444  qv(ke ,:,:), & ! [IN]
445  qc(ke ,:,:), & ! [IN]
446  dz(ke+1,:,:), & ! [IN]
447  ke+1 ) ! [IN]
448 
449  ! density at TOA
450  dens( 1:ks-1,:,:) = 0.d0 ! fill dummy
451  dens(ke+2:ka ,:,:) = 0.d0 ! fill dummy
452 
453  call comm_horizontal_mean( dens_1d(:), dens(:,:,:) )
454  do j = jsb, jeb
455  do i = isb, ieb
456  dens(:,i,j) = dens_1d(:)
457  enddo
458  enddo
459 
460  call atmos_hydrostatic_buildrho_atmos_rev_2d( dens(ke ,:,:), & ! [OUT]
461  temp(ke ,:,:), & ! [OUT]->not used
462  pres(ke ,:,:), & ! [OUT]->not used
463  pott(ke ,:,:), & ! [IN]
464  qv(ke ,:,:), & ! [IN]
465  qc(ke ,:,:), & ! [IN]
466  dens(ke+1,:,:), & ! [IN]
467  pott_toa(:,:), & ! [IN]
468  qv_toa(:,:), & ! [IN]
469  qc_toa(:,:), & ! [IN]
470  dz(ke+1,:,:), & ! [IN]
471  ke+1 ) ! [IN]
472 
473  !--- from top of atmosphere to lowermost atmosphere
474  call atmos_hydrostatic_buildrho_atmos_rev_3d( dens(:,:,:), & ! [INOUT]
475  temp(:,:,:), & ! [OUT]
476  pres(:,:,:), & ! [OUT]
477  pott(:,:,:), & ! [IN]
478  qv(:,:,:), & ! [IN]
479  qc(:,:,:), & ! [IN]
480  dz(:,:,:) ) ! [IN]
481 
482  return
483  end subroutine atmos_hydrostatic_buildrho_3d
484 
485  !-----------------------------------------------------------------------------
487  subroutine atmos_hydrostatic_buildrho_real_3d( &
488  dens, &
489  temp, &
490  pres, &
491  pott, &
492  qv, &
493  qc )
494  use scale_process, only: &
495  prc_mpistop
496  implicit none
497 
498  real(RP), intent(out) :: dens(ka,ia,ja)
499  real(RP), intent(out) :: temp(ka,ia,ja)
500  real(RP), intent(inout) :: pres(ka,ia,ja)
501  real(RP), intent(in) :: pott(ka,ia,ja)
502  real(RP), intent(in) :: qv (ka,ia,ja)
503  real(RP), intent(in) :: qc (ka,ia,ja)
504 
505  real(RP) :: dz(ka,ia,ja)
506 
507  real(RP) :: pott_toa(ia,ja)
508  real(RP) :: qv_toa (ia,ja)
509  real(RP) :: qc_toa (ia,ja)
510 
511  real(RP) :: Rtot
512  real(RP) :: CVtot
513  real(RP) :: CVovCP
514 
515  integer :: kref
516  integer :: k, i, j
517  !---------------------------------------------------------------------------
518 
519  !--- from surface to lowermost atmosphere
520 
521  do j = jsb, jeb
522  do i = isb, ieb
523  dz(ks,i,j) = real_cz(ks,i,j) - real_fz(ks-1,i,j) ! distance from surface to cell center
524  do k = ks+1, ke
525  dz(k,i,j) = real_cz(k,i,j) - real_cz(k-1,i,j) ! distance from cell center to cell center
526  enddo
527  dz(ke+1,i,j) = real_fz(ke,i,j) - real_cz(ke,i,j) ! distance from cell center to TOA
528  enddo
529  enddo
530 
531  do j = jsb, jeb
532  do i = isb, ieb
533  pott_toa(i,j) = pott(ke,i,j)
534  qv_toa(i,j) = qv(ke,i,j)
535  qc_toa(i,j) = qc(ke,i,j)
536  enddo
537  enddo
538 
539  kref = hydrostatic_buildrho_real_kref + ks - 1
540 
541  ! calc density at reference level
542  do j = jsb, jeb
543  do i = isb, ieb
544  rtot = rdry * ( 1.0_rp - qv(kref,i,j) - qc(kref,i,j) ) &
545  + rvap * qv(kref,i,j)
546  cvtot = cvdry * ( 1.0_rp - qv(kref,i,j) - qc(kref,i,j) ) &
547  + cv_qv * qv(kref,i,j) &
548  + cv_qc * qc(kref,i,j)
549  cvovcp = cvtot / ( cvtot + rtot )
550  dens(kref,i,j) = p00 / ( rtot * pott(kref,i,j) ) * ( pres(kref,i,j)/p00 )**cvovcp
551  enddo
552  enddo
553 
554  !--- from lowermost atmosphere to top of atmosphere
555  call atmos_hydrostatic_buildrho_atmos_3d( dens(:,:,:), & ! [INOUT]
556  temp(:,:,:), & ! [OUT]
557  pres(:,:,:), & ! [OUT]
558  pott(:,:,:), & ! [IN]
559  qv(:,:,:), & ! [IN]
560  qc(:,:,:), & ! [IN]
561  dz(:,:,:), & ! [IN]
562  kref ) ! [IN]
563  call atmos_hydrostatic_buildrho_atmos_rev_3d( dens(:,:,:), & ! [INOUT]
564  temp(:,:,:), & ! [OUT]
565  pres(:,:,:), & ! [OUT]
566  pott(:,:,:), & ! [IN]
567  qv(:,:,:), & ! [IN]
568  qc(:,:,:), & ! [IN]
569  dz(:,:,:), & ! [IN]
570  kref ) ! [IN]
571 
572  call atmos_hydrostatic_buildrho_atmos_2d( dens(ke+1,:,:), & ! [OUT]
573  temp(ke+1,:,:), & ! [OUT]->not used
574  pres(ke+1,:,:), & ! [OUT]->not used
575  pott_toa(:,:), & ! [IN]
576  qv_toa(:,:), & ! [IN]
577  qc_toa(:,:), & ! [IN]
578  dens(ke ,:,:), & ! [IN]
579  pott(ke ,:,:), & ! [IN]
580  qv(ke ,:,:), & ! [IN]
581  qc(ke ,:,:), & ! [IN]
582  dz(ke+1,:,:), & ! [IN]
583  ke+1 ) ! [IN]
584 
585  ! density at TOA
586  dens( 1:ks-1,:,:) = 0.d0 ! fill dummy
587  dens(ke+2:ka ,:,:) = 0.d0 ! fill dummy
588 
589  return
590  end subroutine atmos_hydrostatic_buildrho_real_3d
591 
592  !-----------------------------------------------------------------------------
594  subroutine atmos_hydrostatic_buildrho_atmos_0d( &
595  dens_L2, &
596  temp_L2, &
597  pres_L2, &
598  pott_L2, &
599  qv_L2, &
600  qc_L2, &
601  dens_L1, &
602  pott_L1, &
603  qv_L1, &
604  qc_L1, &
605  dz, &
606  k )
607  use scale_process, only: &
608  prc_mpistop
609  implicit none
610 
611  real(RP), intent(out) :: dens_L2
612  real(RP), intent(out) :: temp_L2
613  real(RP), intent(out) :: pres_L2
614  real(RP), intent(in) :: pott_L2
615  real(RP), intent(in) :: qv_L2
616  real(RP), intent(in) :: qc_L2
617  real(RP), intent(in) :: dens_L1
618  real(RP), intent(in) :: pott_L1
619  real(RP), intent(in) :: qv_L1
620  real(RP), intent(in) :: qc_L1
621  real(RP), intent(in) :: dz
622  integer, intent(in) :: k
623 
624  real(RP) :: Rtot_L1 , Rtot_L2
625  real(RP) :: CVtot_L1 , CVtot_L2
626  real(RP) :: CPovCV_L1, CPovCV_L2
627 
628  real(RP) :: RovCV
629  real(RP) :: dens_s, dhyd, dgrd
630  integer :: ite
631  logical :: converged
632  !---------------------------------------------------------------------------
633 
634  rtot_l1 = rdry * ( 1.0_rp - qv_l1 - qc_l1 ) &
635  + rvap * qv_l1
636  cvtot_l1 = cvdry * ( 1.0_rp - qv_l1 - qc_l1 ) &
637  + cv_qv * qv_l1 &
638  + cv_qc * qc_l1
639  cpovcv_l1 = ( cvtot_l1 + rtot_l1 ) / cvtot_l1
640 
641  rtot_l2 = rdry * ( 1.0_rp - qv_l2 - qc_l2 ) &
642  + rvap * qv_l2
643  cvtot_l2 = cvdry * ( 1.0_rp - qv_l2 - qc_l2 ) &
644  + cv_qv * qv_l2 &
645  + cv_qc * qc_l2
646  cpovcv_l2 = ( cvtot_l2 + rtot_l2 ) / cvtot_l2
647 
648  rovcv = rtot_l2 / cvtot_l2
649 
650  dens_s = 0.0_rp
651  dens_l2 = dens_l1 ! first guess
652 
653  converged = .false.
654  do ite = 1, itelim
655  if ( abs(dens_l2-dens_s) <= criteria ) then
656  converged = .true.
657  exit
658  endif
659 
660  dens_s = dens_l2
661 
662  dhyd = + ( p00 * ( dens_l1 * rtot_l1 * pott_l1 / p00 )**cpovcv_l1 &
663  - p00 * ( dens_s * rtot_l2 * pott_l2 / p00 )**cpovcv_l2 ) / dz & ! dp/dz
664  - grav * 0.5_rp * ( dens_l1 + dens_s ) ! rho*g
665 
666  dgrd = - p00 * ( rtot_l2 * pott_l2 / p00 )**cpovcv_l2 / dz &
667  * cpovcv_l2 * dens_s**rovcv &
668  - 0.5_rp * grav
669 
670  dens_l2 = dens_s - dhyd/dgrd
671 
672  if( dens_l2*0.0_rp /= 0.0_rp) exit
673  enddo
674 
675  if ( .NOT. converged ) then
676  write(*,*) 'xxx [buildrho 0D atmos] iteration not converged!', &
677  k,dens_l2,ite,dens_s,dhyd,dgrd
678  call prc_mpistop
679  endif
680 
681  pres_l2 = p00 * ( dens_l2 * rtot_l2 * pott_l2 / p00 )**cpovcv_l2
682  temp_l2 = pres_l2 / ( dens_l2 * rtot_l2 )
683 
684  return
685  end subroutine atmos_hydrostatic_buildrho_atmos_0d
686 
687  !-----------------------------------------------------------------------------
689  subroutine atmos_hydrostatic_buildrho_atmos_1d( &
690  dens, &
691  temp, &
692  pres, &
693  pott, &
694  qv, &
695  qc )
696  use scale_process, only: &
697  prc_mpistop
698  implicit none
699 
700  real(RP), intent(inout) :: dens(ka)
701  real(RP), intent(out) :: temp(ka)
702  real(RP), intent(out) :: pres(ka)
703  real(RP), intent(in) :: pott(ka)
704  real(RP), intent(in) :: qv (ka)
705  real(RP), intent(in) :: qc (ka)
706 
707  real(RP) :: Rtot (ka)
708  real(RP) :: CVtot (ka)
709  real(RP) :: CPovCV(ka)
710 
711  real(RP) :: RovCV
712  real(RP) :: dens_s, dhyd, dgrd
713  integer :: ite
714  logical :: converged
715 
716  integer :: k
717  !---------------------------------------------------------------------------
718 
719  do k = ks, ke
720  rtot(k) = rdry * ( 1.0_rp - qv(k) - qc(k) ) &
721  + rvap * qv(k)
722  cvtot(k) = cvdry * ( 1.0_rp - qv(k) - qc(k) ) &
723  + cv_qv * qv(k) &
724  + cv_qc * qc(k)
725  cpovcv(k) = ( cvtot(k) + rtot(k) ) / cvtot(k)
726  enddo
727 
728  do k = ks+1, ke
729  rovcv = rtot(k) / cvtot(k)
730 
731  dens_s = 0.0_rp
732  dens(k) = dens(k-1) ! first guess
733 
734  converged = .false.
735  do ite = 1, itelim
736  if ( abs(dens(k)-dens_s) <= criteria ) then
737  converged = .true.
738  exit
739  endif
740 
741  dens_s = dens(k)
742 
743  dhyd = + ( p00 * ( dens(k-1) * rtot(k-1) * pott(k-1) / p00 )**cpovcv(k-1) &
744  - p00 * ( dens_s * rtot(k ) * pott(k ) / p00 )**cpovcv(k ) ) / grid_fdz(k-1) & ! dpdz
745  - grav * 0.5_rp * ( dens(k-1) + dens_s ) ! rho*g
746 
747  dgrd = - p00 * ( rtot(k) * pott(k) / p00 )**cpovcv(k) / grid_fdz(k-1) &
748  * cpovcv(k) * dens_s**rovcv &
749  - 0.5_rp * grav
750 
751  dens(k) = dens_s - dhyd/dgrd
752 
753  if( dens(k)*0.0_rp /= 0.0_rp) exit
754  enddo
755 
756  if ( .NOT. converged ) then
757  write(*,*) 'xxx [buildrho 1D atmos] iteration not converged!', &
758  k,dens(k),ite,dens_s,dhyd,dgrd
759  call prc_mpistop
760  endif
761  enddo
762 
763  do k = ks, ke
764  pres(k) = p00 * ( dens(k) * rtot(k) * pott(k) / p00 )**cpovcv(k)
765  temp(k) = pres(k) / ( dens(k) * rtot(k) )
766  enddo
767 
768  dens( 1:ks-1) = dens(ks)
769  dens(ke+1:ka ) = dens(ke)
770  pres( 1:ks-1) = pres(ks)
771  pres(ke+1:ka ) = pres(ke)
772  temp( 1:ks-1) = temp(ks)
773  temp(ke+1:ka ) = temp(ke)
774 
775  return
776  end subroutine atmos_hydrostatic_buildrho_atmos_1d
777 
778  !-----------------------------------------------------------------------------
780  subroutine atmos_hydrostatic_buildrho_atmos_2d( &
781  dens_L2, &
782  temp_L2, &
783  pres_L2, &
784  pott_L2, &
785  qv_L2, &
786  qc_L2, &
787  dens_L1, &
788  pott_L1, &
789  qv_L1, &
790  qc_L1, &
791  dz, &
792  k )
793  use scale_process, only: &
794  prc_mpistop
795  implicit none
796 
797  real(RP), intent(out) :: dens_L2(ia,ja)
798  real(RP), intent(out) :: temp_L2(ia,ja)
799  real(RP), intent(out) :: pres_L2(ia,ja)
800  real(RP), intent(in) :: pott_L2(ia,ja)
801  real(RP), intent(in) :: qv_L2 (ia,ja)
802  real(RP), intent(in) :: qc_L2 (ia,ja)
803  real(RP), intent(in) :: dens_L1(ia,ja)
804  real(RP), intent(in) :: pott_L1(ia,ja)
805  real(RP), intent(in) :: qv_L1 (ia,ja)
806  real(RP), intent(in) :: qc_L1 (ia,ja)
807  real(RP), intent(in) :: dz (ia,ja)
808  integer, intent(in) :: k
809 
810  real(RP) :: Rtot_L1 (ia,ja), Rtot_L2 (ia,ja)
811  real(RP) :: CVtot_L1 (ia,ja), CVtot_L2 (ia,ja)
812  real(RP) :: CPovCV_L1(ia,ja), CPovCV_L2(ia,ja)
813 
814  real(RP) :: RovCV
815  real(RP) :: dens_s, dhyd, dgrd
816  integer :: ite
817  logical :: converged
818 
819  integer :: i, j
820  !---------------------------------------------------------------------------
821 
822  do j = jsb, jeb
823  do i = isb, ieb
824  rtot_l1(i,j) = rdry * ( 1.0_rp - qv_l1(i,j) - qc_l1(i,j) ) &
825  + rvap * qv_l1(i,j)
826  cvtot_l1(i,j) = cvdry * ( 1.0_rp - qv_l1(i,j) - qc_l1(i,j) ) &
827  + cv_qv * qv_l1(i,j) &
828  + cv_qc * qc_l1(i,j)
829  cpovcv_l1(i,j) = ( cvtot_l1(i,j) + rtot_l1(i,j) ) / cvtot_l1(i,j)
830 
831  rtot_l2(i,j) = rdry * ( 1.0_rp - qv_l2(i,j) - qc_l2(i,j) ) &
832  + rvap * qv_l2(i,j)
833  cvtot_l2(i,j) = cvdry * ( 1.0_rp - qv_l2(i,j) - qc_l2(i,j) ) &
834  + cv_qv * qv_l2(i,j) &
835  + cv_qc * qc_l2(i,j)
836  cpovcv_l2(i,j) = ( cvtot_l2(i,j) + rtot_l2(i,j) ) / cvtot_l2(i,j)
837  enddo
838  enddo
839 
840  do j = jsb, jeb
841  do i = isb, ieb
842  rovcv = rtot_l2(i,j) / cvtot_l2(i,j)
843 
844  dens_s = 0.0_rp
845  dens_l2(i,j) = dens_l1(i,j) ! first guess
846 
847  converged = .false.
848  do ite = 1, itelim
849  if ( abs(dens_l2(i,j)-dens_s) <= criteria ) then
850  converged = .true.
851  exit
852  endif
853 
854  dens_s = dens_l2(i,j)
855 
856  dhyd = + ( p00 * ( dens_l1(i,j) * rtot_l1(i,j) * pott_l1(i,j) / p00 )**cpovcv_l1(i,j) &
857  - p00 * ( dens_s * rtot_l2(i,j) * pott_l2(i,j) / p00 )**cpovcv_l2(i,j) ) / dz(i,j) & ! dp/dz
858  - grav * 0.5_rp * ( dens_l1(i,j) + dens_s ) ! rho*g
859 
860  dgrd = - p00 * ( rtot_l2(i,j) * pott_l2(i,j) / p00 )**cpovcv_l2(i,j) / dz(i,j) &
861  * cpovcv_l2(i,j) * dens_s**rovcv &
862  - 0.5_rp * grav
863 
864  dens_l2(i,j) = dens_s - dhyd/dgrd
865 
866  if( dens_l2(i,j)*0.0_rp /= 0.0_rp) exit
867  enddo
868 
869  if ( .NOT. converged ) then
870  write(*,*) 'xxx [buildrho 2D atmos] iteration not converged!', &
871  i,j,k,dens_l2(i,j),ite,dens_s,dhyd,dgrd
872  call prc_mpistop
873  endif
874  enddo
875  enddo
876 
877  do j = jsb, jeb
878  do i = isb, ieb
879  pres_l2(i,j) = p00 * ( dens_l2(i,j) * rtot_l2(i,j) * pott_l2(i,j) / p00 )**cpovcv_l2(i,j)
880  temp_l2(i,j) = pres_l2(i,j) / ( dens_l2(i,j) * rtot_l2(i,j) )
881  enddo
882  enddo
883 
884  return
885  end subroutine atmos_hydrostatic_buildrho_atmos_2d
886 
887  !-----------------------------------------------------------------------------
889  subroutine atmos_hydrostatic_buildrho_atmos_3d( &
890  dens, &
891  temp, &
892  pres, &
893  pott, &
894  qv, &
895  qc, &
896  dz, &
897  kref_in )
898  use scale_process, only: &
899  prc_mpistop
900  implicit none
901 
902  real(RP), intent(inout) :: dens(ka,ia,ja)
903  real(RP), intent(out) :: temp(ka,ia,ja)
904  real(RP), intent(out) :: pres(ka,ia,ja)
905  real(RP), intent(in) :: pott(ka,ia,ja)
906  real(RP), intent(in) :: qv (ka,ia,ja)
907  real(RP), intent(in) :: qc (ka,ia,ja)
908  real(RP), intent(in) :: dz (ka,ia,ja)
909  integer, intent(in), optional :: kref_in
910 
911  real(RP) :: Rtot (ka,ia,ja)
912  real(RP) :: CVtot (ka,ia,ja)
913  real(RP) :: CPovCV(ka,ia,ja)
914 
915  real(RP) :: RovCV
916  real(RP) :: dens_s, dhyd, dgrd
917  integer :: ite
918  logical :: converged
919 
920  integer :: kref
921  integer :: k, i, j
922  !---------------------------------------------------------------------------
923 
924  if ( present(kref_in) ) then
925  kref = kref_in
926  else
927  kref = ks
928  end if
929 
930  do j = jsb, jeb
931  do i = isb, ieb
932  do k = kref, ke
933  rtot(k,i,j) = rdry * ( 1.0_rp - qv(k,i,j) - qc(k,i,j) ) &
934  + rvap * qv(k,i,j)
935  cvtot(k,i,j) = cvdry * ( 1.0_rp - qv(k,i,j) - qc(k,i,j) ) &
936  + cv_qv * qv(k,i,j) &
937  + cv_qc * qc(k,i,j)
938  cpovcv(k,i,j) = ( cvtot(k,i,j) + rtot(k,i,j) ) / cvtot(k,i,j)
939  enddo
940  enddo
941  enddo
942 
943  do j = jsb, jeb
944  do i = isb, ieb
945  do k = kref+1, ke
946  rovcv = rtot(k,i,j) / cvtot(k,i,j)
947 
948  dens_s = 0.0_rp
949  dens(k,i,j) = dens(k-1,i,j) ! first guess
950 
951  converged = .false.
952  do ite = 1, itelim
953  if ( abs(dens(k,i,j)-dens_s) <= criteria ) then
954  converged = .true.
955  exit
956  endif
957 
958  dens_s = dens(k,i,j)
959 
960  dhyd = + ( p00 * ( dens(k-1,i,j) * rtot(k-1,i,j) * pott(k-1,i,j) / p00 )**cpovcv(k-1,i,j) &
961  - p00 * ( dens_s * rtot(k ,i,j) * pott(k ,i,j) / p00 )**cpovcv(k ,i,j) ) / dz(k,i,j) & ! dpdz
962  - grav * 0.5_rp * ( dens(k-1,i,j) + dens_s ) ! rho*g
963 
964  dgrd = - p00 * ( rtot(k,i,j) * pott(k,i,j) / p00 )**cpovcv(k,i,j) / dz(k,i,j) &
965  * cpovcv(k,i,j) * dens_s**rovcv &
966  - 0.5_rp * grav
967 
968  dens(k,i,j) = dens_s - dhyd/dgrd
969 
970  if( dens(k,i,j)*0.0_rp /= 0.0_rp) exit
971  enddo
972 
973  if ( .NOT. converged ) then
974  write(*,*) 'xxx [buildrho 3D atmos] iteration not converged!', &
975  k,i,j,dens(k,i,j),ite,dens_s,dhyd,dgrd
976  call prc_mpistop
977  endif
978  enddo
979  enddo
980  enddo
981 
982  do j = jsb, jeb
983  do i = isb, ieb
984  do k = kref, ke
985  pres(k,i,j) = p00 * ( dens(k,i,j) * rtot(k,i,j) * pott(k,i,j) / p00 )**cpovcv(k,i,j)
986  temp(k,i,j) = pres(k,i,j) / ( dens(k,i,j) * rtot(k,i,j) )
987  enddo
988  enddo
989  enddo
990 
991  return
992  end subroutine atmos_hydrostatic_buildrho_atmos_3d
993 
994  !-----------------------------------------------------------------------------
996  subroutine atmos_hydrostatic_buildrho_atmos_rev_2d( &
997  dens_L1, &
998  temp_L1, &
999  pres_L1, &
1000  pott_L1, &
1001  qv_L1, &
1002  qc_L1, &
1003  dens_L2, &
1004  pott_L2, &
1005  qv_L2, &
1006  qc_L2, &
1007  dz, &
1008  k )
1009  use scale_process, only: &
1010  prc_mpistop
1011  implicit none
1012 
1013  real(RP), intent(out) :: dens_L1(ia,ja)
1014  real(RP), intent(out) :: temp_L1(ia,ja)
1015  real(RP), intent(out) :: pres_L1(ia,ja)
1016  real(RP), intent(in) :: pott_L1(ia,ja)
1017  real(RP), intent(in) :: qv_L1 (ia,ja)
1018  real(RP), intent(in) :: qc_L1 (ia,ja)
1019  real(RP), intent(in) :: dens_L2(ia,ja)
1020  real(RP), intent(in) :: pott_L2(ia,ja)
1021  real(RP), intent(in) :: qv_L2 (ia,ja)
1022  real(RP), intent(in) :: qc_L2 (ia,ja)
1023  real(RP), intent(in) :: dz (ia,ja)
1024  integer, intent(in) :: k
1025 
1026  real(RP) :: Rtot_L1 (ia,ja), Rtot_L2 (ia,ja)
1027  real(RP) :: CVtot_L1 (ia,ja), CVtot_L2 (ia,ja)
1028  real(RP) :: CPovCV_L1(ia,ja), CPovCV_L2(ia,ja)
1029 
1030  real(RP) :: RovCV
1031  real(RP) :: dens_s, dhyd, dgrd
1032  integer :: ite
1033  logical :: converged
1034 
1035  integer :: i, j
1036  !---------------------------------------------------------------------------
1037 
1038  do j = jsb, jeb
1039  do i = isb, ieb
1040  rtot_l1(i,j) = rdry * ( 1.0_rp - qv_l1(i,j) - qc_l1(i,j) ) &
1041  + rvap * qv_l1(i,j)
1042  cvtot_l1(i,j) = cvdry * ( 1.0_rp - qv_l1(i,j) - qc_l1(i,j) ) &
1043  + cv_qv * qv_l1(i,j) &
1044  + cv_qc * qc_l1(i,j)
1045  cpovcv_l1(i,j) = ( cvtot_l1(i,j) + rtot_l1(i,j) ) / cvtot_l1(i,j)
1046 
1047  rtot_l2(i,j) = rdry * ( 1.0_rp - qv_l2(i,j) - qc_l2(i,j) ) &
1048  + rvap * qv_l2(i,j)
1049  cvtot_l2(i,j) = cvdry * ( 1.0_rp - qv_l2(i,j) - qc_l2(i,j) ) &
1050  + cv_qv * qv_l2(i,j) &
1051  + cv_qc * qc_l2(i,j)
1052  cpovcv_l2(i,j) = ( cvtot_l2(i,j) + rtot_l2(i,j) ) / cvtot_l2(i,j)
1053  enddo
1054  enddo
1055 
1056  do j = jsb, jeb
1057  do i = isb, ieb
1058  rovcv = rtot_l1(i,j) / cvtot_l1(i,j)
1059 
1060  dens_s = 0.0_rp
1061  dens_l1(i,j) = dens_l2(i,j) ! first guess
1062 
1063  converged = .false.
1064  do ite = 1, itelim
1065  if ( abs(dens_l1(i,j)-dens_s) <= criteria ) then
1066  converged = .true.
1067  exit
1068  endif
1069 
1070  dens_s = dens_l1(i,j)
1071 
1072  dhyd = + ( p00 * ( dens_s * rtot_l1(i,j) * pott_l1(i,j) / p00 )**cpovcv_l1(i,j) &
1073  - p00 * ( dens_l2(i,j) * rtot_l2(i,j) * pott_l2(i,j) / p00 )**cpovcv_l2(i,j) ) / dz(i,j) & ! dp/dz
1074  - grav * 0.5_rp * ( dens_s + dens_l2(i,j) ) ! rho*g
1075 
1076  dgrd = + p00 * ( rtot_l1(i,j) * pott_l1(i,j) / p00 )**cpovcv_l1(i,j) / dz(i,j) &
1077  * cpovcv_l1(i,j) * dens_s**rovcv &
1078  - 0.5_rp * grav
1079 
1080  dens_l1(i,j) = dens_s - dhyd/dgrd
1081 
1082  if( dens_l1(i,j)*0.0_rp /= 0.0_rp) exit
1083  enddo
1084 
1085  if ( .NOT. converged ) then
1086  write(*,*) 'xxx [buildrho 2D rev atmos] iteration not converged!', &
1087  i,j,k,dens_l1(i,j),ite,dens_s,dhyd,dgrd
1088  call prc_mpistop
1089  endif
1090  enddo
1091  enddo
1092 
1093  do j = jsb, jeb
1094  do i = isb, ieb
1095  pres_l1(i,j) = p00 * ( dens_l1(i,j) * rtot_l1(i,j) * pott_l1(i,j) / p00 )**cpovcv_l1(i,j)
1096  temp_l1(i,j) = pres_l1(i,j) / ( dens_l1(i,j) * rtot_l1(i,j) )
1097  enddo
1098  enddo
1099 
1100  return
1101  end subroutine atmos_hydrostatic_buildrho_atmos_rev_2d
1102 
1103  !-----------------------------------------------------------------------------
1105  subroutine atmos_hydrostatic_buildrho_atmos_rev_3d( &
1106  dens, &
1107  temp, &
1108  pres, &
1109  pott, &
1110  qv, &
1111  qc, &
1112  dz, &
1113  kref_in )
1114  use scale_process, only: &
1115  prc_mpistop
1116  implicit none
1117 
1118  real(RP), intent(inout) :: dens(ka,ia,ja)
1119  real(RP), intent(out) :: temp(ka,ia,ja)
1120  real(RP), intent(out) :: pres(ka,ia,ja)
1121  real(RP), intent(in) :: pott(ka,ia,ja)
1122  real(RP), intent(in) :: qv (ka,ia,ja)
1123  real(RP), intent(in) :: qc (ka,ia,ja)
1124  real(RP), intent(in) :: dz (ka,ia,ja)
1125  integer, intent(in), optional :: kref_in
1126 
1127  real(RP) :: Rtot (ka,ia,ja)
1128  real(RP) :: CVtot (ka,ia,ja)
1129  real(RP) :: CPovCV(ka,ia,ja)
1130 
1131  real(RP) :: RovCV
1132  real(RP) :: dens_s, dhyd, dgrd
1133  integer :: ite
1134  logical :: converged
1135 
1136  integer :: kref
1137  integer :: k, i, j
1138  !---------------------------------------------------------------------------
1139 
1140  if ( present(kref_in) ) then
1141  kref = kref_in
1142  else
1143  kref = ke
1144  end if
1145 
1146  do j = jsb, jeb
1147  do i = isb, ieb
1148  do k = ks, kref
1149  rtot(k,i,j) = rdry * ( 1.0_rp - qv(k,i,j) - qc(k,i,j) ) &
1150  + rvap * qv(k,i,j)
1151  cvtot(k,i,j) = cvdry * ( 1.0_rp - qv(k,i,j) - qc(k,i,j) ) &
1152  + cv_qv * qv(k,i,j) &
1153  + cv_qc * qc(k,i,j)
1154  cpovcv(k,i,j) = ( cvtot(k,i,j) + rtot(k,i,j) ) / cvtot(k,i,j)
1155  enddo
1156  enddo
1157  enddo
1158 
1159  do j = jsb, jeb
1160  do i = isb, ieb
1161  do k = kref-1, ks, -1
1162  rovcv = rtot(k,i,j) / cvtot(k,i,j)
1163 
1164  dens_s = 0.0_rp
1165  dens(k,i,j) = dens(k+1,i,j) ! first guess
1166 
1167  converged = .false.
1168  do ite = 1, itelim
1169  if ( abs(dens(k,i,j)-dens_s) <= criteria ) then
1170  converged = .true.
1171  exit
1172  endif
1173 
1174  dens_s = dens(k,i,j)
1175 
1176  dhyd = + ( p00 * ( dens_s * rtot(k ,i,j) * pott(k ,i,j) / p00 )**cpovcv(k ,i,j) &
1177  - p00 * ( dens(k+1,i,j) * rtot(k+1,i,j) * pott(k+1,i,j) / p00 )**cpovcv(k+1,i,j) ) / dz(k+1,i,j) & ! dpdz
1178  - grav * 0.5_rp * ( dens_s + dens(k+1,i,j) ) ! rho*g
1179 
1180  dgrd = + p00 * ( rtot(k,i,j) * pott(k,i,j) / p00 )**cpovcv(k,i,j) / dz(k+1,i,j) &
1181  * cpovcv(k,i,j) * dens_s**rovcv &
1182  - 0.5_rp * grav
1183 
1184  dens(k,i,j) = dens_s - dhyd/dgrd
1185 
1186  if( dens(k,i,j)*0.0_rp /= 0.0_rp) exit
1187  enddo
1188 
1189  if ( .NOT. converged ) then
1190  write(*,*) 'xxx [buildrho 3D rev atmos] iteration not converged!', &
1191  k,i,j,dens(k,i,j),ite,dens_s,dhyd,dgrd
1192  call prc_mpistop
1193  endif
1194  enddo
1195  enddo
1196  enddo
1197 
1198  do j = jsb, jeb
1199  do i = isb, ieb
1200  do k = ks, kref
1201  pres(k,i,j) = p00 * ( dens(k,i,j) * rtot(k,i,j) * pott(k,i,j) / p00 )**cpovcv(k,i,j)
1202  temp(k,i,j) = pres(k,i,j) / ( dens(k,i,j) * rtot(k,i,j) )
1203  enddo
1204  enddo
1205  enddo
1206 
1207  return
1208  end subroutine atmos_hydrostatic_buildrho_atmos_rev_3d
1209 
1210  !-----------------------------------------------------------------------------
1212  subroutine atmos_hydrostatic_buildrho_bytemp_1d( &
1213  dens, &
1214  pott, &
1215  pres, &
1216  temp, &
1217  qv, &
1218  qc, &
1219  pott_sfc, &
1220  pres_sfc, &
1221  temp_sfc, &
1222  qv_sfc, &
1223  qc_sfc )
1224  use scale_process, only: &
1225  prc_mpistop
1226  implicit none
1227 
1228  real(RP), intent(out) :: dens(ka)
1229  real(RP), intent(out) :: pott(ka)
1230  real(RP), intent(out) :: pres(ka)
1231  real(RP), intent(in) :: temp(ka)
1232  real(RP), intent(in) :: qv (ka)
1233  real(RP), intent(in) :: qc (ka)
1234 
1235  real(RP), intent(out) :: pott_sfc
1236  real(RP), intent(in) :: pres_sfc
1237  real(RP), intent(in) :: temp_sfc
1238  real(RP), intent(in) :: qv_sfc
1239  real(RP), intent(in) :: qc_sfc
1240 
1241  real(RP) :: dens_sfc
1242 
1243  real(RP) :: Rtot_sfc
1244  real(RP) :: CVtot_sfc
1245  real(RP) :: Rtot
1246  real(RP) :: CVtot
1247 
1248  real(RP) :: RovCP_sfc
1249  real(RP) :: dens_s, dhyd, dgrd
1250  integer :: ite
1251  logical :: converged
1252  !---------------------------------------------------------------------------
1253 
1254  !--- from surface to lowermost atmosphere
1255 
1256  rtot_sfc = rdry * ( 1.0_rp - qv_sfc - qc_sfc ) &
1257  + rvap * qv_sfc
1258  cvtot_sfc = cvdry * ( 1.0_rp - qv_sfc - qc_sfc ) &
1259  + cv_qv * qv_sfc &
1260  + cv_qc * qc_sfc
1261 
1262  rtot = rdry * ( 1.0_rp - qv(ks) - qc(ks) ) &
1263  + rvap * qv(ks)
1264  cvtot = cvdry * ( 1.0_rp - qv(ks) - qc(ks) ) &
1265  + cv_qv * qv(ks) &
1266  + cv_qc * qc(ks)
1267 
1268  ! density at surface
1269  rovcp_sfc = rtot_sfc / ( cvtot_sfc + rtot_sfc )
1270  dens_sfc = pres_sfc / ( rtot_sfc * temp_sfc )
1271  pott_sfc = temp_sfc * ( p00/pres_sfc )**rovcp_sfc
1272 
1273  ! make density at lowermost cell center
1274  dens_s = 0.0_rp
1275  dens(ks) = dens_sfc ! first guess
1276 
1277  converged = .false.
1278  do ite = 1, itelim
1279  if ( abs(dens(ks)-dens_s) <= criteria ) then
1280  converged = .true.
1281  exit
1282  endif
1283 
1284  dens_s = dens(ks)
1285 
1286  dhyd = + ( dens_sfc * rtot_sfc * temp_sfc &
1287  - dens_s * rtot * temp(ks) ) / grid_cz(ks) & ! dp/dz
1288  - grav * 0.5_rp * ( dens_sfc + dens_s ) ! rho*g
1289 
1290  dgrd = - rtot * temp(ks) / grid_cz(ks) &
1291  - 0.5_rp * grav
1292 
1293  dens(ks) = dens_s - dhyd/dgrd
1294 
1295  if( dens(ks)*0.0_rp /= 0.0_rp) exit
1296  enddo
1297 
1298  if ( .NOT. converged ) then
1299  write(*,*) 'xxx [buildrho bytemp 1D sfc] iteration not converged!', &
1300  dens(ks),ite,dens_s,dhyd,dgrd
1301  call prc_mpistop
1302  endif
1303 
1304  !--- from lowermost atmosphere to top of atmosphere
1305  call atmos_hydrostatic_buildrho_bytemp_atmos_1d( dens(:), & ! [INOUT]
1306  pott(:), & ! [OUT]
1307  pres(:), & ! [OUT]
1308  temp(:), & ! [IN]
1309  qv(:), & ! [IN]
1310  qc(:) ) ! [IN]
1311 
1312  return
1313  end subroutine atmos_hydrostatic_buildrho_bytemp_1d
1314 
1315  !-----------------------------------------------------------------------------
1317  subroutine atmos_hydrostatic_buildrho_bytemp_3d( &
1318  dens, &
1319  pott, &
1320  pres, &
1321  temp, &
1322  qv, &
1323  qc, &
1324  pott_sfc, &
1325  pres_sfc, &
1326  temp_sfc, &
1327  qv_sfc, &
1328  qc_sfc )
1329  use scale_process, only: &
1330  prc_mpistop
1331  implicit none
1332 
1333  real(RP), intent(out) :: dens(ka,ia,ja)
1334  real(RP), intent(out) :: pott(ka,ia,ja)
1335  real(RP), intent(out) :: pres(ka,ia,ja)
1336  real(RP), intent(in) :: temp(ka,ia,ja)
1337  real(RP), intent(in) :: qv (ka,ia,ja)
1338  real(RP), intent(in) :: qc (ka,ia,ja)
1339 
1340  real(RP), intent(out) :: pott_sfc(1,ia,ja)
1341  real(RP), intent(in) :: pres_sfc(1,ia,ja)
1342  real(RP), intent(in) :: temp_sfc(1,ia,ja)
1343  real(RP), intent(in) :: qv_sfc (1,ia,ja)
1344  real(RP), intent(in) :: qc_sfc (1,ia,ja)
1345 
1346  real(RP) :: dens_sfc (1,ia,ja)
1347 
1348  real(RP) :: Rtot_sfc (ia,ja)
1349  real(RP) :: CVtot_sfc (ia,ja)
1350  real(RP) :: Rtot (ia,ja)
1351  real(RP) :: CVtot (ia,ja)
1352 
1353  real(RP) :: RovCP_sfc
1354  real(RP) :: DZ
1355  real(RP) :: dens_s, dhyd, dgrd
1356  integer :: ite
1357  logical :: converged
1358 
1359  integer :: i, j
1360  !---------------------------------------------------------------------------
1361 
1362  !--- from surface to lowermost atmosphere
1363 
1364  do j = jsb, jeb
1365  do i = isb, ieb
1366  rtot_sfc(i,j) = rdry * ( 1.0_rp - qv_sfc(1,i,j) - qc_sfc(1,i,j) ) &
1367  + rvap * qv_sfc(1,i,j)
1368  cvtot_sfc(i,j) = cvdry * ( 1.0_rp - qv_sfc(1,i,j) - qc_sfc(1,i,j) ) &
1369  + cv_qv * qv_sfc(1,i,j) &
1370  + cv_qc * qc_sfc(1,i,j)
1371  enddo
1372  enddo
1373 
1374  do j = jsb, jeb
1375  do i = isb, ieb
1376  rtot(i,j) = rdry * ( 1.0_rp - qv(ks,i,j) - qc(ks,i,j) ) &
1377  + rvap * qv(ks,i,j)
1378  cvtot(i,j) = cvdry * ( 1.0_rp - qv(ks,i,j) - qc(ks,i,j) ) &
1379  + cv_qv * qv(ks,i,j) &
1380  + cv_qc * qc(ks,i,j)
1381  enddo
1382  enddo
1383 
1384  ! density at surface
1385  do j = jsb, jeb
1386  do i = isb, ieb
1387  rovcp_sfc = rtot_sfc(i,j) / ( cvtot_sfc(i,j) + rtot_sfc(i,j) )
1388  dens_sfc(1,i,j) = pres_sfc(1,i,j) / ( rtot_sfc(i,j) * temp_sfc(1,i,j) )
1389  pott_sfc(1,i,j) = temp_sfc(1,i,j) / ( p00/pres_sfc(1,i,j) )**rovcp_sfc
1390  enddo
1391  enddo
1392 
1393  ! make density at lowermost cell center
1394  do j = jsb, jeb
1395  do i = isb, ieb
1396  dz = real_cz(ks,i,j) - real_fz(ks-1,i,j)
1397 
1398  dens_s = 0.0_rp
1399  dens(ks,i,j) = dens_sfc(1,i,j) ! first guess
1400 
1401  converged = .false.
1402  do ite = 1, itelim
1403  if ( abs(dens(ks,i,j)-dens_s) <= criteria ) then
1404  converged = .true.
1405  exit
1406  endif
1407 
1408  dens_s = dens(ks,i,j)
1409 
1410  dhyd = + ( dens_sfc(1,i,j) * rtot_sfc(i,j) * temp_sfc(1,i,j) &
1411  - dens_s * rtot(i,j) * temp(ks,i,j) ) / dz & ! dp/dz
1412  - grav * 0.5_rp * ( dens_sfc(1,i,j) + dens_s ) ! rho*g
1413 
1414  dgrd = - rtot(i,j) * temp(ks,i,j) / dz &
1415  - 0.5_rp * grav
1416 
1417  dens(ks,i,j) = dens_s - dhyd/dgrd
1418 
1419  if( dens(ks,i,j)*0.0_rp /= 0.0_rp) exit
1420  enddo
1421 
1422  if ( .NOT. converged ) then
1423  write(*,*) 'xxx [buildrho bytemp 3D sfc] iteration not converged!', &
1424  i,j,dens(ks,i,j),ite,dens_s,dhyd,dgrd
1425  call prc_mpistop
1426  endif
1427  enddo
1428  enddo
1429 
1430  !--- from lowermost atmosphere to top of atmosphere
1431  call atmos_hydrostatic_buildrho_bytemp_atmos_3d( dens(:,:,:), & ! [INOUT]
1432  pott(:,:,:), & ! [OUT]
1433  pres(:,:,:), & ! [OUT]
1434  temp(:,:,:), & ! [IN]
1435  qv(:,:,:), & ! [IN]
1436  qc(:,:,:) ) ! [IN]
1437 
1438  return
1439  end subroutine atmos_hydrostatic_buildrho_bytemp_3d
1440 
1441  !-----------------------------------------------------------------------------
1443  subroutine atmos_hydrostatic_buildrho_bytemp_atmos_1d( &
1444  dens, &
1445  pott, &
1446  pres, &
1447  temp, &
1448  qv, &
1449  qc )
1450  use scale_process, only: &
1451  prc_mpistop
1452  implicit none
1453 
1454  real(RP), intent(inout) :: dens(ka)
1455  real(RP), intent(out) :: pott(ka)
1456  real(RP), intent(out) :: pres(ka)
1457  real(RP), intent(in) :: temp(ka)
1458  real(RP), intent(in) :: qv (ka)
1459  real(RP), intent(in) :: qc (ka)
1460 
1461  real(RP) :: Rtot (ka)
1462  real(RP) :: CVtot (ka)
1463 
1464  real(RP) :: RovCP
1465  real(RP) :: dens_s, dhyd, dgrd
1466  integer :: ite
1467  logical :: converged
1468 
1469  integer :: k
1470  !---------------------------------------------------------------------------
1471 
1472  do k = ks, ke
1473  rtot(k) = rdry * ( 1.0_rp - qv(k) - qc(k) ) &
1474  + rvap * qv(k)
1475  cvtot(k) = cvdry * ( 1.0_rp - qv(k) - qc(k) ) &
1476  + cv_qv * qv(k) &
1477  + cv_qc * qc(k)
1478  enddo
1479 
1480  do k = ks+1, ke
1481 
1482  dens_s = 0.0_rp
1483  dens(k) = dens(k-1) ! first guess
1484 
1485  converged = .false.
1486  do ite = 1, itelim
1487  if ( abs(dens(k)-dens_s) <= criteria ) then
1488  converged = .true.
1489  exit
1490  endif
1491 
1492  dens_s = dens(k)
1493 
1494  dhyd = + ( dens(k-1) * rtot(k-1) * temp(k-1) &
1495  - dens_s * rtot(k ) * temp(k ) ) / grid_fdz(k-1) & ! dp/dz
1496  - grav * 0.5_rp * ( dens(k-1) + dens_s ) ! rho*g
1497 
1498  dgrd = - rtot(k) * temp(k) / grid_fdz(k-1) &
1499  - 0.5_rp * grav
1500 
1501  dens(k) = dens_s - dhyd/dgrd
1502 
1503  if( dens(k)*0.0_rp /= 0.0_rp) exit
1504  enddo
1505 
1506  if ( .NOT. converged ) then
1507  write(*,*) 'xxx [buildrho bytemp 1D atmos] iteration not converged!', &
1508  k,dens(k),ite,dens_s,dhyd,dgrd
1509  call prc_mpistop
1510  endif
1511  enddo
1512 
1513  do k = ks, ke
1514  rovcp = rtot(k) / ( cvtot(k) + rtot(k) )
1515  pres(k) = dens(k) * rtot(k) * temp(k)
1516  pott(k) = temp(k) * ( p00 / pres(k) )**rovcp
1517  enddo
1518 
1519  return
1520  end subroutine atmos_hydrostatic_buildrho_bytemp_atmos_1d
1521 
1522  !-----------------------------------------------------------------------------
1524  subroutine atmos_hydrostatic_buildrho_bytemp_atmos_3d( &
1525  dens, &
1526  pott, &
1527  pres, &
1528  temp, &
1529  qv, &
1530  qc )
1531  use scale_process, only: &
1532  prc_mpistop
1533  implicit none
1534 
1535  real(RP), intent(inout) :: dens(ka,ia,ja)
1536  real(RP), intent(out) :: pott(ka,ia,ja)
1537  real(RP), intent(out) :: pres(ka,ia,ja)
1538  real(RP), intent(in) :: temp(ka,ia,ja)
1539  real(RP), intent(in) :: qv (ka,ia,ja)
1540  real(RP), intent(in) :: qc (ka,ia,ja)
1541 
1542  real(RP) :: Rtot (ka,ia,ja)
1543  real(RP) :: CVtot (ka,ia,ja)
1544 
1545  real(RP) :: RovCP
1546  real(RP) :: DZ
1547  real(RP) :: dens_s, dhyd, dgrd
1548  integer :: ite
1549  logical :: converged
1550 
1551  integer :: k, i, j
1552  !---------------------------------------------------------------------------
1553 
1554  do j = jsb, jeb
1555  do i = isb, ieb
1556  do k = ks, ke
1557  rtot(k,i,j) = rdry * ( 1.0_rp - qv(k,i,j) - qc(k,i,j) ) &
1558  + rvap * qv(k,i,j)
1559  cvtot(k,i,j) = cvdry * ( 1.0_rp - qv(k,i,j) - qc(k,i,j) ) &
1560  + cv_qv * qv(k,i,j) &
1561  + cv_qc * qc(k,i,j)
1562  enddo
1563  enddo
1564  enddo
1565 
1566  do j = jsb, jeb
1567  do i = isb, ieb
1568  do k = ks+1, ke
1569  dz = real_cz(k,i,j) - real_cz(k-1,i,j)
1570 
1571  dens_s = 0.0_rp
1572  dens(k,i,j) = dens(k-1,i,j) ! first guess
1573 
1574  converged = .false.
1575  do ite = 1, itelim
1576  if ( abs(dens(k,i,j)-dens_s) <= criteria ) then
1577  converged = .true.
1578  exit
1579  endif
1580 
1581  dens_s = dens(k,i,j)
1582 
1583  dhyd = + ( dens(k-1,i,j) * rtot(k-1,i,j) * temp(k-1,i,j) &
1584  - dens_s * rtot(k ,i,j) * temp(k ,i,j) ) / dz & ! dpdz
1585  - grav * 0.5_rp * ( dens(k-1,i,j) + dens_s ) ! rho*g
1586 
1587  dgrd = - rtot(k,i,j) * temp(k,i,j) / dz &
1588  - 0.5_rp * grav
1589 
1590  dens(k,i,j) = dens_s - dhyd/dgrd
1591 
1592  if( dens(k,i,j)*0.0_rp /= 0.0_rp) exit
1593  enddo
1594 
1595  if ( .NOT. converged ) then
1596  write(*,*) 'xxx [buildrho bytemp 3D atmos] iteration not converged!', &
1597  k,i,j,dens(k,i,j),ite,dens_s,dhyd,dgrd
1598  call prc_mpistop
1599  endif
1600  enddo
1601  enddo
1602  enddo
1603 
1604  do j = jsb, jeb
1605  do i = isb, ieb
1606  do k = ks, ke
1607  rovcp = rtot(k,i,j) / ( cvtot(k,i,j) + rtot(k,i,j) )
1608  pres(k,i,j) = dens(k,i,j) * rtot(k,i,j) * temp(k,i,j)
1609  pott(k,i,j) = temp(k,i,j) * ( p00 / pres(k,i,j) )**rovcp
1610  enddo
1611  enddo
1612  enddo
1613 
1614  return
1615  end subroutine atmos_hydrostatic_buildrho_bytemp_atmos_3d
1616 
1617  !-----------------------------------------------------------------------------
1619  subroutine atmos_hydrostatic_barometric_law_mslp_0d( &
1620  mslp, &
1621  pres, &
1622  temp, &
1623  dz )
1624  implicit none
1625 
1626  real(RP), intent(out) :: mslp
1627  real(RP), intent(in) :: pres
1628  real(RP), intent(in) :: temp
1629  real(RP), intent(in) :: dz
1630 
1631  ! work
1632  real(RP) :: TM
1633  !---------------------------------------------------------------------------
1634 
1635  tm = temp + laps * dz * 0.5_rp ! column-mean air temperature
1636 
1637  ! barometric law
1638  mslp = pres * exp( grav * dz / ( rdry * tm ) )
1639 
1640  return
1641  end subroutine atmos_hydrostatic_barometric_law_mslp_0d
1642 
1643  !-----------------------------------------------------------------------------
1645  subroutine atmos_hydrostatic_barometric_law_mslp_2d( &
1646  mslp, &
1647  pres, &
1648  temp, &
1649  dz )
1650  implicit none
1651 
1652  real(RP), intent(out) :: mslp(ia,ja)
1653  real(RP), intent(in) :: pres(ia,ja)
1654  real(RP), intent(in) :: temp(ia,ja)
1655  real(RP), intent(in) :: dz (ia,ja)
1656 
1657  ! work
1658  real(RP) :: TM
1659 
1660  integer :: i, j
1661  !---------------------------------------------------------------------------
1662 
1663  do j = jsb, jeb
1664  do i = isb, ieb
1665  tm = temp(i,j) + laps * dz(i,j) * 0.5_rp ! column-mean air temperature
1666 
1667  ! barometric law
1668  mslp(i,j) = pres(i,j) * exp( grav * dz(i,j) / ( rdry * tm ) )
1669  enddo
1670  enddo
1671 
1672  return
1673  end subroutine atmos_hydrostatic_barometric_law_mslp_2d
1674 
1675  !-----------------------------------------------------------------------------
1677  subroutine atmos_hydrostatic_barometric_law_pres_0d( &
1678  pres, &
1679  mslp, &
1680  temp, &
1681  dz )
1682  implicit none
1683 
1684  real(RP), intent(out) :: pres
1685  real(RP), intent(in) :: mslp
1686  real(RP), intent(in) :: temp
1687  real(RP), intent(in) :: dz
1688 
1689  ! work
1690  real(RP) :: TM
1691  !---------------------------------------------------------------------------
1692 
1693  tm = temp + laps * dz * 0.5_rp ! column-mean air temperature
1694 
1695  ! barometric law
1696  pres = mslp / exp( grav * dz / ( rdry * tm ) )
1697 
1698  return
1699  end subroutine atmos_hydrostatic_barometric_law_pres_0d
1700 
1701  !-----------------------------------------------------------------------------
1703  subroutine atmos_hydrostatic_barometric_law_pres_2d( &
1704  pres, &
1705  mslp, &
1706  temp, &
1707  dz )
1708  implicit none
1709 
1710  real(RP), intent(out) :: pres(ia,ja)
1711  real(RP), intent(in) :: mslp(ia,ja)
1712  real(RP), intent(in) :: temp(ia,ja)
1713  real(RP), intent(in) :: dz (ia,ja)
1714 
1715  ! work
1716  real(RP) :: TM
1717 
1718  integer :: i, j
1719  !---------------------------------------------------------------------------
1720 
1721  do j = jsb, jeb
1722  do i = isb, ieb
1723  tm = temp(i,j) + laps * dz(i,j) * 0.5_rp ! column-mean air temperature
1724 
1725  ! barometric law
1726  pres(i,j) = mslp(i,j) / exp( grav * dz(i,j) / ( rdry * tm ) )
1727  enddo
1728  enddo
1729 
1730  return
1731  end subroutine atmos_hydrostatic_barometric_law_pres_2d
1732 
1733 end module scale_atmos_hydrostatic
scale_atmos_hydrostatic::atmos_hydrostatic_buildrho_3d
subroutine atmos_hydrostatic_buildrho_3d(dens, temp, pres, pott, qv, qc, temp_sfc, pres_sfc, pott_sfc, qv_sfc, qc_sfc)
Build up density from surface (3D)
Definition: scale_atmos_sub_hydrostatic.F90:305
scale_atmos_hydrostatic::atmos_hydrostatic_buildrho_bytemp_atmos_1d
subroutine atmos_hydrostatic_buildrho_bytemp_atmos_1d(dens, pott, pres, temp, qv, qc)
Build up density from lowermost atmosphere (1D)
Definition: scale_atmos_sub_hydrostatic.F90:1450
scale_const::const_cvdry
real(rp), public const_cvdry
specific heat (dry air,constant volume) [J/kg/K]
Definition: scale_const.F90:59
scale_process::prc_mpistop
subroutine, public prc_mpistop
Abort MPI.
Definition: scale_process.F90:234
scale_grid_index::jeb
integer, public jeb
Definition: scale_grid_index.F90:74
scale_atmos_hydrostatic::atmos_hydrostatic_buildrho_atmos_rev_3d
subroutine, public atmos_hydrostatic_buildrho_atmos_rev_3d(dens, temp, pres, pott, qv, qc, dz, kref_in)
Build up density from lowermost atmosphere (3D)
Definition: scale_atmos_sub_hydrostatic.F90:1114
scale_stdio::io_l
logical, public io_l
output log or not? (this process)
Definition: scale_stdio.F90:59
scale_grid::grid_cz
real(rp), dimension(:), allocatable, public grid_cz
center coordinate [m]: z, local=global
Definition: scale_grid_cartesian.F90:56
scale_const::const_cl
real(rp), parameter, public const_cl
specific heat (liquid water) [J/kg/K]
Definition: scale_const.F90:68
scale_stdio
module STDIO
Definition: scale_stdio.F90:12
scale_grid_index::ke
integer, public ke
end point of inner domain: z, local
Definition: scale_grid_index.F90:59
scale_grid_real::real_fz
real(rp), dimension(:,:,:), allocatable, public real_fz
geopotential height [m] (cell face )
Definition: scale_grid_real.F90:41
scale_grid_real::real_cz
real(rp), dimension(:,:,:), allocatable, public real_cz
geopotential height [m] (cell center)
Definition: scale_grid_real.F90:40
scale_const::const_cvvap
real(rp), public const_cvvap
specific heat (water vapor, constant volume) [J/kg/K]
Definition: scale_const.F90:67
scale_const::const_laps
real(rp), public const_laps
lapse rate of ISA [K/m]
Definition: scale_const.F90:60
scale_const::const_rdry
real(rp), public const_rdry
specific gas constant (dry air) [J/kg/K]
Definition: scale_const.F90:57
scale_grid_index::ieb
integer, public ieb
Definition: scale_grid_index.F90:72
scale_grid_index
module grid index
Definition: scale_grid_index.F90:14
scale_grid_index::ia
integer, public ia
of x whole cells (local, with HALO)
Definition: scale_grid_index.F90:55
scale_comm::comm_horizontal_mean
subroutine, public comm_horizontal_mean(varmean, var)
calculate horizontal mean (global total with communication)
Definition: scale_comm.F90:516
scale_grid::grid_fdz
real(rp), dimension(:), allocatable, public grid_fdz
z-length of grid(k+1) to grid(k) [m]
Definition: scale_grid_cartesian.F90:69
scale_atmos_hydrostatic::atmos_hydrostatic_barometric_law_mslp_0d
subroutine atmos_hydrostatic_barometric_law_mslp_0d(mslp, pres, temp, dz)
Calculate mean sea-level pressure from barometric law (0D)
Definition: scale_atmos_sub_hydrostatic.F90:1624
scale_grid_real
module GRID (real space)
Definition: scale_grid_real.F90:14
scale_grid_index::ka
integer, public ka
of z whole cells (local, with HALO)
Definition: scale_grid_index.F90:54
scale_const::const_lapsdry
real(rp), public const_lapsdry
dry adiabatic lapse rate [K/m]
Definition: scale_const.F90:61
scale_const::const_pre00
real(rp), public const_pre00
pressure reference [Pa]
Definition: scale_const.F90:93
scale_comm
module COMMUNICATION
Definition: scale_comm.F90:23
scale_atmos_hydrostatic
module ATMOSPHERE / Hydrostatic barance
Definition: scale_atmos_sub_hydrostatic.F90:11
scale_const::const_grav
real(rp), public const_grav
standard acceleration of gravity [m/s2]
Definition: scale_const.F90:48
scale_process
module PROCESS
Definition: scale_process.F90:14
scale_atmos_hydrostatic::atmos_hydrostatic_buildrho_atmos_1d
subroutine atmos_hydrostatic_buildrho_atmos_1d(dens, temp, pres, pott, qv, qc)
Build up density from lowermost atmosphere (1D)
Definition: scale_atmos_sub_hydrostatic.F90:696
scale_const::const_rvap
real(rp), parameter, public const_rvap
specific gas constant (water vapor) [J/kg/K]
Definition: scale_const.F90:65
scale_const
module CONSTANT
Definition: scale_const.F90:14
scale_atmos_hydrostatic::atmos_hydrostatic_buildrho_1d
subroutine atmos_hydrostatic_buildrho_1d(dens, temp, pres, pott, qv, qc, temp_sfc, pres_sfc, pott_sfc, qv_sfc, qc_sfc)
Build up density from surface (1D)
Definition: scale_atmos_sub_hydrostatic.F90:180
scale_grid_index::ks
integer, public ks
start point of inner domain: z, local
Definition: scale_grid_index.F90:58
scale_atmos_hydrostatic::atmos_hydrostatic_buildrho_atmos_3d
subroutine atmos_hydrostatic_buildrho_atmos_3d(dens, temp, pres, pott, qv, qc, dz, kref_in)
Build up density from lowermost atmosphere (3D)
Definition: scale_atmos_sub_hydrostatic.F90:898
scale_grid
module GRID (cartesian)
Definition: scale_grid_cartesian.F90:17
scale_atmos_hydrostatic::atmos_hydrostatic_buildrho_bytemp_3d
subroutine atmos_hydrostatic_buildrho_bytemp_3d(dens, pott, pres, temp, qv, qc, pott_sfc, pres_sfc, temp_sfc, qv_sfc, qc_sfc)
Build up density from surface (3D)
Definition: scale_atmos_sub_hydrostatic.F90:1329
scale_prof
module profiler
Definition: scale_prof.F90:10
scale_const::const_eps
real(rp), public const_eps
small number
Definition: scale_const.F90:36
scale_atmos_hydrostatic::atmos_hydrostatic_buildrho_bytemp_1d
subroutine atmos_hydrostatic_buildrho_bytemp_1d(dens, pott, pres, temp, qv, qc, pott_sfc, pres_sfc, temp_sfc, qv_sfc, qc_sfc)
Build up density from surface (1D)
Definition: scale_atmos_sub_hydrostatic.F90:1224
scale_stdio::io_lnml
logical, public io_lnml
output log or not? (for namelist, this process)
Definition: scale_stdio.F90:60
scale_atmos_hydrostatic::atmos_hydrostatic_buildrho_real_3d
subroutine atmos_hydrostatic_buildrho_real_3d(dens, temp, pres, pott, qv, qc)
Build up density from surface (3D), not to reverse from TOA.
Definition: scale_atmos_sub_hydrostatic.F90:494
scale_precision
module PRECISION
Definition: scale_precision.F90:13
scale_atmos_hydrostatic::atmos_hydrostatic_buildrho_bytemp_atmos_3d
subroutine atmos_hydrostatic_buildrho_bytemp_atmos_3d(dens, pott, pres, temp, qv, qc)
Build up density from lowermost atmosphere (3D)
Definition: scale_atmos_sub_hydrostatic.F90:1531
scale_grid_index::isb
integer, public isb
Definition: scale_grid_index.F90:71
scale_stdio::io_fid_conf
integer, public io_fid_conf
Config file ID.
Definition: scale_stdio.F90:55
scale_stdio::io_fid_log
integer, public io_fid_log
Log file ID.
Definition: scale_stdio.F90:56
scale_grid_index::jsb
integer, public jsb
Definition: scale_grid_index.F90:73
scale_const::const_thermodyn_type
character(len=h_short), public const_thermodyn_type
internal energy type
Definition: scale_const.F90:101
scale_atmos_hydrostatic::atmos_hydrostatic_buildrho_atmos_rev_2d
subroutine, public atmos_hydrostatic_buildrho_atmos_rev_2d(dens_L1, temp_L1, pres_L1, pott_L1, qv_L1, qc_L1, dens_L2, pott_L2, qv_L2, qc_L2, dz, k)
Build up density (2D)
Definition: scale_atmos_sub_hydrostatic.F90:1009
scale_atmos_hydrostatic::atmos_hydrostatic_setup
subroutine, public atmos_hydrostatic_setup
Setup.
Definition: scale_atmos_sub_hydrostatic.F90:122
scale_atmos_hydrostatic::atmos_hydrostatic_buildrho_atmos_0d
subroutine, public atmos_hydrostatic_buildrho_atmos_0d(dens_L2, temp_L2, pres_L2, pott_L2, qv_L2, qc_L2, dens_L1, pott_L1, qv_L1, qc_L1, dz, k)
Build up density (0D)
Definition: scale_atmos_sub_hydrostatic.F90:607
scale_grid_index::ja
integer, public ja
of y whole cells (local, with HALO)
Definition: scale_grid_index.F90:56
Generated by   doxygen 1.8.13