Changeset 5601 for branches/2015/dev_r5003_MERCATOR6_CRS/NEMOGCM/NEMO/OPA_SRC/ZDF/zdfmxl.F90

Timestamp:

2015-07-16T11:04:29+02:00 (9 years ago)

Author:

cbricaud

Message:

commit changes/bugfix/... for crs ; ok with time-splitting/fixed volume

File:

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

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

@@ -83 +83 @@
       REAL(wp) ::   zN2_c        ! local scalar
       INTEGER, POINTER, DIMENSION(:,:) ::   imld   ! 2D workspace
+      REAL(wp), POINTER, DIMENSION(:,:) ::   z2d   ! 2D workspace
       !!----------------------------------------------------------------------
       !
@@ -88 +89 @@
       !
       CALL wrk_alloc( jpi,jpj, imld )
+      CALL wrk_alloc( jpi,jpj, z2d )
 
       IF( kt == nit000 ) THEN
@@ -135 +137 @@
          CALL iom_put( "mldr10_1", hmlp )   ! mixed layer depth
          CALL iom_put( "mldkz5"  , hmld )   ! turbocline depth
+         z2d(:,:)=REAL(nmln,wp) 
+         CALL iom_put( "nmln"  , z2d )   ! turbocline depth
       ENDIF
       
@@ -140 +144 @@
       !
       CALL wrk_dealloc( jpi,jpj, imld )
+      CALL wrk_dealloc( jpi,jpj, z2d )
       !
       IF( nn_timing == 1 )  CALL timing_stop('zdf_mxl')