- Timestamp:
- 2010-07-29T17:05:35+02:00 (14 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
branches/DEV_r2006_merge_TRA_TRC/NEMO/OPA_SRC/TRA/traldf.F90
r2024 r2034 45 45 # include "vectopt_loop_substitute.h90" 46 46 !!---------------------------------------------------------------------- 47 !! OPA 9.0 , LOCEAN-IPSL (2006)47 !! NEMO/OPA 3.3 , LOCEAN-IPSL (2010) 48 48 !! $Id$ 49 49 !! Software governed by the CeCILL licence (modipsl/doc/NEMO_CeCILL.txt) … … 74 74 75 75 SELECT CASE ( nldf ) ! compute lateral mixing trend and add it to the general trend 76 CASE ( 0 ) ; CALL tra_ldf_lap ( kt , 'TRA', zgtsu, zgtsv, & 77 & tsb, tsa , jpts ) ! iso-level laplacian 78 CASE ( 1 ) ; CALL tra_ldf_iso ( kt , 'TRA', zgtsu, zgtsv, & 79 & tsb, tsa , jpts , ahtb0 ) ! rotated laplacian 80 CASE ( 2 ) ; CALL tra_ldf_bilap ( kt , 'TRA', zgtsu, zgtsv, & 81 & tsb, tsa , jpts ) ! iso-level bilaplacian 82 CASE ( 3 ) ; CALL tra_ldf_bilapg( kt , 'TRA', tsb, tsa, jpts ) ! s-coord. horizontal bilaplacian 76 CASE ( 0 ) ; CALL tra_ldf_lap ( kt , 'TRA', zgtsu, zgtsv, tsb, tsa, jpts ) ! iso-level laplacian 77 CASE ( 1 ) ; CALL tra_ldf_iso ( kt , 'TRA', zgtsu, zgtsv, tsb, tsa, jpts, ahtb0 ) ! rotated laplacian 78 CASE ( 2 ) ; CALL tra_ldf_bilap ( kt , 'TRA', zgtsu, zgtsv, tsb, tsa, jpts ) ! iso-level bilaplacian 79 CASE ( 3 ) ; CALL tra_ldf_bilapg( kt , 'TRA', tsb, tsa, jpts ) ! s-coord. horizontal bilaplacian 83 80 ! 84 81 CASE ( -1 ) ! esopa: test all possibility with control print 85 CALL tra_ldf_lap ( kt , 'TRA', zgtsu, zgtsv, & 86 & tsb, tsa , jpts ) 87 CALL prt_ctl( tab3d_1=tsa(:,:,:,jp_tem), clinfo1=' ldf0 - Ta: ', mask1=tmask, & 88 & tab3d_2=tsa(:,:,:,jp_sal), clinfo2= ' Sa: ', mask2=tmask, clinfo3='tra' ) 89 CALL tra_ldf_iso ( kt , 'TRA', zgtsu, zgtsv, & 90 & tsb, tsa , jpts , ahtb0 ) 91 CALL prt_ctl( tab3d_1=tsa(:,:,:,jp_tem), clinfo1=' ldf1 - Ta: ', mask1=tmask, & 92 & tab3d_2=tsa(:,:,:,jp_sal), clinfo2= ' Sa: ', mask2=tmask, clinfo3='tra' ) 93 CALL tra_ldf_bilap ( kt , 'TRA', zgtsu, zgtsv, & 94 & tsb, tsa , jpts ) 95 CALL prt_ctl( tab3d_1=tsa(:,:,:,jp_tem), clinfo1=' ldf2 - Ta: ', mask1=tmask, & 96 & tab3d_2=tsa(:,:,:,jp_sal), clinfo2= ' Sa: ', mask2=tmask, clinfo3='tra' ) 97 CALL tra_ldf_bilapg( kt , 'TRA', tsb, tsa, jpts ) 98 CALL prt_ctl( tab3d_1=tsa(:,:,:,jp_tem), clinfo1=' ldf3 - Ta: ', mask1=tmask, & 99 & tab3d_2=tsa(:,:,:,jp_sal), clinfo2= ' Sa: ', mask2=tmask, clinfo3='tra' ) 82 CALL tra_ldf_lap ( kt , 'TRA', zgtsu, zgtsv, tsb, tsa, jpts ) 83 CALL prt_ctl( tab3d_1=tsa(:,:,:,jp_tem), clinfo1=' ldf0 - Ta: ', mask1=tmask, & 84 & tab3d_2=tsa(:,:,:,jp_sal), clinfo2= ' Sa: ', mask2=tmask, clinfo3='tra' ) 85 CALL tra_ldf_iso ( kt , 'TRA', zgtsu, zgtsv, tsb, tsa, jpts, ahtb0 ) 86 CALL prt_ctl( tab3d_1=tsa(:,:,:,jp_tem), clinfo1=' ldf1 - Ta: ', mask1=tmask, & 87 & tab3d_2=tsa(:,:,:,jp_sal), clinfo2= ' Sa: ', mask2=tmask, clinfo3='tra' ) 88 CALL tra_ldf_bilap ( kt , 'TRA', zgtsu, zgtsv, tsb, tsa, jpts ) 89 CALL prt_ctl( tab3d_1=tsa(:,:,:,jp_tem), clinfo1=' ldf2 - Ta: ', mask1=tmask, & 90 & tab3d_2=tsa(:,:,:,jp_sal), clinfo2= ' Sa: ', mask2=tmask, clinfo3='tra' ) 91 CALL tra_ldf_bilapg( kt , 'TRA', tsb, tsa, jpts ) 92 CALL prt_ctl( tab3d_1=tsa(:,:,:,jp_tem), clinfo1=' ldf3 - Ta: ', mask1=tmask, & 93 & tab3d_2=tsa(:,:,:,jp_sal), clinfo2= ' Sa: ', mask2=tmask, clinfo3='tra' ) 100 94 END SELECT 101 95
Note: See TracChangeset
for help on using the changeset viewer.