Commit 8b3d324e authored by pesei's avatar pesei
Browse files

comment out print statements

comment out print statements for the reference height, adjust the i/o
comment line
parent 1a8fbee0
...@@ -20,7 +20,7 @@ ...@@ -20,7 +20,7 @@
! ********************************************************************** ! **********************************************************************
subroutine verttransform_ecmwf(n,uuh,vvh,wwh,pvh) subroutine verttransform_ecmwf(n,uuh,vvh,wwh,pvh)
! i i i i i ! i i i i i
!***************************************************************************** !*****************************************************************************
! * ! *
! This subroutine transforms temperature, dew point temperature and * ! This subroutine transforms temperature, dew point temperature and *
...@@ -167,7 +167,7 @@ subroutine verttransform_ecmwf(n,uuh,vvh,wwh,pvh) ...@@ -167,7 +167,7 @@ subroutine verttransform_ecmwf(n,uuh,vvh,wwh,pvh)
end do end do
3 continue 3 continue
print*,'oldheights at' ,ixref,jyref,ps(ixref,jyref,1,n) ! print*,'oldheights at' ,ixref,jyref,ps(ixref,jyref,1,n)
else ! new reference grid point else ! new reference grid point
! PS: the old version fails if the pressure is <=1000 hPa in the whole ! PS: the old version fails if the pressure is <=1000 hPa in the whole
! domain. Let us find a good replacement, not just a quick fix. ! domain. Let us find a good replacement, not just a quick fix.
...@@ -185,11 +185,11 @@ subroutine verttransform_ecmwf(n,uuh,vvh,wwh,pvh) ...@@ -185,11 +185,11 @@ subroutine verttransform_ecmwf(n,uuh,vvh,wwh,pvh)
!> (reject mountains, accept valleys) !> (reject mountains, accept valleys)
psmean = sum( ps(:,:,1,n), ps(:,:,1,n) > psmean - psstd ) / & psmean = sum( ps(:,:,1,n), ps(:,:,1,n) > psmean - psstd ) / &
count(ps(:,:,1,n) > psmean - psstd) count(ps(:,:,1,n) > psmean - psstd)
print*,'height: std, new psmean',psstd,psmean ! print*,'height: std, new psmean',psstd,psmean
ixyref = minloc( abs( ps(:,:,1,n) - psmean ) ) ixyref = minloc( abs( ps(:,:,1,n) - psmean ) )
ixref = ixyref(1) ixref = ixyref(1)
jyref = ixyref(2) jyref = ixyref(2)
print*,'newheights at' ,ixref,jyref,ps(ixref,jyref,1,n) ! print*,'newheights at' ,ixref,jyref,ps(ixref,jyref,1,n)
endif endif
tvoldref=tt2(ixref,jyref,1,n)* & tvoldref=tt2(ixref,jyref,1,n)* &
......
Supports Markdown
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment