Changeset 521 for trunk/NEMO/OPA_SRC/TRD/trdmld.F90

Timestamp:

2006-10-11T16:45:09+02:00 (18 years ago)

Author:

opalod

Message:

nemo_v1_update_73 : CT : build Mixed Layer restart files using iom

File:

• trunk/NEMO/OPA_SRC/TRD/trdmld.F90 (modified) (4 diffs)

trunk/NEMO/OPA_SRC/TRD/trdmld.F90

@@ -35 +35 @@
    USE lbclnk          ! ocean lateral boundary conditions (or mpp link)
    USE diadimg         ! dimg direct access file format output
-   USE trdmld_rst , ONLY : trd_mld_rst_read  ! restart for diagnosing the ML trends
+   USE trdmld_rst      ! restart for diagnosing the ML trends
    USE prtctl          ! Print control
 
@@ -230 +230 @@
       !!
       INTEGER :: ji, jj, jk, jl, ik, it
-      LOGICAL :: lldebug = .TRUE.
+      LOGICAL :: lldebug = .FALSE.
       REAL(wp) :: zavt, zfn, zfn2
       REAL(wp) ,DIMENSION(jpi,jpj) ::     &
@@ -253 +253 @@
       ! ... These terms can be estimated by flux computation at the lower boundary of the ML 
       !     (we compute (-1/h) * K_z * d_z( T ) and (-1/h) * K_z * d_z( S ))
-      DO jj = 1,jpj
-         DO ji = 1,jpi
-            ik = nmld(ji,jj)
-            zavt = avt(ji,jj,ik)
-            tmltrd(ji,jj,jpmld_zdf) = - zavt / fse3w(ji,jj,ik) * tmask(ji,jj,ik)  &
-               &                      * ( tn(ji,jj,ik-1) - tn(ji,jj,ik) )         &
-               &                      / MAX( 1., rmld(ji,jj) ) * tmask(ji,jj,1)
-            zavt = fsavs(ji,jj,ik)
-            smltrd(ji,jj,jpmld_zdf) = - zavt / fse3w(ji,jj,ik) * tmask(ji,jj,ik)  &
-               &                      * ( sn(ji,jj,ik-1) - sn(ji,jj,ik) )         &
-               &                      / MAX( 1., rmld(ji,jj) ) * tmask(ji,jj,1)
-         END DO
-      END DO
-
-      ! ... Remove this K_z trend from the iso-neutral diffusion term (if any)
       IF( ln_traldf_iso ) THEN
+         DO jj = 1,jpj
+            DO ji = 1,jpi
+               ik = nmld(ji,jj)
+               zavt = avt(ji,jj,ik)
+               tmltrd(ji,jj,jpmld_zdf) = - zavt / fse3w(ji,jj,ik) * tmask(ji,jj,ik)  &
+                  &                      * ( tn(ji,jj,ik-1) - tn(ji,jj,ik) )         &
+                  &                      / MAX( 1., rmld(ji,jj) ) * tmask(ji,jj,1)
+               zavt = fsavs(ji,jj,ik)
+               smltrd(ji,jj,jpmld_zdf) = - zavt / fse3w(ji,jj,ik) * tmask(ji,jj,ik)  &
+                  &                      * ( sn(ji,jj,ik-1) - sn(ji,jj,ik) )         &
+                  &                      / MAX( 1., rmld(ji,jj) ) * tmask(ji,jj,1)
+            END DO
+         END DO
+
+         ! ... Remove this K_z trend from the iso-neutral diffusion term (if any)
          tmltrd(:,:,jpmld_ldf) = tmltrd(:,:,jpmld_ldf) - tmltrd(:,:,jpmld_zdf)
          smltrd(:,:,jpmld_ldf) = smltrd(:,:,jpmld_ldf) - smltrd(:,:,jpmld_zdf)
@@ -698 +698 @@
       END IF
 
+      ! ======================================================================
+      ! V. Write restart file
+      ! ======================================================================
+
+      CALL trd_mld_rst_write( kt ) 
+
    END SUBROUTINE trd_mld