Changeset 2037 for branches

Timestamp:

2010-08-02T12:01:37+02:00 (14 years ago)

Author:

rfurner

Message:

bug fix to sbcrnf, jpkm1 changed to mbathy-1

File:

• branches/DEV_R1821_Rivers/NEMO/OPA_SRC/SBC/sbcrnf.F90 (modified) (1 diff)

branches/DEV_R1821_Rivers/NEMO/OPA_SRC/SBC/sbcrnf.F90

@@ -230 +230 @@
                 IF ( rnf_dep(ji,jj) > 0.e0 ) THEN  
                   jk=2  
-                  DO WHILE ( jk/=jpkm1 .AND. fsdept(ji,jj,jk) < rnf_dep(ji,jj) );  jk=jk+1;   ENDDO  
+                  DO WHILE ( jk/=(mbathy(ji,jj)-1) .AND. fsdept(ji,jj,jk) < rnf_dep(ji,jj) );  jk=jk+1;   ENDDO  
                   rnf_mod_dep(ji,jj)=jk  
                 ELSE IF ( rnf_dep(ji,jj) .eq. -1 ) THEN  
                   rnf_mod_dep(ji,jj)=1  
                 ELSE IF ( rnf_dep(ji,jj) .eq. -999 ) THEN  
-                  rnf_mod_dep(ji,jj)=jpkm1  
+                  rnf_mod_dep(ji,jj)=mbathy(ji,jj)-1
                 ELSE IF ( rnf_dep(ji,jj) /= 0 ) THEN  
                   CALL ctl_stop( 'runoff depth not positive, and not -999 or -1, rnf value in file fort.999'  )