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 2082 for branches/DEV_r2006_merge_TRA_TRC/NEMO/TOP_SRC – NEMO

Ignore:
Timestamp:
2010-09-10T12:32:58+02:00 (14 years ago)
Author:
cetlod
Message:

Improve the merge of TRA-TRC, see ticket #717

Location:
branches/DEV_r2006_merge_TRA_TRC/NEMO/TOP_SRC
Files:
24 edited

Legend:

Unmodified
Added
Removed

  • branches/DEV_r2006_merge_TRA_TRC/NEMO/TOP_SRC/C14b/trcsms_c14b.F90

    r2038 r2082  
    234234            !  Computation of solubility   
    235235            IF (tmask(ji,jj,1) >  0.) THEN 
    236                ztp  = ( tn(ji,jj,1) + 273.16 ) * 0.01 
     236               ztp  = ( tsn(ji,jj,1,jp_tem) + 273.16 ) * 0.01 
    237237               zsk  = 0.023517 + ztp * ( -0.023656 + 0.0047036 * ztp ) 
    238                zsol = EXP( -60.2409 + 93.4517 / ztp  + 23.3585 * LOG( ztp ) + zsk * sn(ji,jj,1) ) 
     238               zsol = EXP( -60.2409 + 93.4517 / ztp  + 23.3585 * LOG( ztp ) + zsk * tsn(ji,jj,1,jp_sal) ) 
    239239               ! convert solubilities [mol/(l * atm)] -> [mol/(m^3 * ppm)] 
    240240               zsol = zsol * 1.0e-03 
    … …  
    247247 
    248248            ! Computes the Schmidt number of CO2 in seawater 
    249             zt   = tn(ji,jj,1) 
     249            zt   = tsn(ji,jj,1,jp_tem) 
    250250            zsch = 2073.1 + zt * ( -125.62 + zt * (3.6276 - 0.043219 * zt ) ) 
    251251 

  • branches/DEV_r2006_merge_TRA_TRC/NEMO/TOP_SRC/CFC/trcsms_cfc.F90

    r2052 r2082  
    129129               ! coefficient for solubility for CFC-11/12 in  mol/l/atm 
    130130               IF( tmask(ji,jj,1) .GE. 0.5 ) THEN 
    131                   ztap  = ( tn(ji,jj,1) + 273.16 ) * 0.01 
     131                  ztap  = ( tsn(ji,jj,1,jp_tem) + 273.16 ) * 0.01 
    132132                  zdtap = sob(1,jl) + ztap * ( sob(2,jl) + ztap * sob(3,jl) )  
    133133                  zsol  =  EXP( soa(1,jl) + soa(2,jl) / ztap + soa(3,jl) * LOG( ztap )   & 
    134                      &                    + soa(4,jl) * ztap * ztap + sn(ji,jj,1) * zdtap )  
     134                     &                    + soa(4,jl) * ztap * ztap + tsn(ji,jj,1,jp_sal) * zdtap )  
    135135               ELSE 
    136136                  zsol  = 0.e0 
    … …  
    143143               ! Computation of speed transfert 
    144144               !    Schmidt number 
    145                zt1  = tn(ji,jj,1) 
     145               zt1  = tsn(ji,jj,1,jp_tem) 
    146146               zt2  = zt1 * zt1  
    147147               zt3  = zt1 * zt2 

  • branches/DEV_r2006_merge_TRA_TRC/NEMO/TOP_SRC/LOBSTER/trcbio.F90

    r2038 r2082  
    8181      !!--------------------------------------------------------------------- 
    8282 
    83       IF( kt == nit000 ) THEN 
     83      IF( kt == nittrc000 ) THEN 
    8484         IF(lwp) WRITE(numout,*) 
    8585         IF(lwp) WRITE(numout,*) ' trc_bio: LOBSTER bio-model' 

  • branches/DEV_r2006_merge_TRA_TRC/NEMO/TOP_SRC/LOBSTER/trcexp.F90

    r2038 r2082  
    6060      !!--------------------------------------------------------------------- 
    6161 
    62       IF( kt == nit000 ) THEN 
     62      IF( kt == nittrc000 ) THEN 
    6363         IF(lwp) WRITE(numout,*) 
    6464         IF(lwp) WRITE(numout,*) ' trc_exp: LOBSTER export' 

  • branches/DEV_r2006_merge_TRA_TRC/NEMO/TOP_SRC/LOBSTER/trcopt.F90

    r2038 r2082  
    6565      !!--------------------------------------------------------------------- 
    6666 
    67       IF( kt == nit000 ) THEN 
     67      IF( kt == nittrc000 ) THEN 
    6868         IF(lwp) WRITE(numout,*) 
    6969         IF(lwp) WRITE(numout,*) ' trc_opt : LOBSTER optic-model' 

  • branches/DEV_r2006_merge_TRA_TRC/NEMO/TOP_SRC/LOBSTER/trcsed.F90

    r2038 r2082  
    6767      !!--------------------------------------------------------------------- 
    6868 
    69       IF( kt == nit000 ) THEN 
     69      IF( kt == nittrc000 ) THEN 
    7070         IF(lwp) WRITE(numout,*) 
    7171         IF(lwp) WRITE(numout,*) ' trc_sed: LOBSTER sedimentation' 

  • branches/DEV_r2006_merge_TRA_TRC/NEMO/TOP_SRC/PISCES/p4zche.F90

    r1953 r2082  
    181181 
    182182            !                             ! SET ABSOLUTE TEMPERATURE 
    183             ztkel = tn(ji,jj,1) + 273.16 
     183            ztkel = tsn(ji,jj,1,jp_tem) + 273.16 
    184184            zqtt  = ztkel * 0.01 
    185185            zqtt2 = zqtt * zqtt 
    186             zsal  = sn(ji,jj,1) + (1.- tmask(ji,jj,1) ) * 35. 
     186            zsal  = tsn(ji,jj,1,jp_sal) + (1.- tmask(ji,jj,1) ) * 35. 
    187187            zlqtt = LOG( zqtt ) 
    188188 
    … …  
    214214 
    215215               ! SET ABSOLUTE TEMPERATURE 
    216                ztkel   = tn(ji,jj,jk) + 273.16 
     216               ztkel   = tsn(ji,jj,jk,jp_tem) + 273.16 
    217217               zqtt    = ztkel * 0.01 
    218                zsal    = sn(ji,jj,jk) + ( 1.-tmask(ji,jj,jk) ) * 35. 
     218               zsal    = tsn(ji,jj,jk,jp_sal) + ( 1.-tmask(ji,jj,jk) ) * 35. 
    219219               zsqrt  = SQRT( zsal ) 
    220220               zsal15  = zsqrt * zsal 
    … …  
    224224               zis2   = zis * zis 
    225225               zisqrt = SQRT( zis ) 
    226                ztc     = tn(ji,jj,jk) + ( 1.- tmask(ji,jj,jk) ) * 20. 
     226               ztc     = tsn(ji,jj,jk,jp_tem) + ( 1.- tmask(ji,jj,jk) ) * 20. 
    227227 
    228228               ! CHLORINITY (WOOSTER ET AL., 1969) 

  • branches/DEV_r2006_merge_TRA_TRC/NEMO/TOP_SRC/PISCES/p4zflx.F90

    r2038 r2082  
    128128!CDIR NOVERRCHK 
    129129         DO ji = 1, jpi 
    130             ztc  = MIN( 35., tn(ji,jj,1) ) 
     130            ztc  = MIN( 35., tsn(ji,jj,1,jp_tem) ) 
    131131            ztc2 = ztc * ztc 
    132132            ztc3 = ztc * ztc2  

  • branches/DEV_r2006_merge_TRA_TRC/NEMO/TOP_SRC/PISCES/p4zint.F90

    r1753 r2082  
    5555      ! ------------------------------------------- 
    5656 
    57       tgfunc (:,:,:) = EXP( 0.063913 * tn(:,:,:) ) 
    58       tgfunc2(:,:,:) = EXP( 0.07608  * tn(:,:,:) ) 
     57      tgfunc (:,:,:) = EXP( 0.063913 * tsn(:,:,:,jp_tem) ) 
     58      tgfunc2(:,:,:) = EXP( 0.07608  * tsn(:,:,:,jp_tem) ) 
    5959 
    6060      ! Computation of the silicon dependant half saturation 

  • branches/DEV_r2006_merge_TRA_TRC/NEMO/TOP_SRC/PISCES/p4zlim.F90

    r1953 r2082  
    161161         DO jj = 1, jpj 
    162162            DO ji = 1, jpi 
    163                ztemp = MAX( 0., tn(ji,jj,jk) ) 
     163               ztemp = MAX( 0., tsn(ji,jj,jk,jp_tem) ) 
    164164               xfracal(ji,jj,jk) = caco3r * xlimphy(ji,jj,jk)   & 
    165165                  &                       * MAX( 0.0001, ztemp / ( 2.+ ztemp ) )   & 

  • branches/DEV_r2006_merge_TRA_TRC/NEMO/TOP_SRC/PISCES/p4zprod.F90

    r2038 r2082  
    147147               ! Computation of the P-I slope for nanos and diatoms 
    148148               IF( etot(ji,jj,jk) > 1.E-3 ) THEN 
    149                    ztn    = MAX( 0., tn(ji,jj,jk) - 15. ) 
     149                   ztn    = MAX( 0., tsn(ji,jj,jk,jp_tem) - 15. ) 
    150150                   zadap  = 0.+ 1.* ztn / ( 2.+ ztn ) 
    151151                   zadap2 = 0.e0 

  • branches/DEV_r2006_merge_TRA_TRC/NEMO/TOP_SRC/PISCES/p4zrem.F90

    r2038 r2082  
    287287               zsatur  = ( sio3eq(ji,jj,jk) - trn(ji,jj,jk,jpsil) ) / ( sio3eq(ji,jj,jk) + rtrn ) 
    288288               zsatur  = MAX( rtrn, zsatur ) 
    289                zsatur2 = zsatur * ( 1. + tn(ji,jj,jk) / 400.)**4 
    290                znusil  = 0.225  * ( 1. + tn(ji,jj,jk) / 15.) * zsatur + 0.775 * zsatur2**9 
     289               zsatur2 = zsatur * ( 1. + tsn(ji,jj,jk,jp_tem) / 400.)**4 
     290               znusil  = 0.225  * ( 1. + tsn(ji,jj,jk,jp_tem) / 15.) * zsatur + 0.775 * zsatur2**9 
    291291#    if defined key_degrad 
    292292               zsiremin = xsirem * xstep * znusil * facvol(ji,jj,jk) 

  • branches/DEV_r2006_merge_TRA_TRC/NEMO/TOP_SRC/PISCES/p4zsed.F90

    r2038 r2082  
    383383      imois2 = nmonth 
    384384 
    385       ! 1. first call kt=nit000 
     385      ! 1. first call kt=nittrc000 
    386386      ! ----------------------- 
    387387 
    388       IF( kt == nit000 ) THEN 
     388      IF( kt == nittrc000 ) THEN 
    389389         ! initializations 
    390390         nflx1  = 0 
    … …  
    402402      ! ---------------- 
    403403 
    404       IF( kt == nit000 .OR. imois /= nflx1 ) THEN 
     404      IF( kt == nittrc000 .OR. imois /= nflx1 ) THEN 
    405405 
    406406         ! Calendar computation 
    … …  
    423423         CALL iom_get ( numdust, jpdom_data, 'dust', dustmo(:,:,1), nflx1 ) 
    424424         CALL iom_get ( numdust, jpdom_data, 'dust', dustmo(:,:,2), nflx2 ) 
    425  
    426          IF(lwp .AND. nitend-nit000 <= 100 ) THEN 
    427             WRITE(numout,*) 
    428             WRITE(numout,*) ' read clio flx ok' 
    429             WRITE(numout,*) 
    430                WRITE(numout,*) 
    431                WRITE(numout,*) 'Clio month: ',nflx1,'  field: dust' 
    432                CALL prihre( dustmo(:,:,1),jpi,jpj,1,jpi,20,1,jpj,10,1e9,numout ) 
    433          ENDIF 
    434425 
    435426      ENDIF 

  • branches/DEV_r2006_merge_TRA_TRC/NEMO/TOP_SRC/SED/par_sed.F90

    r1250 r2082  
    1616      jpjm1    =>   jpjm1 ,  & !: jpj - 1 
    1717      jpij     =>   jpij       !: jpi x jpj 
     18      jp_tem   =>   jp_tem     !: indice of temperature 
     19      jp_sal   =>   jp_sal     !: indice of salintity 
    1820 
    1921#if ! defined key_sed_off 

  • branches/DEV_r2006_merge_TRA_TRC/NEMO/TOP_SRC/SED/sed.F90

    r1715 r2082  
    3434 
    3535   USE oce , ONLY :            & 
    36       tn      =>    tn    ,  & !: pot. temperature (celsius) 
    37       sn      =>    sn         !: salinity (psu) 
     36      tsn      =>    tsn       & !: pot. temperature (celsius) and salinity (psu) 
    3837 
    3938   USE trc, ONLY :  & 
    … …  
    215214   INTEGER, PUBLIC :: & 
    216215     numsed = 27 
    217     
    218216#else 
    219217   !!====================================================================== 

  • branches/DEV_r2006_merge_TRA_TRC/NEMO/TOP_SRC/SED/seddta.F90

    r1264 r2082  
    131131                  trc_data(ji,jj,9 ) = sinking2(ji,jj,ikt) 
    132132                  trc_data(ji,jj,10) = sinkcal (ji,jj,ikt) 
    133                   trc_data(ji,jj,11) = tn      (ji,jj,ikt) 
    134                   trc_data(ji,jj,12) = sn      (ji,jj,ikt) 
     133                  trc_data(ji,jj,11) = tsn     (ji,jj,ikt,jp_tem) 
     134                  trc_data(ji,jj,12) = tsn     (ji,jj,ikt,jp_sal) 
    135135#   else 
    136136                  trc_data(ji,jj,7 ) = sinksil (ji,jj,ikt) 
    137137                  trc_data(ji,jj,8 ) = sinking (ji,jj,ikt) 
    138138                  trc_data(ji,jj,9 ) = sinkcal (ji,jj,ikt) 
    139                   trc_data(ji,jj,10) = tn      (ji,jj,ikt) 
    140                   trc_data(ji,jj,11) = sn      (ji,jj,ikt)        
     139                  trc_data(ji,jj,10) = tsn     (ji,jj,ikt,jp_tem) 
     140                  trc_data(ji,jj,11) = tsn     (ji,jj,ikt,jp_sal)        
    141141#   endif 
    142142               ENDIF 

  • branches/DEV_r2006_merge_TRA_TRC/NEMO/TOP_SRC/TRP/trcadv.F90

    r2034 r2082  
    3333 
    3434   PUBLIC   trc_adv    ! routine called by step module 
     35 
    3536   INTEGER ::   nadv   ! choice of the type of advection scheme 
     37   REAL(wp), DIMENSION(jpk) ::   r2dt  ! vertical profile time-step, = 2 rdttra 
     38      !                                ! except at nit000 (=rdttra) if neuler=0 
    3639 
    3740   !! * Substitutions 
    … …  
    6467      IF( kt == nittrc000 )   CALL trc_adv_ctl          ! initialisation & control of options 
    6568 
     69      IF( neuler == 0 .AND. kt == nittrc000 ) THEN     ! at nit000 
     70         r2dt(:) =  rdttra(:) * FLOAT(nn_dttrc)          ! = rdtra (restarting with Euler time stepping) 
     71      ELSEIF( kt <= nittrc000 + nn_dttrc ) THEN          ! at nit000 or nit000+1 
     72         r2dt(:) = 2. * rdttra(:) * FLOAT(nn_dttrc)      ! = 2 rdttra (leapfrog) 
     73      ENDIF 
     74 
    6675      !                                                   ! effective transport 
    6776      DO jk = 1, jpkm1 
    … …  
    7887      ! 
    7988      SELECT CASE ( nadv )                            !==  compute advection trend and add it to general trend  ==! 
    80       CASE ( 1 )   ;    CALL tra_adv_cen2  ( kt, 'TRC', zun, zvn, zwn, trb, trn, tra, jptra )   !  2nd order centered 
    81       CASE ( 2 )   ;    CALL tra_adv_tvd   ( kt, 'TRC', zun, zvn, zwn, trb, trn, tra, jptra )   !  TVD  
    82       CASE ( 3 )   ;    CALL tra_adv_muscl ( kt, 'TRC', zun, zvn, zwn, trb,      tra, jptra )   !  MUSCL  
    83       CASE ( 4 )   ;    CALL tra_adv_muscl2( kt, 'TRC', zun, zvn, zwn, trb, trn, tra, jptra )   !  MUSCL2  
    84       CASE ( 5 )   ;    CALL tra_adv_ubs   ( kt, 'TRC', zun, zvn, zwn, trb, trn, tra, jptra )   !  UBS  
    85       CASE ( 6 )   ;    CALL tra_adv_qck   ( kt, 'TRC', zun, zvn, zwn, trb, trn, tra, jptra )   !  QUICKEST  
     89      CASE ( 1 )   ;    CALL tra_adv_cen2  ( kt, 'TRC',       zun, zvn, zwn, trb, trn, tra, jptra )   !  2nd order centered 
     90      CASE ( 2 )   ;    CALL tra_adv_tvd   ( kt, 'TRC', r2dt, zun, zvn, zwn, trb, trn, tra, jptra )   !  TVD  
     91      CASE ( 3 )   ;    CALL tra_adv_muscl ( kt, 'TRC', r2dt, zun, zvn, zwn, trb,      tra, jptra )   !  MUSCL  
     92      CASE ( 4 )   ;    CALL tra_adv_muscl2( kt, 'TRC', r2dt, zun, zvn, zwn, trb, trn, tra, jptra )   !  MUSCL2  
     93      CASE ( 5 )   ;    CALL tra_adv_ubs   ( kt, 'TRC', r2dt, zun, zvn, zwn, trb, trn, tra, jptra )   !  UBS  
     94      CASE ( 6 )   ;    CALL tra_adv_qck   ( kt, 'TRC', r2dt, zun, zvn, zwn, trb, trn, tra, jptra )   !  QUICKEST  
    8695      ! 
    8796      CASE (-1 )                                      !==  esopa: test all possibility with control print  ==! 
    88          CALL tra_adv_cen2  ( kt, 'TRC', zun, zvn, zwn, trb, trn, tra, jptra )           
     97         CALL tra_adv_cen2  ( kt, 'TRC',       zun, zvn, zwn, trb, trn, tra, jptra )           
    8998         WRITE(charout, FMT="('adv1')")  ; CALL prt_ctl_trc_info(charout) 
    9099                                           CALL prt_ctl_trc(tab4d=tra, mask=tmask, clinfo=ctrcnm,clinfo2='trd') 
    91          CALL tra_adv_tvd   ( kt, 'TRC', zun, zvn, zwn, trb, trn, tra, jptra )           
     100         CALL tra_adv_tvd   ( kt, 'TRC', r2dt, zun, zvn, zwn, trb, trn, tra, jptra )           
    92101         WRITE(charout, FMT="('adv2')")  ; CALL prt_ctl_trc_info(charout) 
    93102                                           CALL prt_ctl_trc(tab4d=tra, mask=tmask, clinfo=ctrcnm,clinfo2='trd') 
    94          CALL tra_adv_muscl ( kt, 'TRC', zun, zvn, zwn, trb,      tra, jptra )           
     103         CALL tra_adv_muscl ( kt, 'TRC', r2dt, zun, zvn, zwn, trb,      tra, jptra )           
    95104         WRITE(charout, FMT="('adv3')")  ; CALL prt_ctl_trc_info(charout) 
    96105                                           CALL prt_ctl_trc(tab4d=tra, mask=tmask, clinfo=ctrcnm,clinfo2='trd') 
    97          CALL tra_adv_muscl2( kt, 'TRC', zun, zvn, zwn, trb, trn, tra, jptra )           
     106         CALL tra_adv_muscl2( kt, 'TRC', r2dt, zun, zvn, zwn, trb, trn, tra, jptra )           
    98107         WRITE(charout, FMT="('adv4')")  ; CALL prt_ctl_trc_info(charout) 
    99108                                           CALL prt_ctl_trc(tab4d=tra, mask=tmask, clinfo=ctrcnm,clinfo2='trd') 
    100          CALL tra_adv_ubs   ( kt, 'TRC', zun, zvn, zwn, trb, trn, tra, jptra )           
     109         CALL tra_adv_ubs   ( kt, 'TRC', r2dt, zun, zvn, zwn, trb, trn, tra, jptra )           
    101110         WRITE(charout, FMT="('adv5')")  ; CALL prt_ctl_trc_info(charout) 
    102111                                           CALL prt_ctl_trc(tab4d=tra, mask=tmask, clinfo=ctrcnm,clinfo2='trd') 
    103          CALL tra_adv_qck   ( kt, 'TRC', zun, zvn, zwn, trb, trn, tra, jptra )           
     112         CALL tra_adv_qck   ( kt, 'TRC', r2dt, zun, zvn, zwn, trb, trn, tra, jptra )           
    104113         WRITE(charout, FMT="('adv6')")  ; CALL prt_ctl_trc_info(charout) 
    105114                                           CALL prt_ctl_trc(tab4d=tra, mask=tmask, clinfo=ctrcnm,clinfo2='trd') 

  • branches/DEV_r2006_merge_TRA_TRC/NEMO/TOP_SRC/TRP/trcldf.F90

    r2034 r2082  
    7070 
    7171      SELECT CASE ( nldf )                       ! compute lateral mixing trend and add it to the general trend 
    72       CASE ( 0 )   ;   CALL tra_ldf_lap   ( kt , 'TRC', gtru, gtrv, trb, tra, jptra            )  ! iso-level laplacian 
    73       CASE ( 1 )   ;   CALL tra_ldf_iso   ( kt , 'TRC', gtru, gtrv, trb, tra, jptra, rn_ahtb_0 )  ! rotated laplacian  
    74       CASE ( 2 )   ;   CALL tra_ldf_bilap ( kt , 'TRC', gtru, gtrv, trb, tra, jptra            )  ! iso-level bilaplacian 
    75       CASE ( 3 )   ;   CALL tra_ldf_bilapg( kt , 'TRC',             trb, tra, jptra            )  ! s-coord. horizontal bilaplacian 
     72      CASE ( 0 )   ;   CALL tra_ldf_lap   ( kt, 'TRC', gtru, gtrv, trb, tra, jptra            )  ! iso-level laplacian 
     73      CASE ( 1 )   ;   CALL tra_ldf_iso   ( kt, 'TRC', gtru, gtrv, trb, tra, jptra, rn_ahtb_0 )  ! rotated laplacian  
     74      CASE ( 2 )   ;   CALL tra_ldf_bilap ( kt, 'TRC', gtru, gtrv, trb, tra, jptra            )  ! iso-level bilaplacian 
     75      CASE ( 3 )   ;   CALL tra_ldf_bilapg( kt, 'TRC',             trb, tra, jptra            )  ! s-coord. horizontal bilaplacian 
    7676         ! 
    7777      CASE ( -1 )                                     ! esopa: test all possibility with control print 
    78          CALL tra_ldf_lap   ( kt , 'TRC', gtru, gtrv, trb, tra, jptra            ) 
     78         CALL tra_ldf_lap   ( kt, 'TRC', gtru, gtrv, trb, tra, jptra            ) 
    7979         WRITE(charout, FMT="('ldf0 ')") ;  CALL prt_ctl_trc_info(charout) 
    8080                                            CALL prt_ctl_trc( tab4d=tra, mask=tmask, clinfo=ctrcnm, clinfo2='trd' ) 
    81          CALL tra_ldf_iso   ( kt , 'TRC', gtru, gtrv, trb, tra, jptra, rn_ahtb_0 ) 
     81         CALL tra_ldf_iso   ( kt, 'TRC', gtru, gtrv, trb, tra, jptra, rn_ahtb_0 ) 
    8282         WRITE(charout, FMT="('ldf1 ')") ;  CALL prt_ctl_trc_info(charout) 
    8383                                            CALL prt_ctl_trc( tab4d=tra, mask=tmask, clinfo=ctrcnm, clinfo2='trd' ) 
    84          CALL tra_ldf_bilap ( kt , 'TRC', gtru, gtrv, trb, tra, jptra            ) 
     84         CALL tra_ldf_bilap ( kt, 'TRC', gtru, gtrv, trb, tra, jptra            ) 
    8585         WRITE(charout, FMT="('ldf2 ')") ;  CALL prt_ctl_trc_info(charout) 
    8686                                            CALL prt_ctl_trc( tab4d=tra, mask=tmask, clinfo=ctrcnm, clinfo2='trd' ) 
    87          CALL tra_ldf_bilapg( kt , 'TRC',             trb, tra, jptra            ) 
     87         CALL tra_ldf_bilapg( kt, 'TRC',             trb, tra, jptra            ) 
    8888         WRITE(charout, FMT="('ldf3 ')") ;  CALL prt_ctl_trc_info(charout) 
    8989                                            CALL prt_ctl_trc( tab4d=tra, mask=tmask, clinfo=ctrcnm, clinfo2='trd' ) 

  • branches/DEV_r2006_merge_TRA_TRC/NEMO/TOP_SRC/TRP/trcnxt.F90

    r2034 r2082  
    4444   PUBLIC trc_nxt          ! routine called by step.F90 
    4545 
    46   REAL(wp), DIMENSION(jpk) ::   r2dt_t 
     46  REAL(wp), DIMENSION(jpk) ::   r2dt 
    4747   !!---------------------------------------------------------------------- 
    4848   !!   TOP 1.0 , LOCEAN-IPSL (2005)  
    … …  
    109109 
    110110      ! set time step size (Euler/Leapfrog) 
    111       IF( neuler == 0 .AND. kt ==  nittrc000) THEN  ;  r2dt_t(:) =     rdttra(:) * FLOAT( nn_dttrc )  ! at nit000             (Euler) 
    112       ELSEIF( kt <= nittrc000 + 1 )           THEN  ;  r2dt_t(:) = 2.* rdttra(:) * FLOAT( nn_dttrc )  ! at nit000 or nit000+1 (Leapfrog) 
     111      IF( neuler == 0 .AND. kt ==  nittrc000) THEN  ;  r2dt(:) =     rdttra(:) * FLOAT( nn_dttrc )  ! at nit000             (Euler) 
     112      ELSEIF( kt <= nittrc000 + 1 )           THEN  ;  r2dt(:) = 2.* rdttra(:) * FLOAT( nn_dttrc )  ! at nit000 or nit000+1 (Leapfrog) 
    113113      ENDIF 
    114114 
    … …  
    133133         DO jn = 1, jptra 
    134134            DO jk = 1, jpkm1 
    135                zfact = 1.e0 / r2dt_t(jk)   
     135               zfact = 1.e0 / r2dt(jk)   
    136136               ztrdt(:,:,jk,jn) = ( trb(:,:,jk,jn) - ztrdt(:,:,jk,jn) ) * zfact  
    137137               CALL trd_tra( kt, 'TRC', jn, jptra_trd_atf, ztrdt ) 

  • branches/DEV_r2006_merge_TRA_TRC/NEMO/TOP_SRC/TRP/trctrp.F90

    r2052 r2082  
    2828   USE trcrad          ! positivity                          (trc_rad routine) 
    2929   USE trcsbc          ! surface boundary condition          (trc_sbc routine) 
    30    USE zpshde          ! partial step: hor. derivative   (zps_hde_trc routine) 
     30   USE zpshde          ! partial step: hor. derivative       (zps_hde routine) 
    3131 
    3232#if defined key_agrif 
    … …  
    7474#endif 
    7575                                CALL trc_zdf( kstp )            ! vertical mixing and after tracer fields 
    76          IF( ln_zps    )        CALL zps_hde_trc( kstp, jptra, trn, gtru, gtrv )  ! Partial steps: now horizontal gradient of passive 
     76         IF( ln_zps    )        CALL zps_hde( kstp, jptra, trn, gtru, gtrv )  ! Partial steps: now horizontal gradient of passive 
    7777                                                                ! tracers at the bottom ocean level 
    7878                                CALL trc_nxt( kstp )            ! tracer fields at next time step      

  • branches/DEV_r2006_merge_TRA_TRC/NEMO/TOP_SRC/TRP/trczdf.F90

    r2052 r2082  
    8080      SELECT CASE ( nzdf )                       ! compute lateral mixing trend and add it to the general trend 
    8181      CASE ( -1 )                                       ! esopa: test all possibility with control print 
    82          CALL tra_zdf_exp( kt , 'TRC', r2dt, nn_trczdf_exp, trb, tra, jptra )  
     82         CALL tra_zdf_exp( kt, 'TRC', r2dt, nn_trczdf_exp, trb, tra, jptra )  
    8383         WRITE(charout, FMT="('zdf1 ')") ;  CALL prt_ctl_trc_info(charout) 
    8484                                            CALL prt_ctl_trc( tab4d=tra, mask=tmask, clinfo=ctrcnm, clinfo2='trd' ) 
    85          CALL tra_zdf_imp( kt , 'TRC', r2dt,                trb, tra, jptra )  
     85         CALL tra_zdf_imp( kt, 'TRC', r2dt,                trb, tra, jptra )  
    8686         WRITE(charout, FMT="('zdf2 ')") ;  CALL prt_ctl_trc_info(charout) 
    8787                                            CALL prt_ctl_trc( tab4d=tra, mask=tmask, clinfo=ctrcnm, clinfo2='trd' ) 
    88       CASE ( 0 ) ;  CALL tra_zdf_exp( kt , 'TRC', r2dt, nn_trczdf_exp, trb, tra, jptra )    !   explicit scheme  
    89       CASE ( 1 ) ;  CALL tra_zdf_imp( kt , 'TRC', r2dt,                trb, tra, jptra )    !   implicit scheme           
     88      CASE ( 0 ) ;  CALL tra_zdf_exp( kt, 'TRC', r2dt, nn_trczdf_exp, trb, tra, jptra )    !   explicit scheme  
     89      CASE ( 1 ) ;  CALL tra_zdf_imp( kt, 'TRC', r2dt,                trb, tra, jptra )    !   implicit scheme           
    9090 
    9191      END SELECT 

  • branches/DEV_r2006_merge_TRA_TRC/NEMO/TOP_SRC/oce_trc.F90

    r2052 r2082  
    173173   USE oce , ONLY :   l_traldf_rot => l_traldf_rot  !: rotated laplacian operator for lateral diffusion 
    174174#if defined key_offline 
    175    USE oce , ONLY :   gtu     =>    gtu     !: t-, s- and rd horizontal gradient at u- and 
    176    USE oce , ONLY :   gsu     =>    gsu     !: v-points at bottom ocean level 
     175   USE oce , ONLY :   gtsu    =>    gtsu    !: t-, s- and rd horizontal gradient at u- and 
     176   USE oce , ONLY :   gtsv    =>    gtsv    !: 
    177177   USE oce , ONLY :   gru     =>    gru     !: 
    178    USE oce , ONLY :   gtv     =>    gtv     !: 
    179    USE oce , ONLY :   gsv     =>    gsv     !: 
    180178   USE oce , ONLY :   grv     =>    grv     !:  
    181179# if defined key_degrad 

  • branches/DEV_r2006_merge_TRA_TRC/NEMO/TOP_SRC/trc.F90

    r2052 r2082  
    5555   !! passive tracers restart (input and output) 
    5656   !! ------------------------------------------   
    57    INTEGER , PUBLIC  ::  nn_dttrc      !: frequency of step on passive tracers 
    58    INTEGER , PUBLIC  ::  nittrc000     !: first time step of passive tracers model 
    5957   LOGICAL , PUBLIC  ::  ln_rsttr      !: boolean term for restart i/o for passive tracers (namelist) 
    6058   LOGICAL , PUBLIC  ::  lrst_trc      !: logical to control the trc restart write 

  • branches/DEV_r2006_merge_TRA_TRC/NEMO/TOP_SRC/trcini.F90

    r2038 r2082  
    2929   USE daymod 
    3030#endif 
    31    USE zpshde          ! partial step: hor. derivative   (zps_hde_trc routine) 
     31   USE zpshde          ! partial step: hor. derivative   (zps_hde routine) 
    3232   USE in_out_manager  ! I/O manager 
    3333   USE prtctl_trc      ! Print control passive tracers (prt_ctl_trc_init routine) 
    … …  
    4949CONTAINS 
    5050    
    51    SUBROUTINE trc_ini 
     51   SUBROUTINE trc_init 
    5252      !!--------------------------------------------------------------------- 
    53       !!                     ***  ROUTINE trc_ini  *** 
     53      !!                     ***  ROUTINE trc_init  *** 
    5454      !! 
    5555      !! ** Purpose :   Initialization of the passive tracer fields  
    … …  
    6767 
    6868      IF(lwp) WRITE(numout,*) 
    69       IF(lwp) WRITE(numout,*) 'trc_ini : initial set up of the passive tracers' 
     69      IF(lwp) WRITE(numout,*) 'trc_init : initial set up of the passive tracers' 
    7070      IF(lwp) WRITE(numout,*) '~~~~~~~' 
    7171 
    … …  
    138138       
    139139      IF( ln_zps .AND. .NOT. lk_trc_c1d )   &              ! Partial steps: before horizontal gradient of passive 
    140       &                     CALL zps_hde_trc( nittrc000, jptra, trb, gtru, gtrv )       ! tracers at the bottom ocean level 
     140      &                     CALL zps_hde( nittrc000, jptra, trb, gtru, gtrv )       ! tracers at the bottom ocean level 
    141141 
    142142 
    … …  
    170170      ENDIF 
    171171 
    172    END SUBROUTINE trc_ini 
     172   END SUBROUTINE trc_init 
    173173 
    174174#else 
    … …  
    177177   !!---------------------------------------------------------------------- 
    178178CONTAINS 
    179    SUBROUTINE trc_ini                      ! Dummy routine    
    180    END SUBROUTINE trc_ini 
     179   SUBROUTINE trc_init                      ! Dummy routine    
     180   END SUBROUTINE trc_init 
    181181#endif 
    182182 
Note: See TracChangeset for help on using the changeset viewer.