Changeset 1870 for branches/DEV_r1837_MLF/NEMO/OPA_SRC/TRA/trasbc.F90

Timestamp:

2010-05-12T17:36:00+02:00 (14 years ago)

Author:

gm

Message:

ticket: #663 step-1 : introduce the modified forcing term

File:

• branches/DEV_r1837_MLF/NEMO/OPA_SRC/TRA/trasbc.F90 (modified) (5 diffs)

branches/DEV_r1837_MLF/NEMO/OPA_SRC/TRA/trasbc.F90

@@ -4 +4 @@
    !! Ocean active tracers:  surface boundary condition
    !!==============================================================================
-   !! History :  8.2  !  98-10  (G. Madec, G. Roullet, M. Imbard)  Original code
-   !!            8.2  !  01-02  (D. Ludicone)  sea ice and free surface
-   !!            8.5  !  02-06  (G. Madec)  F90: Free form and module
+   !! History :  OPA  !  2998-10  (G. Madec, G. Roullet, M. Imbard)  Original code
+   !!            8.2  !  2001-02  (D. Ludicone)  sea ice and free surface
+   !!   NEMO     1.0  !  2002-06  (G. Madec)  F90: Free form and module
+   !!            3.3  !  2010-04  (M. Leclair, G. Madec)  Forcing averaged over 2 time steps
    !!----------------------------------------------------------------------
 
@@ -31 +32 @@
 #  include "vectopt_loop_substitute.h90"
    !!----------------------------------------------------------------------
-   !!   OPA 9.0 , LOCEAN-IPSL (2005) 
+   !! NEMO/OPA 3.3 , LOCEAN-IPSL (2010) 
    !! $Id$
    !! Software governed by the CeCILL licence (modipsl/doc/NEMO_CeCILL.txt)
@@ -46 +47 @@
       !!      and add it to the general trend of tracer equations.
       !!
-      !! ** Method :
-      !!      Following Roullet and Madec (2000), the air-sea flux can be divided 
-      !!      into three effects: (1) Fext, external forcing; 
+      !! ** Method  :   Following Roullet and Madec (2000), the air-sea flux 
+      !!      can be divided into three effects: 
+      !!      (1) Fext, external forcing; 
       !!      (2) Fwi, concentration/dilution effect due to water exchanged 
       !!         at the surface by evaporation, precipitations and runoff (E-P-R); 
@@ -58 +59 @@
       !!         solar part is added in traqsr.F routine.
       !!            ta = ta + q /(rau0 rcp e3t)  for k=1
-      !!            - salinity    : no salt flux
+      !!            - salinity    : only salt flux exchanged with sea-ice
       !!
       !!      The formulation for Fwb and Fwi vary according to the free 
@@ -123 +124 @@
       ENDIF
 
-      IF( .NOT.ln_traqsr )   qsr(:,:) = 0.e0   ! no solar radiation penetration
+!!gm      IF( .NOT.ln_traqsr )   qsr(:,:) = 0.e0   ! no solar radiation penetration
+      IF( .NOT.ln_traqsr ) THEN     ! no solar radiation penetration
+         qns(:,:) = qns(:,:) + qsr(:,:)      ! total heat flux in qns
+         qsr(:,:) = 0.e0                     ! qsr set to zero
+         IF( kt == nit000 ) THEN             ! idem on before field at nit000
+            qns_b(:,:) = qns_b(:,:) + qsr_b(:,:)  
+            qsr_b(:,:) = 0.e0                     
+         ENDIF
+      ENDIF
 
       ! Concentration dillution effect on (t,s)
-      DO jj = 2, jpj
-         DO ji = fs_2, fs_jpim1   ! vector opt.
-#if ! defined key_zco
-            zse3t = 1. / fse3t(ji,jj,1)
-#endif
-            IF( lk_vvl) THEN
-               zta = ro0cpr * qns(ji,jj) * zse3t &                   ! temperature : heat flux
-                &    - emp(ji,jj) * zsrau * tn(ji,jj,1)  * zse3t     ! & cooling/heating effet of EMP flux
-               zsa = 0.e0                                            ! No salinity concent./dilut. effect
-            ELSE
-               zta = ro0cpr * qns(ji,jj) * zse3t     ! temperature : heat flux
-               zsa = emps(ji,jj) * zsrau * sn(ji,jj,1)   * zse3t     ! salinity :  concent./dilut. effect
-            ENDIF
-            ta(ji,jj,1) = ta(ji,jj,1) + zta                          ! add the trend to the general tracer trend
-            sa(ji,jj,1) = sa(ji,jj,1) + zsa
+!!      DO jj = 2, jpj
+!!         DO ji = fs_2, fs_jpim1   ! vector opt.
+!!#if ! defined key_zco
+!!            zse3t = 1. / fse3t(ji,jj,1)
+!!#endif
+!!            IF( lk_vvl) THEN
+!!               zta = ro0cpr * qns(ji,jj) * zse3t &                   ! temperature : heat flux
+!!                &    - emp(ji,jj) * zsrau * tn(ji,jj,1)  * zse3t     ! & cooling/heating effet of EMP flux
+!!               zsa = 0.e0                                            ! No salinity concent./dilut. effect
+!!            ELSE
+!!               zta = ro0cpr * qns(ji,jj) * zse3t     ! temperature : heat flux
+!!               zsa = emps(ji,jj) * zsrau * sn(ji,jj,1)   * zse3t     ! salinity :  concent./dilut. effect
+!!            ENDIF
+!!            ta(ji,jj,1) = ta(ji,jj,1) + zta                          ! add the trend to the general tracer trend
+!!            sa(ji,jj,1) = sa(ji,jj,1) + zsa
+!!         END DO
+!!      END DO
+
+      !                             ! ---------------------- !
+      IF( lk_vvl ) THEN             !  Variable Volume case  !
+         !                          ! ---------------------- !
+         DO jj = 2, jpj
+            DO ji = fs_2, fs_jpim1   ! vector opt.
+               zse3t = 0.5 / fse3t(ji,jj,1)
+               !                           ! temperature: heat flux + cooling/heating effet of EMP flux
+               ta(ji,jj,1) = ta(ji,jj,1) + (  ro0cpr * ( qns(ji,jj) + qns_b(ji,jj) )                & 
+                  &                         - zsrau  * ( emp(ji,jj) + emp_b(ji,jj) ) * tn(ji,jj,1)  ) * zse3t
+               !                           ! salinity: salt flux
+               sa(ji,jj,1) = sa(ji,jj,1) + ( emps(ji,jj) + emps_b(ji,jj) ) * zse3t
+
+!!gm BUG : in key_vvl emps must be modified to only include the salt flux due to with sea-ice freezing/melting
+!!gm       otherwise this flux will be missing  ==> modification required in limsbc,  limsbc_2 and CICE interface.
+
+            END DO
          END DO
-      END DO
+         !                          ! ---------------------- !
+      ELSE                          !  Constant Volume case  !
+         !                          ! ---------------------- !
+         DO jj = 2, jpj
+            DO ji = fs_2, fs_jpim1   ! vector opt.
+               zse3t = 0.5 / fse3t(ji,jj,1)
+               !                           ! temperature: heat flux 
+               ta(ji,jj,1) = ta(ji,jj,1) + ro0cpr * ( qns (ji,jj) + qns_b (ji,jj) )               * zse3t 
+               !                           ! salinity: salt flux + concent./dilut. effect (both in emps)
+               sa(ji,jj,1) = sa(ji,jj,1) + zsrau  * ( emps(ji,jj) + emps_b(ji,jj) ) * sn(ji,jj,1) * zse3t
+            END DO
+         END DO
+         !
+      ENDIF
 
-      IF( l_trdtra ) THEN      ! save the sbc trends for diagnostic
+      IF( l_trdtra ) THEN        ! save the sbc trends for diagnostic
          ztrdt(:,:,:) = ta(:,:,:) - ztrdt(:,:,:)
          ztrds(:,:,:) = sa(:,:,:) - ztrds(:,:,:)
          CALL trd_mod(ztrdt, ztrds, jptra_trd_nsr, 'TRA', kt)
       ENDIF
-      !
+      !                          ! control print
       IF(ln_ctl)   CALL prt_ctl( tab3d_1=ta, clinfo1=' sbc  - Ta: ', mask1=tmask,   &
          &                       tab3d_2=sa, clinfo2=       ' Sa: ', mask2=tmask, clinfo3='tra' )