Changes in src/verttransform.f90 [6481010:341f4b7] in flexpart.git

File:

• src/verttransform.f90 (modified) (7 diffs)

src/verttransform.f90

@@ -91 +91 @@
 
   logical :: init = .true.
-  logical :: init_w = .false.
-  logical :: init_r = .false.
-
 
   !ZHG SEP 2014 tests  
@@ -106 +103 @@
   ! CHARACTER(LEN=3)  :: aspec
   ! integer :: virr=0
-  !real :: tot_cloud_h
-  !real :: dbg_height(nzmax) 
+  real :: tot_cloud_h
+  real :: dbg_height(nzmax) 
 !ZHG
 
@@ -126 +123 @@
   if (init) then
 
-
-    if (init_r) then
-
-        open(333,file='heights.txt', &
-          form='formatted')
-        do kz=1,nuvz
-            read(333,*) height(kz)
-        end do
-        close(333)
-        write(*,*) 'height read'
-    else
-
-
 ! Search for a point with high surface pressure (i.e. not above significant topography)
 ! Then, use this point to construct a reference z profile, to be used at all times
@@ -177 +161 @@
     end do
 
-    if (init_w) then
-        open(333,file='heights.txt', &
-          form='formatted')
-        do kz=1,nuvz
-              write(333,*) height(kz)
-        end do
-        close(333)
-    endif
-
-    endif ! init
 
 ! Determine highest levels that can be within PBL
@@ -204 +178 @@
     init=.false.
 
-!    dbg_height = height
+    dbg_height = height
 
   endif
@@ -624 +598 @@
         convp=convprec(ix,jy,1,n)
         prec=lsp+convp
-!        tot_cloud_h=0
+        tot_cloud_h=0
 ! Find clouds in the vertical
         do kz=1, nz-1 !go from top to bottom
@@ -630 +604 @@
 ! assuming rho is in kg/m3 and hz in m gives: kg/kg * kg/m3 *m3/kg /m = m2/m3 
             clw(ix,jy,kz,n)=(clwc(ix,jy,kz,n)*rho(ix,jy,kz,n))*(height(kz+1)-height(kz))
-!            tot_cloud_h=tot_cloud_h+(height(kz+1)-height(kz)) 
+            tot_cloud_h=tot_cloud_h+(height(kz+1)-height(kz)) 
             
 !            icloud_stats(ix,jy,4,n)= icloud_stats(ix,jy,4,n)+clw(ix,jy,kz,n)          ! Column cloud water [m3/m3]