New URL for NEMO forge!   http://forge.nemo-ocean.eu

Since March 2022 along with NEMO 4.2 release, the code development moved to a self-hosted GitLab.
This present forge is now archived and remained online for history.
Changeset 1316 for trunk/NEMO/OPA_SRC/DIA/diaptr.F90 – NEMO

Ignore:
Timestamp:
2009-02-17T16:12:46+01:00 (15 years ago)
Author:
smasson
Message:

continue changeset:1312, see ticke:322

File:
1 edited

Legend:

Unmodified
Added
Removed
  • trunk/NEMO/OPA_SRC/DIA/diaptr.F90

    r1312 r1316  
    419419      CHARACTER (len=40)       ::   clhstnam, clop                   ! temporary names 
    420420      INTEGER                  ::   iline, it, ji                    ! 
    421       REAL(wp)                 ::   zsto, zout, zdt, zmax, zjulian   ! temporary scalars 
     421      REAL(wp)                 ::   zsto, zout, zdt, zjulian   ! temporary scalars 
    422422      REAL(wp), DIMENSION(jpj) ::   zphi, zfoo 
    423423      !!---------------------------------------------------------------------- 
    424424       
    425       ! Define frequency of output and means 
    426       zdt = rdt 
    427       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 
    431 #if defined key_diainstant 
    432       zsto = nf_ptr * zdt 
    433       clop = "inst("//TRIM(clop)//")" 
    434 #else 
    435       zsto=zdt 
    436       clop = "ave("//TRIM(clop)//")" 
    437 #endif 
    438       zout = nf_ptr * zdt 
    439       zmax = ( nitend - nit000 + 1 ) * zdt 
    440        
    441           
    442425      ! define time axis 
    443426      it = kt - nit000 + 1 
     
    446429      ! -------------- 
    447430      IF( kt == nit000 ) THEN 
    448        
    449       zdt = rdt 
    450       IF( nacc == 1 ) zdt = rdtmin 
     431          
     432         zdt = rdt 
     433         IF( nacc == 1 ) zdt = rdtmin 
    451434 
    452435         ! Reference latitude 
     
    482465         ! Define frequency of output and means 
    483466         zsto = nf_ptr * zdt 
    484          clop = "ave(x)" 
     467         IF( ln_mskland )   THEN   ;   clop = "ave(only(x))"   ! put 1.e+20 on land (very expensive!!) 
     468         ELSE                      ;   clop = "ave(x)"         ! no use of the mask value (require less cpu time) 
     469         ENDIF 
    485470         zout = nf_ptr * zdt 
    486471         zfoo(:) = 0.e0 
Note: See TracChangeset for help on using the changeset viewer.