Opened 16 years ago
Closed 16 years ago
#351 closed bug (Fixed)
Added missing CPP directives
| Reported by: | arango | Owned by: | arango |
|---|---|---|---|
| Priority: | major | Milestone: | Adjoint Based Algorithms |
| Component: | Adjoint | Version: | 3.3 |
| Keywords: | Cc: |
Description
Corrected missing CPP directives in routines tl_step3d_uv.F and ad_step3d_uv.F. This bug was introduced in src:ticket:346. Many thanks to Gregoire Broquet and Andy Moore for reporting this.
I also corrected a small defect in tl_balance.F and ad_balance.F. We need to have:
DO j=JstrR,JendR
DO i=IstrR,IendR
DO k=1,N(ng)
cff=0.5_r8*(dTdz(i,j,k-1)+dTdz(i,j,k))
IF (ABS(cff).lt.dTdz_min) THEN
dzdT=0.0_r8
ELSE
dzdT=1.0_r8/cff
END IF
dSdT(k)=(0.5_r8*(dSdz(i,j,k-1)+ &
& dSdz(i,j,k )))*dzdT
END DO
instead of:
DO j=JstrR,JendR
DO i=IstrR,IendR
DO k=1,N(ng)
IF (ABS(dTdz(i,j,k)).lt.dTdz_min) THEN
dzdT=0.0_r8
ELSE
dzdT=1.0_r8/(0.5_r8*(dTdz(i,j,k-1)+ &
& dTdz(i,j,k )))
END IF
dSdT(k)=(0.5_r8*(dSdz(i,j,k-1)+ &
& dSdz(i,j,k )))*dzdT
END DO
Note:
See TracTickets
for help on using tickets.
