Changeset 341f4b7 in flexpart.git for src/com_mod.f90

Timestamp:

Apr 20, 2016, 1:39:51 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:

61df8d9

Parents:

32b49c3

Message:

SPECIES files converted to namelist format. Fixed-format SPCIES files renamed with .oldformat extension. Added 2 wet depo parameters to old-format SPECIES. Renamed internal variables and parameters used for wet deposition. incloud_ratio increased to 2.2

File:

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

src/com_mod.f90

@@ -160 +160 @@
   !real xmass(maxpoint,maxspec)
   real :: decay(maxspec)
-  real :: weta(maxspec),wetb(maxspec)
+  real :: weta_gas(maxspec),wetb_gas(maxspec)
+  real :: crain_aero(maxspec),csnow_aero(maxspec)
 ! NIK: 31.01.2013- parameters for in-cloud scavening
-  real :: weta_in(maxspec), wetb_in(maxspec)
+  real :: ccn_aero(maxspec),in_aero(maxspec)
   real :: reldiff(maxspec),henry(maxspec),f0(maxspec)
   real :: density(maxspec),dquer(maxspec),dsigma(maxspec)
@@ -197 +198 @@
 
   ! WET DEPOSITION
-  ! weta, wetb              parameters for determining below-cloud wet scavenging coefficients
-  ! weta_in, wetb_in       parameters for determining in-cloud wet scavenging coefficients
-  ! wetc_in, wetd_in       parameters for determining in-cloud wet scavenging coefficients
+  ! weta_gas, wetb_gas     parameters for below-cloud wet scavenging coefficients (gasses)
+  ! crain_aero, in_aero    parameters for in-cloud wet scavenging coefficients (aerosols)
+  ! wetc_in, wetd_in       parameters for in-cloud wet scavenging coefficients (aerosols)
 
   ! GAS DEPOSITION
@@ -216 +217 @@
   ! dquer [m]               mean diameter of particles
   ! dsigma                  dsigma=10 or dsigma=0.1 means that 68% of the
-  !                    mass are between 0.1*dquer and 10*dquer
+  !                         mass are between 0.1*dquer and 10*dquer
 
   ! fract                   mass fraction of each diameter interval
@@ -371 +372 @@
 
 !ZHG Sep 2015  
-   real :: icloud_stats(0:nxmax-1,0:nymax-1,5,numwfmem)
-   real :: clw4(0:nxmax-1,0:nymax-1,numwfmem) ! eso: =icloud_stats(:,:,4,:)
+!   real :: icloud_stats(0:nxmax-1,0:nymax-1,5,numwfmem)
+   real :: ctwc(0:nxmax-1,0:nymax-1,numwfmem) ! eso: =icloud_stats(:,:,4,:)
 
 
@@ -392 +393 @@
 
   ! pplev for the GFS version
+  ! ctwc                 total cloud water content
 
   ! 2d fields
@@ -489 +491 @@
   real,allocatable,dimension(:,:,:,:,:) :: uun, vvn, wwn, ttn, qvn, pvn,&
        & rhon, drhodzn, tthn, qvhn, clwcn, ciwcn, clwn, clwchn, ciwchn
-  real,allocatable,dimension(:,:,:,:) :: clw4n
+  real,allocatable,dimension(:,:,:,:) :: ctwcn
   integer,allocatable,dimension(:,:,:,:) :: cloudshn
   integer(kind=1),allocatable,dimension(:,:,:,:,:) :: cloudsn
@@ -802 +804 @@
     allocate(clwchn(0:nxmaxn-1,0:nymaxn-1,nuvzmax,numwfmem,numbnests))
     allocate(ciwchn(0:nxmaxn-1,0:nymaxn-1,nuvzmax,numwfmem,numbnests))
-    allocate(clw4n(0:nxmax-1,0:nymax-1,numwfmem,numbnests))
-
-!    clw4n(:,:,:,:)=0.
+    allocate(ctwcn(0:nxmax-1,0:nymax-1,numwfmem,numbnests))
+
     clwcn(:,:,:,:,:)=0.
     ciwcn(:,:,:,:,:)=0.
     clwchn(:,:,:,:,:)=0.
     ciwchn(:,:,:,:,:)=0.
-    ! clwn(:,:,:,:,:)=0.
-
-    ! cloudsn(:,:,:,:,:)=0
-    ! cloudshn(:,:,:,:)=0
     
   end subroutine com_mod_allocate_nests