Changeset 37 for branches/petra/src/com_mod.f90

Timestamp:

Apr 22, 2015, 3:42:35 PM (9 years ago)

Author:

pesei

Message:

Wet dep quick fix and other small changes. Wet depo quick fix not final yet.

File:

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

branches/petra/src/com_mod.f90

@@ -32 +32 @@
 !        Modifications: 15 August 2013 IP, 
 !        2/2015 PS, add incremental deposition switch
+!        2/2015 PS, change variables for cloud diagnostics
 !                                                                              *
 !*******************************************************************************
@@ -365 +366 @@
   real :: qvh(0:nxmax-1,0:nymax-1,nuvzmax,2)
   real :: pplev(0:nxmax-1,0:nymax-1,nuvzmax,2)
-  integer(kind=1) :: clouds(0:nxmax-1,0:nymax-1,nzmax,2)
-  integer :: cloudsh(0:nxmax-1,0:nymax-1,2)
-
-
   ! uu,vv,ww [m/2]       wind components in x,y and z direction
   ! uupol,vvpol [m/s]    wind components in polar stereographic projection
@@ -377 +374 @@
   ! drhodz [kg/m2]       vertical air density gradient
   ! tth,qvh              tth,qvh on original eta levels
-  ! clouds:   no cloud, no precipitation   0
-  !      cloud, no precipitation      1
-  !      rainout  conv/lsp dominated  2/3
-  !      washout  conv/lsp dominated  4/5
 
   ! pplev for the GFS version
@@ -406 +399 @@
   real :: oli(0:nxmax-1,0:nymax-1,1,2)
   real :: diffk(0:nxmax-1,0:nymax-1,1,2)
+  integer, dimension (0:nxmax-1,0:nymax-1,2) :: icloudbot, icloudthck
 
   ! ps                   surface pressure
@@ -426 +420 @@
   ! oli [m]              inverse Obukhov length (1/L)
   ! diffk [m2/s]         diffusion coefficient at reference height
+  ! icloudbot (m)        cloud bottom height
+  ! icloudthck (m)       cloud thickness   
 
 
@@ -484 +480 @@
   real :: qvn(0:nxmaxn-1,0:nymaxn-1,nzmax,2,maxnests)
   real :: pvn(0:nxmaxn-1,0:nymaxn-1,nzmax,2,maxnests)
-  integer(kind=1) :: cloudsn(0:nxmaxn-1,0:nymaxn-1,0:nzmax,2,maxnests)
-  integer :: cloudsnh(0:nxmaxn-1,0:nymaxn-1,2,maxnests)
   real :: rhon(0:nxmaxn-1,0:nymaxn-1,nzmax,2,maxnests)
   real :: drhodzn(0:nxmaxn-1,0:nymaxn-1,nzmax,2,maxnests)
@@ -514 +508 @@
   real :: diffkn(0:nxmaxn-1,0:nymaxn-1,1,2,maxnests)
   real :: vdepn(0:nxmaxn-1,0:nymaxn-1,maxspec,2,maxnests)
+  integer, dimension (0:nxmax-1,0:nymax-1,2,maxnests) :: icloudbotn, icloudthckn