Changeset 3c1da52 in flexpart.git for src/getfields.f90
- Timestamp:
- Jul 15, 2015, 9:45:31 AM (9 years ago)
- Branches:
- svn-petra
- Children:
- c07b09e
- Parents:
- 31674b5
- git-author:
- Ignacio Pisso <Ignacio.Pisso@…> (07/15/15 09:23:56)
- git-committer:
- Ignacio Pisso <Ignacio.Pisso@…> (07/15/15 09:45:31)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
src/getfields.f90
re200b7a r3c1da52 41 41 ! Variables tth,qvh,tthn,qvhn (on eta coordinates) in common block. 42 42 ! Function of nstop extended. 43 ! 44 ! 3/2015, PS: put in verbosity output 45 ! 43 46 !***************************************************************************** 44 47 ! * … … 83 86 integer :: indmin = 1 84 87 88 call verboutput(verbosity,1,' getfields') 85 89 86 90 ! Check, if wind fields are available for the current time step … … 125 129 do indj=indmin,numbwf-1 126 130 if (ldirect*wftime(indj+1).gt.ldirect*itime) then 131 call verboutput(verbosity,1,' call readwind') 127 132 call readwind(indj+1,memind(2),uuh,vvh,wwh) 133 call verboutput(verbosity,1,' call readwind_nests') 128 134 call readwind_nests(indj+1,memind(2),uuhn,vvhn,wwhn) 135 call verboutput(verbosity,1,' call calcpar') 129 136 call calcpar(memind(2),uuh,vvh,pvh) 137 call verboutput(verbosity,1,' call calcpar_nests') 130 138 call calcpar_nests(memind(2),uuhn,vvhn,pvhn) 139 call verboutput(verbosity,1,' call verttransform') 131 140 call verttransform(memind(2),uuh,vvh,wwh,pvh) 141 call verboutput(verbosity,1,' call verttransform_nests') 132 142 call verttransform_nests(memind(2),uuhn,vvhn,wwhn,pvhn) 133 143 memtime(2)=wftime(indj+1) … … 148 158 (ldirect*wftime(indj+1).gt.ldirect*itime)) then 149 159 memind(1)=1 160 call verboutput(verbosity,1,' call readwind') 150 161 call readwind(indj,memind(1),uuh,vvh,wwh) 162 call verboutput(verbosity,1,' call readwind_nests') 151 163 call readwind_nests(indj,memind(1),uuhn,vvhn,wwhn) 164 call verboutput(verbosity,1,' call calcpar') 152 165 call calcpar(memind(1),uuh,vvh,pvh) 166 call verboutput(verbosity,1,' call calcpar_nests') 153 167 call calcpar_nests(memind(1),uuhn,vvhn,pvhn) 168 call verboutput(verbosity,1,' call verttransform') 154 169 call verttransform(memind(1),uuh,vvh,wwh,pvh) 170 call verboutput(verbosity,1,' call verttransform_nests') 155 171 call verttransform_nests(memind(1),uuhn,vvhn,wwhn,pvhn) 156 172 memtime(1)=wftime(indj) 157 173 memind(2)=2 174 call verboutput(verbosity,1,' call readwind') 158 175 call readwind(indj+1,memind(2),uuh,vvh,wwh) 176 call verboutput(verbosity,1,' call readwind_nests') 159 177 call readwind_nests(indj+1,memind(2),uuhn,vvhn,wwhn) 178 call verboutput(verbosity,1,' call calcpar') 160 179 call calcpar(memind(2),uuh,vvh,pvh) 180 call verboutput(verbosity,1,' call calcpar_nests') 161 181 call calcpar_nests(memind(2),uuhn,vvhn,pvhn) 182 call verboutput(verbosity,1,' call verttransform') 162 183 call verttransform(memind(2),uuh,vvh,wwh,pvh) 184 call verboutput(verbosity,1,' call verttransform_nests') 163 185 call verttransform_nests(memind(2),uuhn,vvhn,wwhn,pvhn) 164 186 memtime(2)=wftime(indj+1) … … 175 197 if (lwindinterv.gt.idiffmax) nstop=3 176 198 199 call verboutput(verbosity,1,' leaving getfields') 200 177 201 end subroutine getfields
Note: See TracChangeset
for help on using the changeset viewer.