Changeset 1312


Ignore:
Timestamp:
2009-02-16T17:35:36+01:00 (12 years ago)
Author:
smasson
Message:

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

Location:
trunk
Files:
21 edited

Legend:

Unmodified
Added
Removed
  • trunk/CONFIG/GYRE/EXP00/namelist

    r1278 r1312  
    3939   nwrite      =      60   !  frequency of write in the output file  
    4040   ln_dimgnnn  = .false.   !  DIMG file format: 1 file for all processors (F) or by processor (T) 
     41   ln_mskland  = .false.   !  mask land points (1.e+20) in NetCDF outputs (costly: + ~15%) 
    4142/ 
    4243!!====================================================================== 
  • trunk/CONFIG/GYRE_LOBSTER/EXP00/namelist

    r1278 r1312  
    3939   nwrite      =     360   !  frequency of write in the output file  
    4040   ln_dimgnnn  = .false.   !  DIMG file format: 1 file for all processors (F) or by processor (T) 
     41   ln_mskland  = .false.   !  mask land points in NetCDF outputs (costly: + ~15%) 
    4142/ 
    4243!!====================================================================== 
  • trunk/CONFIG/ORCA2_LIM/EXP00/1_namelist

    r1278 r1312  
    3939   nwrite      =   10950   !  frequency of write in the output file  
    4040   ln_dimgnnn  = .false.   !  DIMG file format: 1 file for all processors (F) or by processor (T) 
     41   ln_mskland  = .false.   !  mask land points in NetCDF outputs (costly: + ~15%) 
    4142/ 
    4243!!====================================================================== 
  • trunk/CONFIG/ORCA2_LIM/EXP00/namelist

    r1278 r1312  
    3939   nwrite      =    5475   !  frequency of write in the output file  
    4040   ln_dimgnnn  = .false.   !  DIMG file format: 1 file for all processors (F) or by processor (T) 
     41   ln_mskland  = .false.   !  mask land points in NetCDF outputs (costly: + ~15%) 
    4142/ 
    4243!!====================================================================== 
  • trunk/CONFIG/ORCA2_LIM_PISCES/EXP00/namelist

    r1278 r1312  
    3939   nwrite      =    5475   !  frequency of write in the output file  
    4040   ln_dimgnnn  = .false.   !  DIMG file format: 1 file for all processors (F) or by processor (T) 
     41   ln_mskland  = .false.   !  mask land points in NetCDF outputs (costly: + ~15%) 
    4142/ 
    4243!!====================================================================== 
  • trunk/CONFIG/ORCA2_OFF_PISCES/EXP00/namelist

    r1292 r1312  
    3030   nstock      =     7300  !  frequency of creation of a restart file 
    3131   nwrite      =     7300  !  frequency of write in the output file  
     32   ln_mskland  = .false.   !  mask land points in NetCDF outputs (costly: + ~15%) 
    3233/ 
    3334!----------------------------------------------------------------------- 
  • trunk/NEMO/C1D_SRC/diawri_c1d.F90

    r1310 r1312  
    9090      zdt = rdt 
    9191      IF( nacc == 1 ) zdt = rdtmin 
     92      IF( ln_mskland )   THEN   ;   clop = "only(x)"   ! put 1.e+20 on land (very expensive!!) 
     93      ELSE                      ;   clop = "x"         ! no use of the mask value (require less cpu time) 
     94      ENDIF 
    9295#if defined key_diainstant 
    9396      zsto = nwrite * zdt 
    94       clop = "inst(x)"           ! no use of the mask value (require less cpu time) 
    95       !!! clop="inst(only(x))"   ! put 1.e+20 on land (very expensive!!) 
     97      clop = "inst("//TRIM(clop)//")" 
    9698#else 
    9799      zsto=zdt 
    98       clop="ave(x)"              ! no use of the mask value (require less cpu time) 
    99       !!! clop="ave(only(x))"    ! put 1.e+20 on land (very expensive!!) 
     100      clop = "ave("//TRIM(clop)//")" 
    100101#endif 
    101102      zout = nwrite * zdt 
  • trunk/NEMO/LIM_SRC_2/limwri_2.F90

    r1310 r1312  
    9999                            
    100100         zsto     = rdt_ice 
    101 !!Chris         clop     = "ave(only(x))"      !ibug  namelist parameter a ajouter 
    102          clop     = "ave(x)" 
     101         IF( ln_mskland )   THEN   ;   clop = "ave(only(x))"   ! put 1.e+20 on land (very expensive!!) 
     102         ELSE                      ;   clop = "ave(x)"         ! no use of the mask value (require less cpu time) 
     103         ENDIF 
    103104         zout     = nwrite * rdt_ice / nn_fsbc 
    104105         niter    = 0 
  • trunk/NEMO/LIM_SRC_3/limwri.F90

    r1310 r1312  
    137137 
    138138         zsto     = rdt_ice 
    139          clop     = "ave(x)" 
     139         IF( ln_mskland )   THEN   ;   clop = "ave(only(x))"   ! put 1.e+20 on land (very expensive!!) 
     140         ELSE                      ;   clop = "ave(x)"         ! no use of the mask value (require less cpu time) 
     141         ENDIF 
    140142         zout     = nwrite * rdt_ice / nn_fsbc 
    141143         niter    = 0 
  • trunk/NEMO/OFF_SRC/DOM/domain.F90

    r1152 r1312  
    110110      NAMELIST/namrun/ no    , cexper   , ln_rstart , nrstdt , nit000,         & 
    111111         &             nitend, ndate0   , nleapy    , ninist , nstock,         & 
    112          &             nwrite  
     112         &             nwrite, ln_mskland  
    113113 
    114114      NAMELIST/namctl/ ln_ctl , nprint, nictls, nictle,   & 
     
    147147         WRITE(numout,*) '           frequency of restart file       nstock    = ', nstock 
    148148         WRITE(numout,*) '           frequency of output file        nwrite    = ', nwrite 
     149         WRITE(numout,*) '           mask land points             ln_mskland   = ', ln_mskland 
    149150      ENDIF 
    150151 
  • trunk/NEMO/OFF_SRC/IOM/in_out_manager.F90

    r1152 r1312  
    3535   LOGICAL            ::   ln_dimgnnn = .FALSE.   !: type of dimgout. (F): 1 file for all proc 
    3636                                                  !:                  (T): 1 file per proc 
     37   LOGICAL            ::   ln_mskland = .FALSE.   !: mask land points in NetCDF outputs (costly: + ~15%) 
    3738   !!---------------------------------------------------------------------- 
    3839   !! was in restart but moved here because of the OFF line... better solution should be found... 
  • trunk/NEMO/OPA_SRC/DIA/diagap.F90

    r1310 r1312  
    153153         ! Define frequency of output and means 
    154154         zsto = ngap * zdt 
    155          clop = "ave(x)" 
     155         IF( ln_mskland )   THEN   ;   clop = "ave(only(x))"   ! put 1.e+20 on land (very expensive!!) 
     156         ELSE                      ;   clop = "ave(x)"         ! no use of the mask value (require less cpu time) 
     157         ENDIF 
    156158         zout = ngap * zdt 
    157159         zmax = FLOAT( nitend - nit000 + 1 ) * zdt 
  • trunk/NEMO/OPA_SRC/DIA/diaptr.F90

    r1310 r1312  
    426426      zdt = rdt 
    427427      IF( nacc == 1 ) zdt = rdtmin 
     428      IF( ln_mskland )   THEN   ;   clop = "only(x)"   ! put 1.e+20 on land (very expensive!!) 
     429      ELSE                      ;   clop = "x"         ! no use of the mask value (require less cpu time) 
     430      ENDIF 
    428431#if defined key_diainstant 
    429          zsto = nf_ptr * zdt 
    430          clop = "inst(x)"               ! no use of the mask value (require less cpu time) 
    431          !!! clop="inst(only(x))"       ! put 1.e+20 on land (very expensive!!) 
     432      zsto = nf_ptr * zdt 
     433      clop = "inst("//TRIM(clop)//")" 
    432434#else 
    433          zsto = zdt 
    434          clop = "ave(x)"                ! no use of the mask value (require less cpu time) 
    435          !!! clop="ave(only(x))"        ! put 1.e+20 on land (very expensive!!) 
     435      zsto=zdt 
     436      clop = "ave("//TRIM(clop)//")" 
    436437#endif 
    437438      zout = nf_ptr * zdt 
  • trunk/NEMO/OPA_SRC/DIA/diawri.F90

    r1310 r1312  
    129129      zdt = rdt 
    130130      IF( nacc == 1 ) zdt = rdtmin 
     131      IF( ln_mskland )   THEN   ;   clop = "only(x)"   ! put 1.e+20 on land (very expensive!!) 
     132      ELSE                      ;   clop = "x"         ! no use of the mask value (require less cpu time) 
     133      ENDIF 
    131134#if defined key_diainstant 
    132135      zsto = nwrite * zdt 
    133       clop = "inst(x)"           ! no use of the mask value (require less cpu time) 
    134       !!! clop="inst(only(x))"   ! put 1.e+20 on land (very expensive!!) 
     136      clop = "inst("//TRIM(clop)//")" 
    135137#else 
    136138      zsto=zdt 
    137       clop="ave(x)"              ! no use of the mask value (require less cpu time) 
    138       !!! clop="ave(only(x))"    ! put 1.e+20 on land (very expensive!!) 
     139      clop = "ave("//TRIM(clop)//")" 
    139140#endif 
    140141      zout = nwrite * zdt 
  • trunk/NEMO/OPA_SRC/DOM/domain.F90

    r1241 r1312  
    142142      NAMELIST/namrun/ no    , cexper, cn_ocerst_in, cn_ocerst_out, ln_rstart, nrstdt,   & 
    143143         &             nit000, nitend, ndate0      , nleapy       , ninist   , nstock,   & 
    144          &             nwrite, ln_dimgnnn 
     144         &             nwrite, ln_dimgnnn, ln_mskland 
    145145 
    146146      NAMELIST/namdom/ ntopo , e3zps_min, e3zps_rat, nmsh   ,   & 
     
    174174         WRITE(numout,*) '           frequency of output file        nwrite    = ', nwrite 
    175175         WRITE(numout,*) '           multi file dimgout           ln_dimgnnn   = ', ln_dimgnnn 
     176         WRITE(numout,*) '           mask land points             ln_mskland   = ', ln_mskland 
    176177      ENDIF 
    177178 
  • trunk/NEMO/OPA_SRC/IOM/in_out_manager.F90

    r1239 r1312  
    3939   LOGICAL            ::   ln_dimgnnn    = .FALSE.     !: type of dimgout. (F): 1 file for all proc 
    4040                                                       !:                  (T): 1 file per proc 
     41   LOGICAL            ::   ln_mskland    = .FALSE.     !: mask land points in NetCDF outputs (costly: + ~15%) 
    4142   !!---------------------------------------------------------------------- 
    4243   !! was in restart but moved here because of the OFF line... better solution should be found... 
  • trunk/NEMO/OPA_SRC/TRD/trdmld.F90

    r1310 r1312  
    824824      ! II.1 Define frequency of output and means 
    825825      ! ----------------------------------------- 
     826      IF( ln_mskland )   THEN   ;   clop = "only(x)"   ! put 1.e+20 on land (very expensive!!) 
     827      ELSE                      ;   clop = "x"         ! no use of the mask value (require less cpu time) 
     828      ENDIF 
    826829#  if defined key_diainstant 
    827830      IF( .NOT. ln_trdmld_instant ) THEN 
     
    829832      END IF 
    830833      zsto = ntrd * rdt 
    831       clop ="inst(only(x))" 
     834      clop = "inst("//TRIM(clop)//")" 
    832835#  else 
    833836      IF( ln_trdmld_instant ) THEN 
     
    836839         zsto = ntrd * rdt          ! mean  diags : we DO NOT use any IOIPSL time averaging 
    837840      END IF 
    838       clop ="ave(only(x))" 
     841      clop = "ave("//TRIM(clop)//")" 
    839842#  endif 
    840843      zout = ntrd * rdt 
  • trunk/NEMO/OPA_SRC/TRD/trdvor.F90

    r1310 r1312  
    539539      ! II.1 Define frequency of output and means 
    540540      ! ----------------------------------------- 
     541      IF( ln_mskland )   THEN   ;   clop = "only(x)"   ! put 1.e+20 on land (very expensive!!) 
     542      ELSE                      ;   clop = "x"         ! no use of the mask value (require less cpu time) 
     543      ENDIF 
    541544#if defined key_diainstant 
    542545      zsto = nwrite*rdt 
    543       clop ="inst(x)" 
     546      clop = "inst("//TRIM(clop)//")" 
    544547#else 
    545548      zsto = rdt 
    546       clop ="ave(x)" 
     549      clop = "ave("//TRIM(clop)//")" 
    547550#endif 
    548551      zout = ntrd*rdt 
  • trunk/NEMO/TOP_SRC/SED/sedwri.F90

    r1310 r1312  
    5858      ! Define frequency of output and means 
    5959      zdt = dtsed 
     60      IF( ln_mskland )   THEN   ;   clop = "only(x)"   ! put 1.e+20 on land (very expensive!!) 
     61      ELSE                      ;   clop = "x"         ! no use of the mask value (require less cpu time) 
     62      ENDIF 
    6063#if defined key_diainstant 
    6164      zsto = nwrised * zdt 
    62 !!      clop = "inst(x)"           ! no use of the mask value (require less cpu time) 
    63       clop="inst(only(x))"   ! put 1.e+20 on land (very expensive!!) 
     65      clop = "inst("//TRIM(clop)//")" 
    6466#else 
    6567      zsto = zdt 
    66 !!      clop = "ave(x)"              ! no use of the mask value (require less cpu time) 
    67       clop="ave(only(x))"    ! put 1.e+20 on land (very expensive!!) 
     68      clop = "ave("//TRIM(clop)//")" 
    6869#endif 
    6970      zout = nwrised * zdt 
  • trunk/NEMO/TOP_SRC/TRP/trdmld_trc.F90

    r1310 r1312  
    13911391      ! II.1 Define frequency of output and means 
    13921392      ! ----------------------------------------- 
     1393      IF( ln_mskland )   THEN   ;   clop = "only(x)"   ! put 1.e+20 on land (very expensive!!) 
     1394      ELSE                      ;   clop = "x"         ! no use of the mask value (require less cpu time) 
     1395      ENDIF 
    13931396#  if defined key_diainstant 
    13941397      IF( .NOT. ln_trdmld_trc_instant ) THEN 
     
    13961399      ENDIF 
    13971400      zsto = ntrd_trc * rdt 
    1398       clop ="inst(only(x))" 
     1401      clop = "inst("//TRIM(clop)//")" 
    13991402#  else 
    14001403      IF( ln_trdmld_trc_instant ) THEN 
     
    14031406         zsto = ntrd_trc * rdt                                    ! mean  diags : we DO NOT use any IOIPSL time averaging 
    14041407      ENDIF 
    1405       clop ="ave(only(x))" 
     1408      clop = "ave("//TRIM(clop)//")" 
    14061409#  endif 
    14071410      zout = ntrd_trc * rdt 
  • trunk/NEMO/TOP_SRC/trcdia.F90

    r1310 r1312  
    123123      ! Define frequency of output and means 
    124124      zdt = rdt 
     125      IF( ln_mskland )   THEN   ;   clop = "only(x)"   ! put 1.e+20 on land (very expensive!!) 
     126      ELSE                      ;   clop = "x"         ! no use of the mask value (require less cpu time) 
     127      ENDIF 
    125128# if defined key_diainstant 
    126129      zsto = nwritetrc * rdt 
    127       clop = 'inst(only(x))' 
     130      clop = "inst("//TRIM(clop)//")" 
    128131# else 
    129132      zsto = zdt 
    130       clop = 'ave(only(x))' 
     133      clop = "ave("//TRIM(clop)//")" 
    131134# endif 
    132135      zout = nwritetrc * zdt 
Note: See TracChangeset for help on using the changeset viewer.