diff --git a/sorc/ncep_post.fd/GFSPOST.F b/sorc/ncep_post.fd/GFSPOST.F index 823aaf7ed..ae43a9f1b 100644 --- a/sorc/ncep_post.fd/GFSPOST.F +++ b/sorc/ncep_post.fd/GFSPOST.F @@ -419,7 +419,7 @@ subroutine tpause(km,p,u,v,t,h,ptp,utp,vtp,ttp,htp,shrtp) integer klim(2),k,kd,ktp ! - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - ! find tropopause level - call rsearch1(km-2,p(2),2,ptplim(1),klim(1)) + call rsearch1(km-2,p(2),2,ptplim,klim) klim(1)=klim(1)+1 klim(2)=klim(2)+2 ! klim(1) > klim(2) or loops does not run ; klim(2) has a @@ -519,7 +519,7 @@ subroutine mxwind(km,p,u,v,t,h,pmw,umw,vmw,tmw,hmw) real spd(km),spdmw,wmw,dhd,dhu,shrd,shru,dhmw,ub,vb,spdb ! - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - ! find maximum wind level - call rsearch1(km,p(1),2,pmwlim(1),klim(1)) + call rsearch1(km,p(1),2,pmwlim,klim) ! klim(1)=klim(1)+1 klim(2)=klim(2)+1 ! spd(klim(1):klim(2))=sqrt(u(klim(1):klim(2))**2+v(klim(1):klim(2))**2) diff --git a/sorc/ncep_post.fd/INITPOST_GFS_NETCDF.f b/sorc/ncep_post.fd/INITPOST_GFS_NETCDF.f index 36557b75a..914606704 100644 --- a/sorc/ncep_post.fd/INITPOST_GFS_NETCDF.f +++ b/sorc/ncep_post.fd/INITPOST_GFS_NETCDF.f @@ -135,7 +135,7 @@ SUBROUTINE INITPOST_GFS_NETCDF(ncid3d) ! DECLARE VARIABLES. ! ! REAL fhour - integer nfhour ! forecast hour from nems io file +! integer nfhour ! forecast hour from nems io file integer fhzero !bucket real dtp !physics time step REAL RINC(5) @@ -573,7 +573,7 @@ SUBROUTINE INITPOST_GFS_NETCDF(ncid3d) deallocate(glat1d,glon1d) print*,'idate = ',(idate(i),i=1,7) - print*,'nfhour = ',nfhour +! print*,'nfhour = ',nfhour ! sample print point ii = im/2 diff --git a/sorc/ncep_post.fd/INITPOST_NETCDF.f b/sorc/ncep_post.fd/INITPOST_NETCDF.f index 9984b2388..72f51b64a 100644 --- a/sorc/ncep_post.fd/INITPOST_NETCDF.f +++ b/sorc/ncep_post.fd/INITPOST_NETCDF.f @@ -2795,7 +2795,7 @@ subroutine read_netcdf_3d_scatter(me,ncid,ifhr,im,jm,jsta,jsta_2l & ! iret = nf90_get_var(ncid,varid,dummy2,start=(/1,1,l,ifhr/), & ! count=(/im,jm,1,1/)) if (iret /= 0) then - print*,VarName,l," not found -Assigned missing values" + print*,VarName," not found -Assigned missing values" do l=1,lm do j=1,jm do i=1,im