Changeset 05cf28d in flexpart.git for src/readpartpositions.f90

Timestamp:

Apr 22, 2016, 9:35:28 AM (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:

356dcb9

Parents:

61df8d9

Message:

Updated checks and warning messages for wet deposition parameters given in SPECIES files

File:

• src/readpartpositions.f90 (modified) (1 diff)

src/readpartpositions.f90

@@ -105 +105 @@
   open(unitpartin,file=path(2)(1:length(2))//'partposit_end', &
        form='unformatted',err=998)
+  
 
+100 read(unitpartin,end=99) itimein
+  i=0
+200 i=i+1
+  read(unitpartin) npoint(i),xlonin,ylatin,ztra1(i),itramem(i), &
+       topo,pvi,qvi,rhoi,hmixi,tri,tti,(xmass1(i,j),j=1,nspec)
+  
+  if (xlonin.eq.-9999.9) goto 100
+  xtra1(i)=(xlonin-xlon0)/dx
+  ytra1(i)=(ylatin-ylat0)/dy
+  numparticlecount=max(numparticlecount,npoint(i))
+  goto 200
 
-100   read(unitpartin,end=99) itimein
-    i=0
-200   i=i+1
-    read(unitpartin) npoint(i),xlonin,ylatin,ztra1(i),itramem(i), &
-         topo,pvi,qvi,rhoi,hmixi,tri,tti,(xmass1(i,j),j=1,nspec)
-
-    if (xlonin.eq.-9999.9) goto 100
-    xtra1(i)=(xlonin-xlon0)/dx
-    ytra1(i)=(ylatin-ylat0)/dy
-    numparticlecount=max(numparticlecount,npoint(i))
-    goto 200
-
-99   numpart=i-1
+99 numpart=i-1
 
   close(unitpartin)