Changeset 6772 for branches/2015/dev_r5003_MERCATOR6_CRS/NEMOGCM/NEMO/OPA_SRC/ZDF/zdfmxl_crs.F90

Timestamp:

2016-07-01T18:02:45+02:00 (8 years ago)

Author:

cbricaud

Message:

clean in coarsening branch

File:

• branches/2015/dev_r5003_MERCATOR6_CRS/NEMOGCM/NEMO/OPA_SRC/ZDF/zdfmxl_crs.F90 (modified) (4 diffs)

branches/2015/dev_r5003_MERCATOR6_CRS/NEMOGCM/NEMO/OPA_SRC/ZDF/zdfmxl_crs.F90

@@ -24 +24 @@
    USE trc_oce, ONLY : lk_offline ! offline flag
    USE crs
+   USE ieee_arithmetic
 
    IMPLICIT NONE
@@ -64 +65 @@
       REAL(wp) ::   zN2_c        ! local scalar
       INTEGER, POINTER, DIMENSION(:,:) ::   imld   ! 2D workspace
-      REAL(wp), POINTER, DIMENSION(:,:) ::  z2d   ! 2D workspace
       !!----------------------------------------------------------------------
       !
@@ -70 +70 @@
       !
       CALL wrk_alloc( jpi_crs,jpj_crs, imld )
-      CALL wrk_alloc( jpi_crs,jpj_crs, z2d )
 
       IF( kt == nit000 ) THEN
@@ -95 +94 @@
             iiki = imld(ji,jj)
             iikn = nmln_crs(ji,jj)
-            IF( iiki .NE. 0 ) hmld_crs (ji,jj) = ( gdepw_crs(ji,jj,iiki  ) - gdepw_crs(ji,jj,1 )            )   * tmask_crs(ji,jj,1)  ! Turbocline depth 
-            IF( iiki .NE. 0 ) hmlpt_crs(ji,jj) = ( gdept_crs(ji,jj,iikn-1) - gdepw_crs(ji,jj,1 )            )   * tmask_crs(ji,jj,1)  ! depth of the last T-point inside the mixed layer
+            hmld_crs (ji,jj) = ( fsdepw_crs(ji,jj,iiki  ) - fsdepw_crs(ji,jj,1    )  ) * tmask_crs(ji,jj,1)    ! Turbocline depth 
+            hmlp_crs (ji,jj) = ( fsdepw_crs(ji,jj,iikn  ) - fsdepw_crs(ji,jj,nla10)  ) * tmask_crs(ji,jj,1)    ! Mixed layer depth
+            hmlpt_crs(ji,jj) = ( fsdept_crs(ji,jj,iikn-1) - fsdepw_crs(ji,jj,1    )  ) * tmask_crs(ji,jj,1)    ! depth of the last T-point inside the mixed layer
          END DO
       END DO
       !
-      z2d=REAL(nmln_crs,wp)
-      CALL iom_put("nmln_crs",z2d)
-      CALL iom_put("hmlpt_crs",hmlpt_crs)
-      !
       CALL wrk_dealloc( jpi_crs,jpj_crs, imld )
-      CALL wrk_dealloc( jpi_crs,jpj_crs, z2d )
       !
       IF( nn_timing == 1 )  CALL timing_stop('zdf_mxl_crs')