Changeset 11362 for NEMO/branches/2019/dev_r10984_HPC-13_IRRMANN_BDY_optimization/src/ICE/icedyn.F90

Timestamp:

2019-07-29T13:54:58+02:00 (5 years ago)

Author:

clem

Message:

debug the ice output by adding a missing value to the outputed fields. Unfortunately this method imposes a value of 1.e20 where there is no ice, which is annoying when using ncview for example, but this is the only way (from what I know) to output averages

File:

• NEMO/branches/2019/dev_r10984_HPC-13_IRRMANN_BDY_optimization/src/ICE/icedyn.F90 (modified) (4 diffs)

NEMO/branches/2019/dev_r10984_HPC-13_IRRMANN_BDY_optimization/src/ICE/icedyn.F90

@@ -75 +75 @@
       !!
       INTEGER  ::   ji, jj        ! dummy loop indices
+      REAL(wp) ::   zmiss_val     ! missing value retrieved from xios 
       REAL(wp) ::   zcoefu, zcoefv
-      REAL(wp), ALLOCATABLE, DIMENSION(:,:) ::   zdivu_i
+      REAL(wp), ALLOCATABLE, DIMENSION(:,:) ::   zdivu_i, zmsk00
       !!--------------------------------------------------------------------
       !
@@ -88 +89 @@
       ENDIF
       !                      
+      ! get missing value from xml
+      CALL iom_miss_val( "icethic", zmiss_val )
+      !
       ! retrieve thickness from volume for landfast param. and UMx advection scheme
       WHERE( a_i(:,:,:) >= epsi20 )
@@ -155 +159 @@
          CASE ( np_dynADV1D , np_dynADV2D )
 
-            ALLOCATE( zdivu_i(jpi,jpj) )
+            ALLOCATE( zdivu_i(jpi,jpj), zmsk00(jpi,jpj) )
+
+            DO jj = 1, jpj
+               DO ji = 1, jpi
+                  zmsk00(ji,jj) = MAX( 0._wp , SIGN( 1._wp , at_i(ji,jj) - epsi06  ) ) ! 1 if ice    , 0 if no ice
+               END DO
+            END DO
             DO jj = 2, jpjm1
                DO ji = 2, jpim1
@@ -163 +173 @@
             END DO
             CALL lbc_lnk( 'icedyn', zdivu_i, 'T', 1. )
-            CALL iom_put( "icediv" , zdivu_i(:,:) )
-            DEALLOCATE( zdivu_i )
+            ! output
+            CALL iom_put( "icediv" , zdivu_i * zmsk00 + zmiss_val * ( 1._wp - zmsk00 ) )
+
+            DEALLOCATE( zdivu_i, zmsk00 )
 
          END SELECT