Changeset 440 for trunk/NEMO/OPA_SRC/SBC/ocesbc.F90

Timestamp:

2006-04-19T16:43:17+02:00 (18 years ago)

Author:

opalod

Message:

nemo_v1_update_046 : CT : - light modifications related to the way to make SST/SSS damping

• add a flx_init subroutine to read the namflx namelist to get feedback coefficients for SST(dqdt0)/SSS(deds0) damping
• replace the ztrp variable by the dqdt0 read in the namflx namelist

File:

• trunk/NEMO/OPA_SRC/SBC/ocesbc.F90 (modified) (20 diffs)

trunk/NEMO/OPA_SRC/SBC/ocesbc.F90

@@ -87 +87 @@
       !! * Local declarations
       INTEGER ::   ji, jj                   ! dummy loop indices
-      REAL(wp) ::   ztx, ztaux, zty, ztauy, ztrp
+      REAL(wp) ::   ztx, ztaux, zty, ztauy
       REAL(wp) ::   ztdta, ztgel, zqrp
       !!----------------------------------------------------------------------
@@ -113 +113 @@
          ! Surface heat flux (W/m2)
          ! -----------------------
-         ztrp = 0.e0
 
          ! restoring heat flux
@@ -124 +123 @@
                ztdta = MAX( t_dta(ji,jj,1), ztgel )
 #endif
-               zqrp = ztrp * ( tb(ji,jj,1) - ztdta )
+               zqrp = dqdt0 * ( tb(ji,jj,1) - ztdta )
 
                qrp(ji,jj) = (1.0-freeze(ji,jj) ) * zqrp
@@ -306 +305 @@
       !! * Local declarations
       INTEGER  ::   ji, jj, jf         ! dummy loop indices
-      REAL(wp) ::   ztrp, ztgel,    &  ! temporary scalars
+      REAL(wp) ::   ztgel,          &  ! temporary scalars
          zice, zhemis, zqrp, zqri,  &  !    "         "
          zq, zqi, zerp, ze, zei, zro   !    "         "
@@ -313 +312 @@
       ! Compute fluxes
       ! --------------
-
-      ! constant initialization
-      ztrp = -40.   ! restoring term for temperature (w/m2/k)
 
       DO jj = 1, jpj
@@ -343 +339 @@
 
             zqrp = 0.
-            zqri = ztrp*MIN( 0., tb(ji,jj,1)-ztgel )
+            zqri = dqdt0*MIN( 0., tb(ji,jj,1)-ztgel )
             qrp(ji,jj) = ( ( 1. - zice ) * zqrp + zice * zqri ) * tmask(ji,jj,1)
 
@@ -417 +413 @@
       REAL(wp) ::  zxy
       REAL(wp) ::  zsice, zqri, zqrp, ztdta, zqrj
-      REAL(wp) ::  zq, zqi, zhemis, ztrp
+      REAL(wp) ::  zq, zqi, zhemis
       REAL(wp), DIMENSION(jpi,jpj) :: zeri, zerps, ziclim
       REAL(wp), DIMENSION(jpi,jpj) :: zqt, zqsr, zemp  
@@ -456 +452 @@
          CALL oce_sbc_dmp
 
-         ztrp = -40.             ! restoring terme for temperature (w/m2/k)   
          zsice = - 0.04 / 0.8    ! ratio of isohaline compressibility over isotherme compressibility 
                                  ! ( d rho / dt ) / ( d rho / ds )      ( s = 34, t = -1.8 )
@@ -486 +481 @@
                ! - gamma*min(0,t-tgel) if climatological and opa ice (ziclim=1 zicopa=1)
 
-               zqri = ztrp * ( tb(ji,jj,1) - ( fzptn(ji,jj) - 1.) )
-               zqrj = ztrp * MIN( 0., tb(ji,jj,1) - fzptn(ji,jj) )
+               zqri = dqdt0 * ( tb(ji,jj,1) - ( fzptn(ji,jj) - 1.) )
+               zqrj = dqdt0 * MIN( 0., tb(ji,jj,1) - fzptn(ji,jj) )
 
                qrp(ji,jj) =  ( ziclim(ji,jj) * ( (1 - freeze(ji,jj)) * zqri    &
@@ -493 +488 @@
 
 #if ! defined key_flx_bulk_monthly || ! defined key_flx_bulk_daily
-               zqrp = ztrp * ( tb(ji,jj,1) - ztdta )
+               zqrp = dqdt0 * ( tb(ji,jj,1) - ztdta )
                qrp(ji,jj) = qrp(ji,jj) + (1. - ziclim(ji,jj)) * zqrp
 # endif
@@ -579 +574 @@
       USE flxrnf                       ! ocean runoffs
       USE daymod, ONLY : nyear         ! calendar
+      USE dtasss                       ! sea surface salinity data
 
       !! * arguments
@@ -588 +584 @@
          qsr0 = 0.e0,               &  ! solar heat flux
          emp0 = 0.e0                   ! net freshwater flux
-      REAL(wp) ::   ztrp, zemp_S, zemp_N, zemp_sais, zTstar, zcos_sais1, zconv, zcos_sais2
+      REAL(wp) ::   zsrp,           &
+         zemp_S, zemp_N, zemp_sais, &
+         zTstar, zcos_sais1, zconv, &
+         zcos_sais2
       REAL(wp) ::           &
          zsumemp,           &          ! tampon used for the emp sum
@@ -643 +642 @@
           zcos_sais2 = COS( (ztime - ztimemax2) / (ztimemax2 - ztimemin2) * rpi )
 
-          ztrp= - 40.          ! retroaction term (W/m2/K)
           zconv = 3.16e-5      ! convert 1m/yr->3.16e-5mm/s
           DO jj = 1, jpj
@@ -653 +651 @@
                    &             * COS( rpi * (gphit(ji,jj) - 5.)        &                  
                    &                        / (53.5 * ( 1 + 11 / 53.5 * zcos_sais2 ) * 2.) ) 
-                qt  (ji,jj) = ztrp * ( tb(ji,jj,1) - t_star(ji,jj) )
+                qt  (ji,jj) = dqdt0 * ( tb(ji,jj,1) - t_star(ji,jj) )
                 IF( gphit(ji,jj) >= 14.845 .AND. 37.2 >= gphit(ji,jj)) THEN
                    ! zero at 37.8 deg, max at 24.6 deg
@@ -693 +691 @@
           ENDIF
 
-          !salinity terms
+          IF( lk_dtasss ) THEN           ! Sea surface salinity damping
+             ! deds0 from (mm/day) => (mm/s)
+             zsrp = deds0 / rday 
+             erp(:,:) = zsrp * ( sss(:,:) - sb(:,:,1) ) / ( sn(:,:,1) + 1.e-20 )
+          ELSE                           ! NO Sea surface salinity damping 
+             zsrp = 0.e0
+             erp(:,:) = 0.e0
+          ENDIF
+
+          ! freshwater fluxes terms
           DO jj = 1, jpj
              DO ji = 1, jpi
-                erp(ji,jj) = 0.e0 
                 emp(ji, jj) =  emp(ji,jj) - zsumemp * tmask(ji,jj,1)
-                emps(ji, jj)=  emp(ji,jj)
+                emps(ji, jj)=  emp(ji,jj) + erp(ji,jj)
              END DO
           END DO
@@ -722 +728 @@
             WRITE(numout,*)'           zsumemp    = ',zsumemp
             WRITE(numout,*)'           zsurf      = ',zsurf 
-            WRITE(numout,*)'           ztrp       = ',ztrp
+            WRITE(numout,*)'           dqdt0      = ',dqdt0
+            WRITE(numout,*)'           deds0      = ',deds0
             WRITE(numout,*)'           zconv      = ',zconv 
          ENDIF
@@ -728 +735 @@
       ELSE
 
+         ! Constant surface fluxes
+
          IF( kt == nit000 ) THEN
-
-            ! Read Namelist namflx : surface thermohaline fluxes
-            ! --------------------
-            REWIND ( numnam )
-            READ   ( numnam, namflx )
-   
-            IF(lwp) WRITE(numout,*)' '
-            IF(lwp) WRITE(numout,*)' ocesbc  : Constant surface fluxes read in namelist'
-            IF(lwp) WRITE(numout,*)' ~~~~~~~ '
-            IF(lwp) WRITE(numout,*)'           Namelist namflx: set the constant flux values'
-            IF(lwp) WRITE(numout,*)'              net heat flux          q0   = ', q0  , ' W/m2'
-            IF(lwp) WRITE(numout,*)'              solar heat flux        qsr0 = ', qsr0, ' W/m2'
-            IF(lwp) WRITE(numout,*)'              net heat flux          emp0 = ', emp0, ' W/m2'
+            IF(lwp) THEN
+               WRITE(numout,*)' '
+               WRITE(numout,*)' ocesbc  : Constant surface fluxes read in namelist'
+               WRITE(numout,*)' ~~~~~~~ '
+               WRITE(numout,*)'           Namelist namflx: set the constant flux values'
+               WRITE(numout,*)'              net heat flux          q0   = ', q0  , ' W/m2'
+               WRITE(numout,*)'              solar heat flux        qsr0 = ', qsr0, ' W/m2'
+               WRITE(numout,*)'              net heat flux          emp0 = ', emp0, ' W/m2'
+            ENDIF
 
             qt    (:,:) = q0
@@ -752 +757 @@
             runoff(:,:) = 0.e0
          ENDIF
+
       ENDIF
 
@@ -777 +783 @@
       INTEGER ::   ji, jj                   ! dummy loop indices
       REAL(wp), DIMENSION(jpi,jpj)  :: zsss, zfreeze
-      REAL(wp) ::   zerp, ztrp, zsrp
+      REAL(wp) ::   zerp, zsrp
       CHARACTER (len=71) :: charout
 #if ! defined key_dynspg_rl
@@ -807 +813 @@
       ! --------------
       ! Restoring coefficients on SST and SSS   
-      IF( lk_cpl ) THEN
-         ztrp = 0.e0
-         zsrp = 0.e0
-      ELSE
-         ztrp  = -40.                   ! (W/m2/K)    
-         zsrp  = ztrp * ro0cpr * rauw   ! (Kg/m2/s2) 
-      ENDIF
+      zsrp = dqdt0 * ro0cpr * rauw   ! (Kg/m2/s) 
 
 #if ! defined key_dynspg_rl