Changeset 1938 for branches/DEV_R1821_Rivers/NEMO/OPA_SRC/DOM/closea.F90

Timestamp:

2010-06-16T16:34:29+02:00 (14 years ago)

Author:

rfurner

Message:

rnf has been separated from emp and emps. Also temperature and salinity of runoff can be specified, and runoff can be added to a user specified depth

File:

• branches/DEV_R1821_Rivers/NEMO/OPA_SRC/DOM/closea.F90 (modified) (7 diffs)

branches/DEV_R1821_Rivers/NEMO/OPA_SRC/DOM/closea.F90

@@ -180 +180 @@
       INTEGER                     ::   ji, jj, jc, jn   ! dummy loop indices
       REAL(wp)                    ::   zze2
-      REAL(wp), DIMENSION (jpncs) ::   zemp
+      REAL(wp), DIMENSION (jpncs) ::   zfwf 
+ 
+      REAL(wp), DIMENSION (jpi,jpj)       ::   fwf              !: fresh water flux,for volume: emp-rnf  
+      REAL(wp), DIMENSION (jpi,jpj)       ::   fwfs             !: fresh water flux, for conc/dilution effect: emps-rnf  
       !!----------------------------------------------------------------------
       !
@@ -216 +219 @@
       !                                                   !--------------------!
       !                                                   !  update emp, emps  !
-      zemp = 0.e0                                         !--------------------!
+      zfwf = 0.e0                                         !--------------------!
       DO jc = 1, jpncs
          DO jj = ncsj1(jc), ncsj2(jc)
             DO ji = ncsi1(jc), ncsi2(jc)
-               zemp(jc) = zemp(jc) + e1t(ji,jj) * e2t(ji,jj) * emp(ji,jj) * tmask_i(ji,jj)
+               zfwf(jc) = zfwf(jc) + e1t(ji,jj) * e2t(ji,jj) * ( emp(ji,jj)-rnf(ji,jj) ) * tmask_i(ji,jj) 
             END DO  
          END DO 
       END DO
-      IF( lk_mpp )   CALL mpp_sum ( zemp(:) , jpncs )       ! mpp: sum over all the global domain
+      IF( lk_mpp )   CALL mpp_sum ( zfwf(:) , jpncs )       ! mpp: sum over all the global domain
 
       IF( cp_cfg == "orca" .AND. jp_cfg == 2 ) THEN      ! Black Sea case for ORCA_R2 configuration
-         zze2    = ( zemp(3) + zemp(4) ) / 2.
-         zemp(3) = zze2
-         zemp(4) = zze2
+         zze2    = ( zfwf(3) + zfwf(4) ) / 2.
+         zfwf(3) = zze2
+         zfwf(4) = zze2
       ENDIF
 
@@ -236 +239 @@
          IF( ncstt(jc) == 0 ) THEN 
             ! water/evap excess is shared by all open ocean
-            emp (:,:) = emp (:,:) + zemp(jc) / surf(jpncs+1)
-            emps(:,:) = emps(:,:) + zemp(jc) / surf(jpncs+1)
+            emp (:,:) = emp (:,:) + zfwf(jc) / surf(jpncs+1)
+            emps(:,:) = emps(:,:) + zfwf(jc) / surf(jpncs+1)
          ELSEIF( ncstt(jc) == 1 ) THEN 
             ! Excess water in open sea, at outflow location, excess evap shared
-            IF ( zemp(jc) <= 0.e0 ) THEN 
+            IF ( zfwf(jc) <= 0.e0 ) THEN 
                 DO jn = 1, ncsnr(jc)
                   ji = mi0(ncsir(jc,jn))
@@ -246 +249 @@
                   IF (      ji > 1 .AND. ji < jpi   &
                       .AND. jj > 1 .AND. jj < jpj ) THEN 
-                      emp (ji,jj) = emp (ji,jj) + zemp(jc) /   &
+                      emp (ji,jj) = emp (ji,jj) + zfwf(jc) /   &
                          (FLOAT(ncsnr(jc)) * e1t(ji,jj) * e2t(ji,jj))
-                      emps(ji,jj) = emps(ji,jj) + zemp(jc) /   &
+                      emps(ji,jj) = emps(ji,jj) + zfwf(jc) /   &
                           (FLOAT(ncsnr(jc)) * e1t(ji,jj) * e2t(ji,jj))
                   END IF 
                 END DO 
             ELSE 
-                emp (:,:) = emp (:,:) + zemp(jc) / surf(jpncs+1)
-                emps(:,:) = emps(:,:) + zemp(jc) / surf(jpncs+1)
+                emp (:,:) = emp (:,:) + zfwf(jc) / surf(jpncs+1)
+                emps(:,:) = emps(:,:) + zfwf(jc) / surf(jpncs+1)
             ENDIF
          ELSEIF( ncstt(jc) == 2 ) THEN 
@@ -263 +266 @@
                   ji = mi0(ncsir(jc,jn))
                   jj = mj0(ncsjr(jc,jn)) ! Location of outflow in open ocean
-                  emp (ji,jj) = emp (ji,jj) + zemp(jc)   &
+                  emp (ji,jj) = emp (ji,jj) + zfwf(jc)   &
                       / (FLOAT(ncsnr(jc)) *  e1t(ji,jj) * e2t(ji,jj) )
-                  emps(ji,jj) = emps(ji,jj) + zemp(jc)   &
+                  emps(ji,jj) = emps(ji,jj) + zfwf(jc)   &
                       / (FLOAT(ncsnr(jc)) *  e1t(ji,jj) * e2t(ji,jj) )
                 END DO 
@@ -273 +276 @@
          DO jj = ncsj1(jc), ncsj2(jc)
             DO ji = ncsi1(jc), ncsi2(jc)
-               emp (ji,jj) = emp (ji,jj) - zemp(jc) / surf(jc)
-               emps(ji,jj) = emps(ji,jj) - zemp(jc) / surf(jc)
+               emp (ji,jj) = emp (ji,jj) - zfwf(jc) / surf(jc)
+               emps(ji,jj) = emps(ji,jj) - zfwf(jc) / surf(jc)
             END DO  
          END DO 
@@ -280 +283 @@
       END DO 
       !
-      CALL lbc_lnk( emp , 'T', 1. )
-      CALL lbc_lnk( emps, 'T', 1. )
+      fwf(:,:)  = emp(:,:)  - rnf(:,:)  
+      fwfs(:,:) = emps(:,:) - rnf(:,:)  
+      CALL lbc_lnk( fwf , 'T', 1. ) 
+      CALL lbc_lnk( fwfs, 'T', 1. ) 
       !
    END SUBROUTINE sbc_clo