Changeset 8b3d324 in flexpart.git
- Timestamp:
- Jun 25, 2018, 11:28:52 AM (6 years ago)
- Branches:
- univie
- Children:
- f251e57
- Parents:
- 1a8fbee
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
src/verttransform_ecmwf.f90
r1a8fbee r8b3d324 21 21 22 22 subroutine verttransform_ecmwf(n,uuh,vvh,wwh,pvh) 23 ! i i i i i23 ! i i i i i 24 24 !***************************************************************************** 25 25 ! * … … 168 168 3 continue 169 169 170 print*,'oldheights at' ,ixref,jyref,ps(ixref,jyref,1,n)170 ! print*,'oldheights at' ,ixref,jyref,ps(ixref,jyref,1,n) 171 171 else ! new reference grid point 172 172 ! PS: the old version fails if the pressure is <=1000 hPa in the whole … … 186 186 psmean = sum( ps(:,:,1,n), ps(:,:,1,n) > psmean - psstd ) / & 187 187 count(ps(:,:,1,n) > psmean - psstd) 188 print*,'height: std, new psmean',psstd,psmean188 ! print*,'height: std, new psmean',psstd,psmean 189 189 ixyref = minloc( abs( ps(:,:,1,n) - psmean ) ) 190 190 ixref = ixyref(1) 191 191 jyref = ixyref(2) 192 print*,'newheights at' ,ixref,jyref,ps(ixref,jyref,1,n)192 ! print*,'newheights at' ,ixref,jyref,ps(ixref,jyref,1,n) 193 193 endif 194 194
Note: See TracChangeset
for help on using the changeset viewer.