Changeset 5866 for branches/2015/dev_r5836_NOC3_vvl_by_default/NEMOGCM/NEMO/OPA_SRC/TRA/trazdf_exp.F90
- Timestamp:
- 2015-11-06T11:42:35+01:00 (9 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
branches/2015/dev_r5836_NOC3_vvl_by_default/NEMOGCM/NEMO/OPA_SRC/TRA/trazdf_exp.F90
r5845 r5866 66 66 !! (N.B. bottom condition is applied through the masked field avt). 67 67 !! - the after tracer fields due to the whole trend is 68 !! obtained in leap-frog environment by : 69 !! pta = zwx + p2dt pta 70 !! - in case of variable level thickness (lk_vvl=T) the 71 !! the leap-frog is applied on thickness weighted tracer. That is: 68 !! obtained in leap-frog environment applied on thickness weighted tracer by : 72 69 !! pta = [ ptb*e3tb + e3tn*( zwx - ptb + p2dt pta ) ] / e3tn 73 70 !! 74 71 !! ** Action : - after tracer fields pta 75 72 !!--------------------------------------------------------------------- 76 !77 73 INTEGER , INTENT(in ) :: kt ! ocean time-step index 78 74 INTEGER , INTENT(in ) :: kit000 ! first time step index … … 92 88 IF( nn_timing == 1 ) CALL timing_start('tra_zdf_exp') 93 89 ! 94 CALL wrk_alloc( jpi, jpj, jpk,zwx, zwy )90 CALL wrk_alloc( jpi,jpj,jpk, zwx, zwy ) 95 91 ! 96 92 … … 121 117 DO jj = 2, jpjm1 122 118 DO ji = fs_2, fs_jpim1 ! vector opt. 123 zave3r = 1.e0 / e3w_n(ji,jj,jk) 119 zave3r = 1.e0 / e3w_n(ji,jj,jk) !!gm here it should be e3w_b ???? 124 120 IF( cdtype == 'TRA' .AND. jn == jp_tem ) THEN ! temperature : use of avt 125 121 zwy(ji,jj,jk) = avt(ji,jj,jk) * ( zwx(ji,jj,jk-1) - zwx(ji,jj,jk) ) * zave3r … … 144 140 ! After tracer due to all trends 145 141 ! ------------------------------ 146 IF( lk_vvl ) THEN! variable level thickness : leap-frog on tracer*e3t142 IF( .NOT.ln_linssh ) THEN ! variable level thickness : leap-frog on tracer*e3t 147 143 DO jk = 1, jpkm1 148 144 DO jj = 2, jpjm1 … … 154 150 END DO 155 151 END DO 156 ELSE ! fixed level thickness : leap-frog on tracers152 ELSE ! fixed level thickness : leap-frog on tracers 157 153 DO jk = 1, jpkm1 158 154 DO jj = 2, jpjm1 … … 166 162 END DO 167 163 ! 168 CALL wrk_dealloc( jpi, jpj, jpk,zwx, zwy )164 CALL wrk_dealloc( jpi,jpj,jpk, zwx, zwy ) 169 165 ! 170 166 IF( nn_timing == 1 ) CALL timing_stop('tra_zdf_exp')
Note: See TracChangeset
for help on using the changeset viewer.