Changeset 1177


Ignore:
Timestamp:
2008-09-11T18:27:52+02:00 (12 years ago)
Author:
cetlod
Message:

update modules to take into account new trends organization, see ticket:248

Location:
trunk/NEMO/TOP_SRC
Files:
4 edited

Legend:

Unmodified
Added
Removed
  • trunk/NEMO/TOP_SRC/oce_trc.F90

    r1152 r1177  
    154154   USE oce , ONLY :   rhop    =>    rhop    !: potential volumic mass (kg m-3)  
    155155   USE oce , ONLY :   rhd     =>    rhd     !: in situ density anomalie rhd=(rho-rau0)/rau0 (no units) 
    156 # if defined key_trc_diatrd 
    157156   USE oce , ONLY :   hdivn   =>    hdivn   !: horizontal divergence (1/s) 
    158 # endif 
    159157 
    160158 
     
    206204 
    207205   !* mixing & mixed layer depth * 
     206   USE zdfmxl , ONLY :   nmln        =>   nmln        !: number of level in the mixed layer 
    208207   USE zdfmxl , ONLY :   hmld        =>   hmld        !: mixing layer depth (turbocline) 
    209208   USE zdfmxl , ONLY :   hmlp        =>   hmlp        !: mixed layer depth  (rho=rho0+zdcrit) (m) 
  • trunk/NEMO/TOP_SRC/trcini.F90

    r1146 r1177  
    2626   USE trcini_my_trc   ! MY_TRC  initialisation 
    2727   USE trcdta          !  
     28   USE trdmld_trc_oce 
     29   USE trdmld_trc 
     30 
    2831   USE in_out_manager  ! I/O manager 
    2932   USE prtctl_trc      ! Print control passive tracers (prt_ctl_trc_init routine) 
     
    118121      ENDIF 
    119122 
     123      IF( lk_trdmld_trc  )  CALL trd_mld_trc_init     ! trends: Mixed-layer 
     124 
    120125      !                 ! Computation content of all tracers 
    121126      trai = 0.e0 
  • trunk/NEMO/TOP_SRC/trclec.F90

    r1152 r1177  
    2727   USE trclsm_my_trc     ! MY_TRC SMS namelist 
    2828   USE in_out_manager    ! I/O manager 
     29   USE trdmld_trc_oce 
    2930 
    3031   IMPLICIT NONE 
     
    6970      !! 
    7071      NAMELIST/namtoptrc/ ndttrc, nwritetrc, lrsttr, nrsttr, tracer 
    71 #if defined key_trc_diatrd 
    72       NAMELIST/namtoptrd/ luttrd, nwritetrd                                             ! dynamical trends 
     72#if defined key_trdmld_trc || defined key_trc_diatrd 
     73      NAMELIST/namtoptrd/ ntrd_trc, nctls_trc, ucf_trc, & 
     74                         ln_trdmld_trc_restart, ln_trdmld_trc_instant, & 
     75                         luttrd 
    7376#endif 
    7477 
     
    131134      ENDIF 
    132135 
    133 #if defined key_trc_diatrd 
     136#if defined key_trdmld_trc || defined key_trc_diatrd 
    134137 
    135138      ! Namelist natrtd (transport trends) 
    136139      ! ---------------------------------------------- 
    137       nwritetrd = 10                 ! default values (no dynamical trend recording) 
     140      ntrd_trc  = 20 
     141      nctls_trc =  6 
     142      ucf_trc   =  1. 
     143      ln_trdmld_trc_instant = .TRUE. 
     144      ln_trdmld_trc_restart =.FALSE. 
    138145      DO jn = 1, jptra 
    139146         luttrd(jn) = .FALSE. 
    140147      END DO 
    141148 
    142       REWIND( numnat )               ! read natrtd 
     149      REWIND( numnat )               !  namelist namtoptrd : passive tracer trends diagnostic 
    143150      READ  ( numnat, namtoptrd ) 
    144151 
    145       nkeep=0 
    146       ikeep(:)=0 
     152#if defined key_trc_diatrd 
     153      nkeep    = 0 
     154      ikeep(:) = 0 
    147155      DO jn = 1, jptra 
    148156         IF( luttrd(jn) ) THEN  
    149              nkeep    = nkeep + 1 
    150              ikeep(jn)=nkeep 
     157             nkeep     = nkeep + 1 
     158             ikeep(jn) = nkeep 
    151159         END IF  
    152160      END DO 
     
    156164      ENDIF  
    157165 
    158       IF(lwp) THEN                   ! control print 
     166      nwritetrd = ntrd_trc 
     167#endif 
     168 
     169     IF(lwp) THEN 
    159170         WRITE(numout,*) 
    160          WRITE(numout,*) ' Namelist : namtoptrd' 
    161          WRITE(numout,*) '    frequency of outputs for dynamical trends nwritetrd = ', nwritetrd 
     171         WRITE(numout,*) ' trd_mld_trc_init : read namelist namtoptrd                        ' 
     172         WRITE(numout,*) ' ~~~~~~~~~~~~~~~~                                               ' 
     173         WRITE(numout,*) '   * frequency of trends diagnostics   ntrd_trc               = ', ntrd_trc 
     174         WRITE(numout,*) '   * control surface type              nctls_trc              = ', nctls_trc 
     175         WRITE(numout,*) '   * restart for ML diagnostics        ln_trdmld_trc_restart  = ', ln_trdmld_trc_restart 
     176         WRITE(numout,*) '   * flag to diagnose trends of                                 ' 
     177         WRITE(numout,*) '     instantantaneous or mean ML T/S   ln_trdmld_trc_instant  = ', ln_trdmld_trc_instant 
     178         WRITE(numout,*) '   * unit conversion factor            ucf_trc                = ', ucf_trc 
     179#if defined key_trdmld_trc 
     180         DO jn = 1, jptra 
     181            WRITE(numout,*) '    keep dynamical trends for tracer number :', jn, luttrd(jn) 
     182         END DO 
     183#else 
    162184         DO jn = 1, jptra 
    163185            WRITE(numout,*) '    keep dynamical trends for tracer number :', jn, luttrd(jn), ikeep(jn) 
    164186         END DO 
     187 
     188         WRITE(numout,*) 
    165189         WRITE(numout,*) '    total = ', nkeep, ' tracers dyn trends saved' 
    166190         WRITE(numout,*) '    size of trtrd = ', jpi*jpj*jpk*nkeep*jpdiatrc 
    167       ENDIF 
     191#endif 
     192      ENDIF 
     193 
    168194#endif 
    169195 
  • trunk/NEMO/TOP_SRC/trcrst.F90

    r1146 r1177  
    230230          CALL trc_rst_stat            ! statistics 
    231231          CALL iom_close( numrtw )     ! close the restart file (only at last time step) 
     232#if ! defined key_trdmld_trc 
    232233          lrst_trc = .FALSE. 
     234#endif 
    233235       ENDIF 
    234236      ! 
Note: See TracChangeset for help on using the changeset viewer.