Changeset de4c5e9 in flexpart.git for src/verttransform_ecmwf.f90
- Timestamp:
- Oct 5, 2017, 9:09:23 AM (7 years ago)
- Branches:
- master, 10.4.1_pesei, GFS_025, bugfixes+enhancements, dev, release-10, release-10.4.1, scaling-bug, univie
- Children:
- 7233c03
- Parents:
- 08a38b5
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
src/verttransform_ecmwf.f90
r6ecb30a rde4c5e9 53 53 ! Marian Harustak, 12.5.2017 54 54 ! - Renamed from verttransform to verttransform_ecmwf 55 ! 55 !***************************************************************************** 56 ! Date: 2017-05-30 modification of a bug in ew. Don Morton (CTBTO project) * 56 57 !***************************************************************************** 57 58 ! * … … 96 97 97 98 logical :: init = .true. 99 logical :: init_w = .false. 100 logical :: init_r = .true. 101 98 102 99 103 !ZHG SEP 2014 tests … … 108 112 ! CHARACTER(LEN=3) :: aspec 109 113 ! integer :: virr=0 110 real :: tot_cloud_h111 real :: dbg_height(nzmax)114 !real :: tot_cloud_h 115 !real :: dbg_height(nzmax) 112 116 !ZHG 113 117 … … 128 132 if (init) then 129 133 134 135 if (init_r) then 136 137 open(333,file='heights.txt', & 138 form='formatted') 139 do kz=1,nuvz 140 read(333,*) height(kz) 141 end do 142 close(333) 143 write(*,*) 'height read' 144 else 145 146 130 147 ! Search for a point with high surface pressure (i.e. not above significant topography) 131 148 ! Then, use this point to construct a reference z profile, to be used at all times … … 166 183 end do 167 184 185 if (init_w) then 186 open(333,file='heights.txt', & 187 form='formatted') 188 do kz=1,nuvz 189 write(333,*) height(kz) 190 end do 191 close(333) 192 endif 193 194 endif ! init 168 195 169 196 ! Determine highest levels that can be within PBL … … 183 210 init=.false. 184 211 185 dbg_height = height212 ! dbg_height = height 186 213 187 214 endif … … 603 630 convp=convprec(ix,jy,1,n) 604 631 prec=lsp+convp 605 tot_cloud_h=0632 ! tot_cloud_h=0 606 633 ! Find clouds in the vertical 607 634 do kz=1, nz-1 !go from top to bottom … … 609 636 ! assuming rho is in kg/m3 and hz in m gives: kg/kg * kg/m3 *m3/kg /m = m2/m3 610 637 clw(ix,jy,kz,n)=(clwc(ix,jy,kz,n)*rho(ix,jy,kz,n))*(height(kz+1)-height(kz)) 611 tot_cloud_h=tot_cloud_h+(height(kz+1)-height(kz))638 ! tot_cloud_h=tot_cloud_h+(height(kz+1)-height(kz)) 612 639 613 640 ! icloud_stats(ix,jy,4,n)= icloud_stats(ix,jy,4,n)+clw(ix,jy,kz,n) ! Column cloud water [m3/m3]
Note: See TracChangeset
for help on using the changeset viewer.