Changes in src/getfields.f90 [2eefa58:6ecb30a] in flexpart.git

File:

• src/getfields.f90 (modified) (3 diffs)

src/getfields.f90

@@ -87 +87 @@
   real :: pvhn(0:nxmaxn-1,0:nymaxn-1,nuvzmax,maxnests)
   real :: wwhn(0:nxmaxn-1,0:nymaxn-1,nwzmax,maxnests)
-! RLT added partial pressure water vapor 
-  real :: pwater(0:nxmax-1,0:nymax-1,nzmax,numwfmem)
-  integer :: kz, ix
-  character(len=100) :: rowfmt
 
   integer :: indmin = 1
@@ -157 +153 @@
 40  indmin=indj
 
-    if (WETBKDEP) then
-      call writeprecip(itime,memind(1))
-    endif
+   if (WETBKDEP) then
+        call writeprecip(itime,memind(1))
+   endif
 
   else
@@ -208 +204 @@
 60  indmin=indj
 
-    if (WETBKDEP) then
-      call writeprecip(itime,memind(1))
-    endif
-
-  end if
-
-  ! RLT calculate dry air density
-  pwater=qv*prs/((r_air/r_water)*(1.-qv)+qv)
-  rho_dry=(prs-pwater)/(r_air*tt)
-
-  ! test density
-!  write(rowfmt,'(A,I6,A)') '(',nymax,'(E11.4,1X))'
-!  if(itime.eq.0) then
-!    open(500,file=path(2)(1:length(2))//'rho_dry.txt',status='replace',action='write')
-!    do kz=1,nzmax
-!      do ix=1,nxmax
-!        write(500,fmt=rowfmt) rho_dry(ix,:,kz,1)
-!      end do
-!    end do
-!    close(500)
-!    open(500,file=path(2)(1:length(2))//'rho.txt',status='replace',action='write')
-!    do kz=1,nzmax
-!      do ix=1,nxmax
-!        write(500,fmt=rowfmt) rho(ix,:,kz,1)
-!      end do
-!    end do
-!    close(500)
-!  endif
+   if (WETBKDEP) then
+        call writeprecip(itime,memind(1))
+   endif
+
+  endif
 
   lwindinterv=abs(memtime(2)-memtime(1))