Changeset 732


Ignore:
Timestamp:
2007-11-19T16:19:15+01:00 (14 years ago)
Author:
cetlod
Message:

Update modules for passive tracers transport trends computation, see ticket:13

Location:
tags/nemo_v2_3/NEMO/TOP_SRC
Files:
9 edited

Legend:

Unmodified
Added
Removed
  • tags/nemo_v2_3/NEMO/TOP_SRC/TRP/trcadv_muscl.F90

    r501 r732  
    2727   !!---------------------------------------------------------------------- 
    2828   !!   TOP 1.0 , LOCEAN-IPSL (2005)  
    29    !! $Header$  
     29   !! $Header: /home/opalod/NEMOCVSROOT/NEMO/TOP_SRC/TRP/trcadv_muscl.F90,v 1.13 2007/10/12 09:26:30 opalod Exp $  
    3030   !! This software is governed by the CeCILL licence see modipsl/doc/NEMO_CeCILL.txt  
    3131   !!---------------------------------------------------------------------- 
     
    103103 
    104104      DO jn = 1, jptra 
    105  
     105#if defined key_trc_diatrd 
     106        DO jk = 1,jpk 
     107           DO jj = 1,jpj 
     108              DO ji = 1,jpi 
     109                 IF (luttrd(jn)) trtrd(ji,jj,jk,ikeep(jn),1) = 0. 
     110                 IF (luttrd(jn)) trtrd(ji,jj,jk,ikeep(jn),2) = 0. 
     111                 IF (luttrd(jn)) trtrd(ji,jj,jk,ikeep(jn),3) = 0. 
     112              END DO 
     113            END DO 
     114          END DO 
     115#endif 
    106116         ! I. Horizontal advective fluxes 
    107117         ! ------------------------------ 
  • tags/nemo_v2_3/NEMO/TOP_SRC/TRP/trcdmp.F90

    r494 r732  
    3737   !!---------------------------------------------------------------------- 
    3838   !!   TOP 1.0 , LOCEAN-IPSL (2005)  
    39    !! $Header$  
     39   !! $Header: /home/opalod/NEMOCVSROOT/NEMO/TOP_SRC/TRP/trcdmp.F90,v 1.12 2007/10/12 09:26:30 opalod Exp $  
    4040   !! This software is governed by the CeCILL licence see modipsl/doc/NEMO_CeCILL.txt  
    4141   !!---------------------------------------------------------------------- 
     
    108108#    if defined key_trc_diatrd 
    109109                        ! save the trends for diagnostics 
    110                         IF (luttrd(jn)) trtrd(ji,jj,jk,ikeep(jn),jpdiatrc) = ztra 
     110                        IF (luttrd(jn)) trtrd(ji,jj,jk,ikeep(jn),jpdiatrc-1) = ztra 
    111111#    endif 
    112112                     END DO 
     
    129129#    if defined key_trc_diatrd 
    130130                        ! save the trends for diagnostics 
    131                         IF (luttrd(jn)) trtrd(ji,jj,jk,ikeep(jn),jpdiatrc) = ztra 
     131                        IF (luttrd(jn)) trtrd(ji,jj,jk,ikeep(jn),jpdiatrc-1) = ztra 
    132132#    endif 
    133133                     END DO 
     
    149149#    if defined key_trc_diatrd 
    150150                        ! save the trends for diagnostics 
    151                         IF (luttrd(jn)) trtrd(ji,jj,jk,ikeep(jn),jpdiatrc) = ztra 
     151                        IF (luttrd(jn)) trtrd(ji,jj,jk,ikeep(jn),jpdiatrc-1) = ztra 
    152152#    endif 
    153153                     END DO 
  • tags/nemo_v2_3/NEMO/TOP_SRC/TRP/trcldf_bilap.F90

    r668 r732  
    2525   !!---------------------------------------------------------------------- 
    2626   !!   TOP 1.0 , LOCEAN-IPSL (2005)  
    27    !! $Header: /home/opalod/NEMOCVSROOT/NEMO/TOP_SRC/TRP/trcldf_bilap.F90,v 1.11 2007/05/28 02:51:54 opalod Exp $  
     27   !! $Header: /home/opalod/NEMOCVSROOT/NEMO/TOP_SRC/TRP/trcldf_bilap.F90,v 1.12 2007/10/12 09:26:30 opalod Exp $  
    2828   !! This software is governed by the CeCILL licence see modipsl/doc/NEMO_CeCILL.txt  
    2929   !!---------------------------------------------------------------------- 
     
    148148                     ikv = MIN ( mbathy(ji,jj), mbathy(ji  ,jj+1) ) - 1 
    149149                     IF( iku == jk ) THEN 
    150                         ztu(ji,jj,jk) = zeeu(ji,jj) * gtru(ji,jj) 
     150                        ztu(ji,jj,jk) = zeeu(ji,jj) * gtru(ji,jj,jn) 
    151151                     ENDIF 
    152152                     IF( ikv == jk ) THEN 
    153                         ztv(ji,jj,jk) = zeev(ji,jj) * gtrv(ji,jj) 
     153                        ztv(ji,jj,jk) = zeev(ji,jj) * gtrv(ji,jj,jn) 
    154154                     ENDIF 
    155155                  END DO 
  • tags/nemo_v2_3/NEMO/TOP_SRC/TRP/trcldf_iso.F90

    r433 r732  
    2727   !!---------------------------------------------------------------------- 
    2828   !!   TOP 1.0 , LOCEAN-IPSL (2005)  
    29    !! $Header$  
     29   !! $Header: /home/opalod/NEMOCVSROOT/NEMO/TOP_SRC/TRP/trcldf_iso.F90,v 1.10 2007/10/12 09:26:30 opalod Exp $  
    3030   !! This software is governed by the CeCILL licence see modipsl/doc/NEMO_CeCILL.txt  
    3131   !!---------------------------------------------------------------------- 
     
    9595      REAL(wp), DIMENSION(jpi,jpj) ::   & 
    9696         zftug, zftvg 
     97 
     98#if defined key_trc_diatrd 
     99      REAL(wp) ::   & 
     100         ztagu, ztagv 
     101#endif 
     102 
    97103#endif 
    98104 
     
    202208               END DO 
    203209            END DO 
     210 
     211#if defined key_trc_diatrd 
     212#   if defined key_trcldf_eiv 
     213            DO jj = 2 , jpjm1 
     214               DO ji = fs_2, fs_jpim1   ! vector opt. 
     215                  zbtr= 1. / ( e1t(ji,jj) * e2t(ji,jj) * fse3t(ji,jj,jk) ) 
     216                  ztagu = ( zftug(ji,jj) - zftug(ji-1,jj  ) ) * zbtr 
     217                  ztagv = ( zftvg(ji,jj) - zftvg(ji  ,jj-1) ) * zbtr 
     218                  IF (luttrd(jn)) trtrd (ji,jj,jk,ikeep(jn),4) = trtrd(ji,jj,jk,ikeep(jn),4) - ztagu 
     219                  IF (luttrd(jn)) trtrd (ji,jj,jk,ikeep(jn),5) = trtrd(ji,jj,jk,ikeep(jn),5) - ztagv 
     220                  IF (luttrd(jn)) trtrd (ji,jj,jk,ikeep(jn),7) = ztagu 
     221                  IF (luttrd(jn)) trtrd (ji,jj,jk,ikeep(jn),8) = ztagv 
     222               END DO 
     223             END DO 
     224#   endif 
     225#endif 
     226 
    204227            !                                          ! =============== 
    205228         END DO                                        !   End of slab   
  • tags/nemo_v2_3/NEMO/TOP_SRC/TRP/trcldf_iso_zps.F90

    r501 r732  
    2626   !!---------------------------------------------------------------------- 
    2727   !!   TOP 1.0 , LOCEAN-IPSL (2005)  
    28    !! $Header$  
     28   !! $Header: /home/opalod/NEMOCVSROOT/NEMO/TOP_SRC/TRP/trcldf_iso_zps.F90,v 1.11 2007/10/12 09:26:30 opalod Exp $  
    2929   !! This software is governed by the CeCILL licence see modipsl/doc/NEMO_CeCILL.txt  
    3030   !!---------------------------------------------------------------------- 
     
    8585      REAL(wp) ::   & 
    8686         zabe1, zabe2, zcof1, zcof2,   &  ! temporary scalars 
    87          zmsku, zmskv, zbtr, ztra,     & 
    88          ztagu, ztagv 
     87         zmsku, zmskv, zbtr, ztra        
    8988 
    9089      REAL(wp), DIMENSION(jpi,jpj) ::   & 
     
    102101         zuwk1, zvwk1,                 & 
    103102         zcg1,zcg2 
     103 
     104#if defined key_trc_diatrd 
     105      REAL(wp) ::   & 
     106         ztagu, ztagv 
     107#endif 
     108 
    104109#endif 
    105110      CHARACTER (len=22) :: charout 
     
    118123      DO jn = 1, jptra 
    119124 
     125#if defined key_trcldf_eiv && key_trc_diatrd 
    120126         ztagu = 0.e0 
    121127         ztagv = 0.e0 
     128#endif 
    122129 
    123130         ! Horizontal passive tracer gradient  
     
    234241                  IF (luttrd(jn)) trtrd (ji,jj,jk,ikeep(jn),4) = trtrd(ji,jj,jk,ikeep(jn),4) - ztagu 
    235242                  IF (luttrd(jn)) trtrd (ji,jj,jk,ikeep(jn),5) = trtrd(ji,jj,jk,ikeep(jn),5) - ztagv 
     243                  IF (luttrd(jn)) trtrd (ji,jj,jk,ikeep(jn),7) = ztagu 
     244                  IF (luttrd(jn)) trtrd (ji,jj,jk,ikeep(jn),8) = ztagv 
    236245               END DO 
    237246            END DO 
  • tags/nemo_v2_3/NEMO/TOP_SRC/TRP/trczdf_iso.F90

    r615 r732  
    3333   !!---------------------------------------------------------------------- 
    3434   !!   TOP 1.0 , LOCEAN-IPSL (2005)  
    35    !! $Header$  
     35   !! $Header: /home/opalod/NEMOCVSROOT/NEMO/TOP_SRC/TRP/trczdf_iso.F90,v 1.13 2007/10/12 09:26:30 opalod Exp $  
    3636   !! This software is governed by the CeCILL licence see modipsl/doc/NEMO_CeCILL.txt  
    3737   !!---------------------------------------------------------------------- 
     
    343343                  !  WARNING trtrd(ji,jj,jk,6) used for vertical gent velocity trend 
    344344                  !                           not for damping !!! 
    345                   IF (luttrd(jn)) trtrd(ji,jj,jk,ikeep(jn),6) = ztavg 
     345                  IF (luttrd(jn)) trtrd(ji,jj,jk,ikeep(jn),9) = ztavg 
    346346#   endif 
    347347                  IF (luttrd(jn)) trtrd(ji,jj,jk,ikeep(jn),6) = ztav - ztavg 
  • tags/nemo_v2_3/NEMO/TOP_SRC/TRP/trczdf_iso_vopt.F90

    r615 r732  
    3636   !!---------------------------------------------------------------------- 
    3737   !!   TOP 1.0 , LOCEAN-IPSL (2005)  
    38    !! $Header$  
     38   !! $Header: /home/opalod/NEMOCVSROOT/NEMO/TOP_SRC/TRP/trczdf_iso_vopt.F90,v 1.12 2007/10/12 09:26:30 opalod Exp $  
    3939   !! This software is governed by the CeCILL licence see modipsl/doc/NEMO_CeCILL.txt  
    4040   !!---------------------------------------------------------------------- 
     
    548548                  ztavg = ( zwx(ji,jj,jk) - zwx(ji,jj,jk+1) ) * zbtr 
    549549                  !  WARNING trtrd(ji,jj,jk,7) used for vertical gent velocity trend  not for damping !!! 
    550                   IF (luttrd(jn)) trtrd(ji,jj,jk,ikeep(jn),7) = ztavg 
     550                  IF (luttrd(jn)) trtrd(ji,jj,jk,ikeep(jn),9) = ztavg 
    551551#   endif 
    552552                  IF (luttrd(jn)) trtrd(ji,jj,jk,ikeep(jn),6) = ztav - ztavg 
  • tags/nemo_v2_3/NEMO/TOP_SRC/par_trc.F90

    r560 r732  
    1111   !!---------------------------------------------------------------------- 
    1212   !!  TOP 1.0,  LOCEAN-IPSL (2005) 
    13    !! $Header$ 
     13   !! $Header: /home/opalod/NEMOCVSROOT/NEMO/TOP_SRC/par_trc.F90,v 1.6 2007/10/12 09:22:19 opalod Exp $ 
    1414   !! This software is governed by the CeCILL licence see modipsl/doc/NEMO_CeCILL.txt 
    1515   !!---------------------------------------------------------------------- 
     
    2929!! we keep 3 more trends for eddy induced flux (gent velocity) 
    3030#    if defined key_trcdmp 
     31   INTEGER , PARAMETER :: jpdiatrc = 11 
     32#    else 
    3133   INTEGER , PARAMETER :: jpdiatrc = 10 
    32 #    else 
    33    INTEGER , PARAMETER :: jpdiatrc = 9 
    3434#    endif 
    3535#  else 
    3636#    if defined key_trcdmp 
     37   INTEGER , PARAMETER :: jpdiatrc = 8 
     38#    else 
    3739   INTEGER , PARAMETER :: jpdiatrc = 7 
    38 #    else 
    39    INTEGER , PARAMETER :: jpdiatrc = 6 
    4040#    endif 
    4141#  endif 
  • tags/nemo_v2_3/NEMO/TOP_SRC/trcdit.F90

    r560 r732  
    22   !!---------------------------------------------------------------------- 
    33   !!  TOP 1.0,  LOCEAN-IPSL (2005)  
    4    !! $Header$  
     4   !! $Header: /home/opalod/NEMOCVSROOT/NEMO/TOP_SRC/trcdit.F90,v 1.9 2007/10/12 09:22:19 opalod Exp $  
    55   !! This software is governed by the CeCILL licence see modipsl/doc/NEMO_CeCILL.txt  
    66   !!---------------------------------------------------------------------- 
     
    418418                  IF (jl.eq.7) THEN 
    419419! short and long title for x gent velocity for tracer 
    420                       WRITE (cltra,'("XGV",16a)') ctrcnm(jn) 
     420                      WRITE (cltra,'("XGV_",16a)') ctrcnm(jn) 
    421421                      WRITE (cltral,'("X gent velocity trend for ",53a)')  & 
    422422                      &      ctrcnl(jn)(1:53) 
     
    436436# endif 
    437437# if defined key_trcdmp 
    438                   IF (jl.eq.jpdiatrc)  THEN 
     438                  IF (jl.eq.jpdiatrc-1)  THEN 
    439439! last trends for tracer damping : short and long title 
    440440                      WRITE (cltra,'("TDM_",16a)') ctrcnm(jn) 
     
    443443                  END IF 
    444444# endif 
     445                  IF (jl.eq.jpdiatrc)  THEN 
     446! last trends for tracer damping : short and long title 
     447                      WRITE (cltra,'("SBC_",16a)') ctrcnm(jn) 
     448                      WRITE (cltral,'("Surface boundary flux ",58a)')  & 
     449                      &      ctrcnl(jn)(1:58) 
     450                  END IF 
     451 
    445452                  call flush(numout) 
    446453                  cltrau=ctrcun(jn) ! UNIT for tracer /trends 
     
    518525# endif 
    519526# if defined key_trcdmp 
    520                   IF (jl.eq.jpdiatrc) THEN 
     527                  IF (jl.eq.jpdiatrc-1) THEN 
    521528! short for x gent velocity for tracer 
    522529                      WRITE (cltra,'("TDM_",16a)') ctrcnm(jn) 
    523530                  END IF 
    524531# endif 
     532                  IF (jl.eq.jpdiatrc) THEN 
     533! short for surface boundary conditions for tracer 
     534                      WRITE (cltra,'("SBC_",a)') ctrcnm(jn) 
     535                  END IF 
    525536 
    526537                  CALL histwrite(nit6(jn), cltra, it, trtrd(:,:,:,ikeep(jn),jl)  & 
Note: See TracChangeset for help on using the changeset viewer.