Changeset d6a0977 in flexpart.git for src/com_mod.f90

Timestamp:

Dec 14, 2015, 3:10:04 PM (8 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, univie

Children:

f75967d

Parents:

88d8c3d

Message:

Updates to Henrik's wet depo scheme

File:

• src/com_mod.f90 (modified) (6 diffs)

src/com_mod.f90

@@ -130 +130 @@
   ! gdomainfill             .T., if domain-filling is global, .F. if not
 
-!hg
+!ZHG SEP 2015 wheather or not to read clouds from GRIB
   logical :: readclouds
 
@@ -347 +347 @@
   real :: tt(0:nxmax-1,0:nymax-1,nzmax,numwfmem)
   real :: qv(0:nxmax-1,0:nymax-1,nzmax,numwfmem)
-!hg adding cloud water 
-  real :: clwc(0:nxmax-1,0:nymax-1,nzmax,numwfmem)
-  real :: ciwc(0:nxmax-1,0:nymax-1,nzmax,numwfmem)
+!ZHG adding cloud water 
+  real :: clwc(0:nxmax-1,0:nymax-1,nzmax,numwfmem) !liquid   [kg/kg]
+  real :: ciwc(0:nxmax-1,0:nymax-1,nzmax,numwfmem) !ice      [kg/kg]
+  real :: clw(0:nxmax-1,0:nymax-1,nzmax,numwfmem)  !combined [m3/m3]
+
   real :: pv(0:nxmax-1,0:nymax-1,nzmax,numwfmem)
   real :: rho(0:nxmax-1,0:nymax-1,nzmax,numwfmem)
@@ -362 +364 @@
   integer(kind=1) :: clouds(0:nxmax-1,0:nymax-1,nzmax,numwfmem)
   integer :: cloudsh(0:nxmax-1,0:nymax-1,numwfmem)
-  ! integer :: icloudbot(0:nxmax-1,0:nymax-1,numwfmem)
-  ! integer :: icloudthck(0:nxmax-1,0:nymax-1,numwfmem)
+
+!ZHG Sep 2015  
+   real :: icloud_stats(0:nxmax-1,0:nymax-1,5,numwfmem)
+   real :: clh(0:nxmax-1,0:nymax-1,nuvzmax,numwfmem)
 
 
@@ -405 +409 @@
   real :: tropopause(0:nxmax-1,0:nymax-1,1,numwfmem)
   real :: oli(0:nxmax-1,0:nymax-1,1,numwfmem)
-! real :: diffk(0:nxmax-1,0:nymax-1,1,numwfmem) this is not in use?
-
+! real :: diffk(0:nxmax-1,0:nymax-1,1,numwfmem) ESO: this is not in use?
+! logical :: beneath_cloud=.true.
   ! ps                   surface pressure
   ! sd                   snow depth
@@ -667 +671 @@
 
 
+  !CGZ-lifetime
+  real, allocatable, dimension(:,:) ::checklifetime, species_lifetime
+  !CGZ-lifetime
+
   ! numpart                 actual number of particles in memory
   ! itra1 (maxpart) [s]     temporal positions of the particles
@@ -759 +767 @@
          & idt(nmpart),itramem(nmpart),itrasplit(nmpart),&
          & xtra1(nmpart),ytra1(nmpart),ztra1(nmpart),&
-         & xmass1(nmpart, maxspec))
+         & xmass1(nmpart, maxspec),&
+         & checklifetime(nmpart,maxspec), species_lifetime(maxspec,2))!CGZ-lifetime
+
 
     allocate(uap(nmpart),ucp(nmpart),uzp(nmpart),us(nmpart),&