Changeset 1316


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

continue changeset:1312, see ticke:322

Location:
trunk/NEMO
Files:
2 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 
  • trunk/NEMO/TOP_SRC/trcdia.F90

    r1312 r1316  
    263263      ! Define frequency of output and means 
    264264      zdt = rdt 
     265      IF( ln_mskland )   THEN   ;   clop = "only(x)"   ! put 1.e+20 on land (very expensive!!) 
     266      ELSE                      ;   clop = "x"         ! no use of the mask value (require less cpu time) 
     267      ENDIF 
    265268#  if defined key_diainstant 
    266269      zsto = nwritetrd * rdt 
    267       clop = 'inst(only(x))' 
     270      clop = "inst("//TRIM(clop)//")" 
    268271#  else 
    269272      zsto = zdt 
    270       clop = 'ave(only(x))' 
     273      clop = "ave("//TRIM(clop)//")" 
    271274#  endif 
    272275      zout = nwritetrd * zdt 
     
    506509      ! Define frequency of output and means 
    507510      zdt = rdt 
     511      IF( ln_mskland )   THEN   ;   clop = "only(x)"   ! put 1.e+20 on land (very expensive!!) 
     512      ELSE                      ;   clop = "x"         ! no use of the mask value (require less cpu time) 
     513      ENDIF 
    508514#  if defined key_diainstant 
    509515      zsto=nwritedia*zdt 
    510       clop='inst(only(x))' 
     516      clop = "inst("//TRIM(clop)//")" 
    511517#  else 
    512518      zsto=zdt 
    513       clop='ave(only(x))' 
     519      clop = "ave("//TRIM(clop)//")" 
    514520#  endif 
    515521      zout=nwritedia*zdt 
     
    656662      ! Define frequency of output and means 
    657663      zdt = rdt 
     664      IF( ln_mskland )   THEN   ;   clop = "only(x)"   ! put 1.e+20 on land (very expensive!!) 
     665      ELSE                      ;   clop = "x"         ! no use of the mask value (require less cpu time) 
     666      ENDIF 
    658667#        if defined key_diainstant 
    659668      zsto=nwritebio*zdt 
    660       clop='inst(only(x))' 
     669      clop = "inst("//TRIM(clop)//")" 
    661670#        else 
    662671      zsto=zdt 
    663       clop='ave(only(x))' 
     672      clop = "ave("//TRIM(clop)//")" 
    664673#        endif 
    665674      zout=nwritebio*zdt 
Note: See TracChangeset for help on using the changeset viewer.