Changeset db712a8 in flexpart.git for src/calcpar_nests.f90

Timestamp:

Mar 3, 2016, 12:34:56 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:

38b7917

Parents:

b0434e1

Message:

Completed handling of nested wind fields with cloud water (for wet deposition).

File:

• src/calcpar_nests.f90 (modified) (4 diffs)

src/calcpar_nests.f90

@@ -103 +103 @@
       ustarn(ix,jy,1,n,l)=scalev(psn(ix,jy,1,n,l),tt2n(ix,jy,1,n,l), &
            td2n(ix,jy,1,n,l),surfstrn(ix,jy,1,n,l))
+      if (ustarn(ix,jy,1,n,l).le.1.e-8) ustarn(ix,jy,1,n,l)=1.e-8
 
   ! 2) Calculation of inverse Obukhov length scale
@@ -135 +136 @@
         subsceff=min(excessoron(ix,jy,l),hmixplus)
       else
-        subsceff=0
+        subsceff=0.0
       endif
   !
@@ -149 +150 @@
 
       if (DRYDEP) then
-        z0(4)=0.016*ustarn(ix,jy,1,n,l)*ustarn(ix,jy,1,n,l)/ga
-        z0(9)=0.016*ustarn(ix,jy,1,n,l)*ustarn(ix,jy,1,n,l)/ga
+        ! z0(4)=0.016*ustarn(ix,jy,1,n,l)*ustarn(ix,jy,1,n,l)/ga
+        ! z0(9)=0.016*ustarn(ix,jy,1,n,l)*ustarn(ix,jy,1,n,l)/ga
+        z0(7)=0.016*ustarn(ix,jy,1,n,l)*ustarn(ix,jy,1,n,l)/ga
 
   ! Calculate relative humidity at surface
@@ -228 +230 @@
   end do
 
-
-    call calcpv_nests(l,n,uuhn,vvhn,pvhn)
+  ! Calculation of potential vorticity on 3-d grid
+  !***********************************************
+
+  call calcpv_nests(l,n,uuhn,vvhn,pvhn)
 
   end do