Changeset 3952 for branches/2013/dev_r3406_CNRS_LIM3/NEMOGCM/NEMO/OPA_SRC/LDF/ldfdyn_c3d.h90

Timestamp:

2013-07-02T15:20:13+02:00 (11 years ago)

Author:

flavoni

Message:

ok for LIM2, still reproducibility error on LIM3 case

File:

• branches/2013/dev_r3406_CNRS_LIM3/NEMOGCM/NEMO/OPA_SRC/LDF/ldfdyn_c3d.h90 (modified) (3 diffs)

branches/2013/dev_r3406_CNRS_LIM3/NEMOGCM/NEMO/OPA_SRC/LDF/ldfdyn_c3d.h90

@@ -35 +35 @@
       REAL(wp) ::   zd_max       ! maximum grid spacing over the global domain
       REAL(wp) ::   za00, zc, zd, zetmax, zefmax, zeumax, zevmax   ! local scalars
-      REAL(wp), POINTER, DIMENSION(:)   :: zcoef   
-      INTEGER,  PARAMETER               :: zrim = 5  ! number of grid points over which diffusion is increased linearly
-      REAL(wp), PARAMETER               :: zinc_coef = 0.0 ! coef of diffusion increase  
-      INTEGER  ::   ii0, ii1, ij0, ij1   ! temporary integers
+      REAL(wp), POINTER, DIMENSION(:) :: zcoef   
       !!----------------------------------------------------------------------
       !
@@ -91 +88 @@
             CALL ldf_dyn_c3d_orca( ld_print )
          ENDIF
-    
-          ! Control print
-         IF(lwp .AND. ld_print ) THEN
-            WRITE(numout,*)
-            WRITE(numout,*) '         3D ahm1 array (k=1)'
-            CALL prihre( ahm1(:,:,1), jpi, jpj, 1, jpi, 1, 1, jpj, 1, 1.e-3, numout )
-            WRITE(numout,*)
-            WRITE(numout,*) '         3D ahm2 array (k=1)'
-            CALL prihre( ahm2(:,:,1), jpi, jpj, 1, jpi, 1, 1, jpj, 1, 1.e-3, numout )
-         ENDIF
-
-      ENDIF
-      
+
+      ENDIF
+      
+      ! Control print
+      IF(lwp .AND. ld_print ) THEN
+         WRITE(numout,*)
+         WRITE(numout,*) '         3D ahm1 array (k=1)'
+         CALL prihre( ahm1(:,:,1), jpi, jpj, 1, jpi, 20, 1, jpj, 20, 1.e-3, numout )
+         WRITE(numout,*)
+         WRITE(numout,*) '         3D ahm2 array (k=1)'
+         CALL prihre( ahm2(:,:,1), jpi, jpj, 1, jpi, 20, 1, jpj, 20, 1.e-3, numout )
+      ENDIF
+
 
       ! ahm3 and ahm4 at U- and V-points (used for bilaplacian operator
@@ -175 +172 @@
             WRITE(numout,*)
             WRITE(numout,*) 'inildf: ahm3 array at level 1'
-            CALL prihre(ahm3(:,:,1  ),jpi,jpj,1,jpi,1,1,jpj,1,1.e-12,numout)
+            CALL prihre(ahm3(:,:,1  ),jpi,jpj,1,jpi,1,1,jpj,1,1.e-3,numout)
             WRITE(numout,*)
             WRITE(numout,*) 'inildf: ahm4 array at level 1'
-            CALL prihre(ahm4(:,:,jpk),jpi,jpj,1,jpi,1,1,jpj,1,1.e-12,numout)
+            CALL prihre(ahm4(:,:,jpk),jpi,jpj,1,jpi,1,1,jpj,1,1.e-3,numout)
          ENDIF
       ENDIF