Changeset 34f1452 in flexpart.git for src/plumetraj.f90
- Timestamp:
- Sep 10, 2018, 10:27:33 AM (6 years ago)
- Branches:
- univie
- Children:
- d7935de
- Parents:
- 2870133
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
src/plumetraj.f90
r6a678e3 r34f1452 32 32 ! * 33 33 ! 24 January 2002 * 34 ! * 35 ! 2018-09-10, Petra Seibert: increase output width for time in sec * 36 ! small improvements in code layout, name big loops * 34 37 ! * 35 38 ! Variables: * … … 77 80 78 81 79 ! Loop aboutall release points80 !***************************** *81 82 do j=1,numpoint82 ! Loop over all release points 83 !***************************** 84 85 point_loop: do j=1,numpoint 83 86 if (abs(ireleasestart(j)-itime).gt.lage(nageclass)) goto 10 84 87 topocenter=0. … … 93 96 94 97 n=0 95 do i=1,numpart98 particle_loop: do i=1,numpart 96 99 if (itra1(i).ne.itime) goto 20 97 100 if (npoint(i).ne.j) goto 20 … … 125 128 126 129 topo=p1*oro(ix ,jy) & 127 128 129 130 + p2*oro(ixp,jy) & 131 + p3*oro(ix ,jyp) & 132 + p4*oro(ixp,jyp) 130 133 topocenter=topocenter+topo 131 134 … … 151 154 indexh=memind(m) 152 155 pv1(m)=p1*pv(ix ,jy ,ind,indexh) & 153 +p2*pv(ixp,jy ,ind,indexh) &154 +p3*pv(ix ,jyp,ind,indexh) &155 +p4*pv(ixp,jyp,ind,indexh)156 +p2*pv(ixp,jy ,ind,indexh) & 157 +p3*pv(ix ,jyp,ind,indexh) & 158 +p4*pv(ixp,jyp,ind,indexh) 156 159 end do 157 160 pvprof(ind-indz+1)=(pv1(1)*dt2+pv1(2)*dt1)*dtt … … 173 176 174 177 tr(m)=p1*tropopause(ix ,jy ,1,indexh) & 175 176 177 178 + p2*tropopause(ixp,jy ,1,indexh) & 179 + p3*tropopause(ix ,jyp,1,indexh) & 180 + p4*tropopause(ixp,jyp,1,indexh) 178 181 179 182 hm(m)=p1*hmix(ix ,jy ,1,indexh) & 180 181 182 183 + p2*hmix(ixp,jy ,1,indexh) & 184 + p3*hmix(ix ,jyp,1,indexh) & 185 + p4*hmix(ixp,jyp,1,indexh) 183 186 end do 184 187 … … 191 194 hmixcenter=hmixcenter+hmixi 192 195 193 194 196 20 continue 195 end do 197 end do particle_loop 196 198 197 199 … … 234 236 !****************************************** 235 237 236 write(unitouttraj,'(i5,i8,2f9.4,4f8.1,f8.2,4f8.1,3f6.1,& 237 &5(2f8.3,f7.0,f6.1,f8.1))')& 238 &j,itime-(ireleasestart(j)+ireleaseend(j))/2, & 238 write (unitouttraj,901) j,itime-(ireleasestart(j)+ireleaseend(j))/2, & 239 239 xcenter,ycenter,zcenter,topocenter,hmixcenter,tropocenter, & 240 pvcenter,rmsdist,rms,zrmsdist,zrms,hmixfract,pvfract, & 241 tropofract, & 242 (xclust(k),yclust(k),zclust(k),fclust(k),rmsclust(k), & 243 k=1,ncluster) 240 pvcenter,rmsdist,rms,zrmsdist,zrms,hmixfract,pvfract,tropofract, & 241 (xclust(k),yclust(k),zclust(k),fclust(k),rmsclust(k), k=1,ncluster) 242 901 format(i5,i10,2f9.4,4f8.1,f8.2,4f8.1,3f6.1,5(2f8.3,f7.0,f6.1,f8.1)) 244 243 endif 245 244 246 245 247 246 10 continue 248 end do 247 end do point_loop 249 248 250 249
Note: See TracChangeset
for help on using the changeset viewer.