Changeset 5f9d14a in flexpart.git for src/readwind.f90

Timestamp:

Apr 8, 2015, 2:23:27 PM (9 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:

1585284

Parents:

cd85138

Message:

Updated wet depo scheme

File:

• src/readwind.f90 (modified) (5 diffs)

src/readwind.f90

@@ -77 +77 @@
 !HSO  end
 
-  real(kind=4) :: uuh(0:nxmax-1,0:nymax-1,nuvzmax)
-  real(kind=4) :: vvh(0:nxmax-1,0:nymax-1,nuvzmax)
-  real(kind=4) :: wwh(0:nxmax-1,0:nymax-1,nwzmax)
+  real(sp) :: uuh(0:nxmax-1,0:nymax-1,nuvzmax)
+  real(sp) :: vvh(0:nxmax-1,0:nymax-1,nuvzmax)
+  real(sp) :: wwh(0:nxmax-1,0:nymax-1,nwzmax)
   integer :: indj,i,j,k,n,levdiff2,ifield,iumax,iwmax
 
@@ -91 +91 @@
 
   integer :: isec1(56),isec2(22+nxmax+nymax)
-  real(kind=4) :: zsec4(jpunp)
-  real(kind=4) :: xaux,yaux,xaux0,yaux0
-  real(kind=8) :: xauxin,yauxin
-  real,parameter :: eps=1.e-4
-  real(kind=4) :: nsss(0:nxmax-1,0:nymax-1),ewss(0:nxmax-1,0:nymax-1)
-  real :: plev1,pmean,tv,fu,hlev1,ff10m,fflev1,conversion_factor
-
-  logical :: hflswitch,strswitch, readcloud
+  real(sp) :: zsec4(jpunp)
+  real(sp) :: xaux,yaux
+  real(dp) :: xauxin,yauxin
+  real(sp),parameter :: eps=1.e-4
+  real(sp) :: nsss(0:nxmax-1,0:nymax-1),ewss(0:nxmax-1,0:nymax-1)
+  real(sp) :: plev1,pmean,tv,fu,hlev1,ff10m,fflev1,conversion_factor
+
+  logical :: hflswitch,strswitch!,readcloud
 
 !HSO  grib api error messages
@@ -107 +107 @@
   strswitch=.false.
 !hg
-  readcloud=.false.
+!  readcloud=.false.
 !hg end
   levdiff2=nlev_ec-nwz+1
@@ -362 +362 @@
         !write(*,*) 'found water!'
       endif
-
       if(isec1(6).eq.247) then  !! CIWC  Cloud ice water content
         ciwch(i,j,nlev_ec-k+2,n)=zsec4(nxfield*(ny-j-1)+i+1)
@@ -379 +378 @@
 
 50 call grib_close_file(ifile)
-
 
 !error message if no fields found with correct first longitude in it