We read every piece of feedback, and take your input very seriously.
To see all available qualifiers, see our documentation.
There was an error while loading. Please reload this page.
1 parent 5b24c88 commit 6615030Copy full SHA for 6615030
sorc/ncep_post.fd/INITPOST_NETCDF.f
@@ -2074,8 +2074,10 @@ SUBROUTINE INITPOST_NETCDF(ncid2d,ncid3d)
2074
!$omp parallel do private(i,j,tlmh)
2075
Do j=jsta,jend
2076
Do i=ista,iend
2077
- TLMH = T(I,J,LM) * T(I,J,LM)
2078
- Sigt4(i,j) = 5.67E-8 * TLMH * TLMH
+ if ( T(I,J,LM) < spval ) then
+ TLMH = T(I,J,LM) * T(I,J,LM)
2079
+ Sigt4(i,j) = 5.67E-8 * TLMH * TLMH
2080
+ endif
2081
End do
2082
2083
@@ -2624,7 +2626,9 @@ SUBROUTINE INITPOST_NETCDF(ncid2d,ncid3d)
2624
2626
if(ext550(i,j,l)<spval)then
2625
2627
taod5503d ( i, j, l) = ext550 ( i, j, l )
2628
dz = ZINT( i, j, l ) - ZINT( i, j, l+1 )
- aextc55 ( i, j, l ) = taod5503d ( i, j, l ) / dz
2629
+ if ( dz /= 0.0) then
2630
+ aextc55 ( i, j, l ) = taod5503d ( i, j, l ) / dz
2631
2632
endif
2633
if(debugprint.and.i==im/2.and.j==(jsta+jend)/2)print*,'sample taod5503d= ', &
2634
i,j,l,taod5503d ( i, j, l )
0 commit comments