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

Timestamp:

2015-12-17T16:48:41+01:00 (8 years ago)

Author:

cbricaud

Message:

correction of bugs from last update and improvments for CRS

File:

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

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

@@ -64 +64 @@
       REAL(wp) ::   zN2_c        ! local scalar
       INTEGER, POINTER, DIMENSION(:,:) ::   imld   ! 2D workspace
+      REAL(wp), POINTER, DIMENSION(:,:) ::  z2d   ! 2D workspace
       !!----------------------------------------------------------------------
       !
@@ -69 +70 @@
       !
       CALL wrk_alloc( jpi_crs,jpj_crs, imld )
+      CALL wrk_alloc( jpi_crs,jpj_crs, z2d )
 
       IF( kt == nit000 ) THEN
@@ -98 +100 @@
       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')