Changeset 1756 for trunk/NEMO/LIM_SRC_2/limsbc_2.F90

Timestamp:

2009-11-25T15:15:20+01:00 (14 years ago)

Author:

smasson

Message:

implement AR5 diagnostics, see ticket:610

File:

• trunk/NEMO/LIM_SRC_2/limsbc_2.F90 (modified) (8 diffs)

trunk/NEMO/LIM_SRC_2/limsbc_2.F90

@@ -24 +24 @@
    USE lbclnk           ! ocean lateral boundary condition
    USE in_out_manager   ! I/O manager
+   USE diaar5, ONLY :   lk_diaar5
    USE iom              ! 
    USE albedo           ! albedo parameters
@@ -80 +81 @@
       INTEGER  ::   iflt, ial, iadv, ifral, ifrdv
       INTEGER  ::   ii0, ii1, ij0, ij1  ! temporary integers
+      REAL(wp) ::   zrdtir           ! 1. / rdt_ice
       REAL(wp) ::   zqsr  , zqns     ! solar & non solar heat flux
       REAL(wp) ::   zinda            ! switch for testing the values of ice concentration
@@ -97 +99 @@
       !!---------------------------------------------------------------------
      
+      zrdtir = 1. / rdt_ice
+      
       IF( kt == nit000 ) THEN
          IF(lwp) WRITE(numout,*)
@@ -180 +184 @@
             zqns    =  - ( 1. - thcm(ji,jj) ) * zqsr   &   ! part of the solar energy used in leads
                &       + iflt    * ( fscmbq(ji,jj) + ffltbif(ji,jj) )                            &
-               &       + ifral   * ( ial * qcmif(ji,jj) + (1 - ial) * qldif(ji,jj) ) / rdt_ice   &
-               &       + ifrdv   * ( qfvbq(ji,jj) + qdtcn(ji,jj) ) / rdt_ice
+               &       + ifral   * ( ial * qcmif(ji,jj) + (1 - ial) * qldif(ji,jj) ) * zrdtir    &
+               &       + ifrdv   * ( qfvbq(ji,jj) + qdtcn(ji,jj) ) * zrdtir
 
             fsbbq(ji,jj) = ( 1.0 - ( ifvt + iflt ) ) * fscmbq(ji,jj)     ! ???
@@ -190 +194 @@
       END DO
 
+      CALL iom_put( 'hflx_ice_cea', - fdtcn(:,:) )      
       CALL iom_put( 'qns_io_cea', qns(:,:) - zqnsoce(:,:) * pfrld(:,:) )      
       CALL iom_put( 'qsr_io_cea', fstric(:,:) * (1. - pfrld(:,:)) )
@@ -206 +211 @@
 #if defined key_coupled
           zemp = emp_tot(ji,jj) - emp_ice(ji,jj) * ( 1. - pfrld(ji,jj) )    &   ! 
-             &   + rdmsnif(ji,jj) / rdt_ice                                     !  freshwaterflux due to snow melting 
+             &   + rdmsnif(ji,jj) * zrdtir                                      !  freshwaterflux due to snow melting 
 #else
 !!$            !  computing freshwater exchanges at the ice/ocean interface
@@ -217 +222 @@
                &   - tprecip(ji,jj) * ( 1. -  frld(ji,jj) )    &   !  liquid precipitation reaches directly the ocean
                &   + sprecip(ji,jj) * ( 1. - pfrld(ji,jj) )    &   !  taking into account change in ice cover within the time step
-               &   + rdmsnif(ji,jj) / rdt_ice                      !  freshwaterflux due to snow melting 
+               &   + rdmsnif(ji,jj) * zrdtir                       !  freshwaterflux due to snow melting 
                !                                                   !  ice-covered fraction:
 #endif            
 
             !  computing salt exchanges at the ice/ocean interface
-            zfons =  ( soce_r(ji,jj) - sice_r(ji,jj) ) * ( rdmicif(ji,jj) / rdt_ice ) 
+            zfons =  ( soce_r(ji,jj) - sice_r(ji,jj) ) * ( rdmicif(ji,jj) * zrdtir ) 
             
             !  converting the salt flux from ice to a freshwater flux from ocean
@@ -232 +237 @@
          END DO
       END DO
+
+      IF( lk_diaar5 ) THEN
+         CALL iom_put( 'isnwmlt_cea'  ,                 rdmsnif(:,:) * zrdtir )
+         CALL iom_put( 'fsal_virt_cea',   soce_r(:,:) * rdmicif(:,:) * zrdtir )
+         CALL iom_put( 'fsal_real_cea', - sice_r(:,:) * rdmicif(:,:) * zrdtir )
+      ENDIF
 
       !------------------------------------------!