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 786 for branches/dev_001_GM/NEMO/OPA_SRC/TRA/traldf.F90 – NEMO

Ignore:
Timestamp:
2008-01-10T18:11:23+01:00 (17 years ago)
Author:
gm
Message:

dev_001_GM - merge TRC-TRA on OPA only, trabbl & zpshde not done and trdmld not OK - compilation OK

File:
1 edited

Legend:

Unmodified
Added
Removed

  • branches/dev_001_GM/NEMO/OPA_SRC/TRA/traldf.F90

    r719 r786  
    7070 
    7171      SELECT CASE ( nldf )                       ! compute lateral mixing trend and add it to the general trend 
    72       CASE ( 0 )   ;   CALL tra_ldf_lap   ( kt )      ! iso-level laplacian 
    73       CASE ( 1 )   ;   CALL tra_ldf_iso   ( kt )      ! rotated laplacian (except dk[ dk[.] ] part) 
    74       CASE ( 2 )   ;   CALL tra_ldf_bilap ( kt )      ! iso-level bilaplacian 
    75       CASE ( 3 )   ;   CALL tra_ldf_bilapg( kt )      ! s-coord. horizontal bilaplacian 
     72      CASE ( 0 )   ;   CALL tra_ldf_lap   ( kt, 'TRA', jp_tem, gtu, gtv, tb, ta )      ! iso-level laplacian 
     73                       CALL tra_ldf_lap   ( kt, 'TRA', jp_sal, gsu, gsv, sb, sa )      ! iso-level laplacian 
     74      CASE ( 1 )   ;   CALL tra_ldf_iso   ( kt, 'TRA', jp_tem, gtu, gtv, tb, ta )      ! rotated laplacian except dk2 
     75      CASE ( 2 )   ;   CALL tra_ldf_bilap ( kt, 'TRA', jp_tem, gtu, gtv, tb, ta )      ! rotated laplacian except dk2 
     76                       CALL tra_ldf_bilap ( kt, 'TRA', jp_sal, gsu, gsv, sb, sa )      ! iso-level laplacian 
     77      CASE ( 3 )   ;   CALL tra_ldf_bilapg( kt, 'TRA', jp_tem, gtu, gtv, tb, ta )      ! s-coord. horizontal bilaplacian 
     78                       CALL tra_ldf_bilapg( kt, 'TRA', jp_sal, gsu, gsv, sb, sa )      ! s-coord. horizontal bilaplacian 
    7679         ! 
    7780      CASE ( -1 )                                     ! esopa: test all possibility with control print 
    78          CALL tra_ldf_lap    ( kt ) 
    79          CALL prt_ctl( tab3d_1=ta, clinfo1=' ldf0 - Ta: ', mask1=tmask,               & 
    80             &          tab3d_2=sa, clinfo2=       ' Sa: ', mask2=tmask, clinfo3='tra' ) 
    81          CALL tra_ldf_iso    ( kt ) 
    82          CALL prt_ctl( tab3d_1=ta, clinfo1=' ldf1 - Ta: ', mask1=tmask,               & 
    83             &          tab3d_2=sa, clinfo2=       ' Sa: ', mask2=tmask, clinfo3='tra' ) 
    84          CALL tra_ldf_bilap  ( kt ) 
    85          CALL prt_ctl( tab3d_1=ta, clinfo1=' ldf2 - Ta: ', mask1=tmask,               & 
    86             &          tab3d_2=sa, clinfo2=       ' Sa: ', mask2=tmask, clinfo3='tra' ) 
    87          CALL tra_ldf_bilapg ( kt ) 
    88          CALL prt_ctl( tab3d_1=ta, clinfo1=' ldf3 - Ta: ', mask1=tmask,               & 
    89             &          tab3d_2=sa, clinfo2=       ' Sa: ', mask2=tmask, clinfo3='tra' ) 
     81                       CALL tra_ldf_lap   ( kt, 'TRA', jp_tem, gtu, gtv, tb, ta )      ! iso-level laplacian 
     82                       CALL tra_ldf_lap   ( kt, 'TRA', jp_sal, gsu, gsv, sb, sa )      ! iso-level laplacian 
     83                       CALL tra_ldf_iso   ( kt, 'TRA', jp_tem, gtu, gtv, tb, ta )      ! rotated laplacian except dk2 
     84                       CALL tra_ldf_iso   ( kt, 'TRA', jp_sal, gsu, gsv, sb, sa )      ! rotated laplacian except dk2 
     85                       CALL tra_ldf_bilap ( kt, 'TRA', jp_tem, gtu, gtv, tb, ta )      ! iso-level bilaplacian 
     86                       CALL tra_ldf_bilap ( kt, 'TRA', jp_sal, gsu, gsv, sb, sa )      ! iso-level laplacian 
     87                       CALL tra_ldf_bilapg( kt, 'TRA', jp_tem, gtu, gtv, tb, ta )      ! s-coord. horizontal bilaplacian 
     88                       CALL tra_ldf_bilapg( kt, 'TRA', jp_sal, gsu, gsv, sb, sa )      ! s-coord. horizontal bilaplacian 
    9089      END SELECT 
    9190 
    … …  
    9796         ztrdt(:,:,:) = ta(:,:,:) - ztrdt(:,:,:) 
    9897         ztrds(:,:,:) = sa(:,:,:) - ztrds(:,:,:) 
    99          CALL trd_mod( ztrdt, ztrds, jptra_trd_ldf, 'TRA', kt ) 
     98         CALL trd_tra( kt, jp_tem, jpt_trd_ldf, 'TRA', ptrd3d=ztrdt) 
     99         CALL trd_tra( kt, jp_sal, jpt_trd_ldf, 'TRA', ptrd3d=ztrds) 
    100100      ENDIF 
    101101      !                                          ! print mean trends (used for debugging) 
Note: See TracChangeset for help on using the changeset viewer.