Changeset 789 for trunk/NEMO/OPA_SRC/DYN/dynspg.F90
- Timestamp:
- 2008-01-11T19:04:56+01:00 (16 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
trunk/NEMO/OPA_SRC/DYN/dynspg.F90
r719 r789 20 20 USE dynspg_flt ! surface pressure gradient (dyn_spg_flt routine) 21 21 USE dynspg_rl ! surface pressure gradient (dyn_spg_rl routine) 22 USE dynspg_exp_jki ! surface pressure gradient (dyn_spg_exp_jki routine)23 USE dynspg_ts_jki ! surface pressure gradient (dyn_spg_ts_jki routine)24 USE dynspg_flt_jki ! surface pressure gradient (dyn_spg_flt_jki routine)25 22 USE trdmod ! ocean dynamics trends 26 23 USE trdmod_oce ! ocean variables trends … … 68 65 69 66 SELECT CASE ( nspg ) ! compute surf. pressure gradient trend and add it to the general trend 70 ! ! k-j-i loops67 ! 71 68 CASE ( 0 ) ; CALL dyn_spg_exp ( kt ) ! explicit 72 69 CASE ( 1 ) ; CALL dyn_spg_ts ( kt ) ! time-splitting 73 70 CASE ( 2 ) ; CALL dyn_spg_flt ( kt, kindic ) ! filtered 74 71 CASE ( 3 ) ; CALL dyn_spg_rl ( kt, kindic ) ! rigid lid 75 ! ! j-k-i loops 76 CASE ( 10 ) ; CALL dyn_spg_exp_jki( kt ) ! explicit with j-k-i loop 77 CASE ( 11 ) ; CALL dyn_spg_ts_jki ( kt ) ! time-splitting with j-k-i loop 78 CASE ( 12 ) ; CALL dyn_spg_flt_jki( kt, kindic ) ! filtered with j-k-i loop 79 ! 72 ! 80 73 CASE ( -1 ) ! esopa: test all possibility with control print 81 74 ; CALL dyn_spg_exp ( kt ) … … 87 80 ; CALL dyn_spg_flt ( kt, kindic ) 88 81 ; CALL prt_ctl( tab3d_1=ua, clinfo1=' spg2 - Ua: ', mask1=umask, & 89 & tab3d_2=va, clinfo2= ' Va: ', mask2=vmask, clinfo3='dyn' )90 ; CALL dyn_spg_exp_jki( kt )91 ; CALL prt_ctl( tab3d_1=ua, clinfo1=' spg10- Ua: ', mask1=umask, &92 & tab3d_2=va, clinfo2= ' Va: ', mask2=vmask, clinfo3='dyn' )93 ; CALL dyn_spg_ts_jki ( kt )94 ; CALL prt_ctl( tab3d_1=ua, clinfo1=' spg12- Ua: ', mask1=umask, &95 & tab3d_2=va, clinfo2= ' Va: ', mask2=vmask, clinfo3='dyn' )96 ; CALL dyn_spg_flt_jki( kt, kindic )97 ; CALL prt_ctl( tab3d_1=ua, clinfo1=' spg13- Ua: ', mask1=umask, &98 82 & tab3d_2=va, clinfo2= ' Va: ', mask2=vmask, clinfo3='dyn' ) 99 83 END SELECT … … 159 143 IF( lk_dynspg_flt) nspg = 2 160 144 IF( lk_dynspg_rl ) nspg = 3 161 IF( lk_jki ) nspg = nspg + 10162 145 IF( nspg == 13 ) nspg = 3 163 146 … … 171 154 IF( nspg == 2 ) WRITE(numout,*) ' filtered free surface' 172 155 IF( nspg == 3 ) WRITE(numout,*) ' rigid-lid' 173 IF( nspg == 10 ) WRITE(numout,*) ' explicit free surface with j-k-i loop'174 IF( nspg == 11 ) WRITE(numout,*) ' time splitting free surface with j-k-i loop'175 IF( nspg == 12 ) WRITE(numout,*) ' filtered free surface with j-k-i loop'176 156 ENDIF 177 157
Note: See TracChangeset
for help on using the changeset viewer.