Changeset 1312 for trunk/NEMO/OPA_SRC/DIA/diaptr.F90

Timestamp:

2009-02-16T17:35:36+01:00 (15 years ago)

Author:

smasson

Message:

add a namelist logical to mask land points in NetCDF outputs, see ticket:322

File:

• trunk/NEMO/OPA_SRC/DIA/diaptr.F90 (modified) (1 diff)

trunk/NEMO/OPA_SRC/DIA/diaptr.F90

@@ -426 +426 @@
       zdt = rdt
       IF( nacc == 1 ) zdt = rdtmin
+      IF( ln_mskland )   THEN   ;   clop = "only(x)"   ! put 1.e+20 on land (very expensive!!)
+      ELSE                      ;   clop = "x"         ! no use of the mask value (require less cpu time)
+      ENDIF
 #if defined key_diainstant
-         zsto = nf_ptr * zdt
-         clop = "inst(x)"               ! no use of the mask value (require less cpu time)
-         !!! clop="inst(only(x))"       ! put 1.e+20 on land (very expensive!!)
+      zsto = nf_ptr * zdt
+      clop = "inst("//TRIM(clop)//")"
 #else
-         zsto = zdt
-         clop = "ave(x)"                ! no use of the mask value (require less cpu time)
-         !!! clop="ave(only(x))"        ! put 1.e+20 on land (very expensive!!)
+      zsto=zdt
+      clop = "ave("//TRIM(clop)//")"
 #endif
       zout = nf_ptr * zdt