Changeset 2eefa58 in flexpart.git for src/getfields.f90

Timestamp:

May 27, 2019, 3:28:44 PM (5 years ago)

Author:

Espen Sollum ATMOS <eso@…>

Branches:

master, 10.4.1_pesei, GFS_025, bugfixes+enhancements, dev, release-10, release-10.4.1, scaling-bug

Children:

6741557

Parents:

f963113

Message:

Added Ronas changes for inversion output

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
@@ -153 +157 @@
 40  indmin=indj
 
-   if (WETBKDEP) then
-        call writeprecip(itime,memind(1))
-   endif
+    if (WETBKDEP) then
+      call writeprecip(itime,memind(1))
+    endif
 
   else
@@ -204 +208 @@
 60  indmin=indj
 
-   if (WETBKDEP) then
-        call writeprecip(itime,memind(1))
-   endif
-
-  endif
+    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
 
   lwindinterv=abs(memtime(2)-memtime(1))