- Timestamp:
- 2015-02-09T14:39:07+01:00 (9 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
branches/2015/dev_r5044_CNRS_LIM3CLEAN/NEMOGCM/NEMO/LIM_SRC_3/limhdf.F90
r5047 r5070 72 72 DO jj = 2, jpjm1 73 73 DO ji = fs_2 , fs_jpim1 ! vector opt. 74 efact(ji,jj) = ( e2u(ji,jj) + e2u(ji-1,jj) + e1v(ji,jj) + e1v(ji,jj-1) ) / ( e1t(ji,jj) * e2t(ji,jj))74 efact(ji,jj) = ( e2u(ji,jj) + e2u(ji-1,jj) + e1v(ji,jj) + e1v(ji,jj-1) ) * r1_e12t(ji,jj) 75 75 END DO 76 76 END DO … … 98 98 DO jj = 1, jpjm1 ! diffusive fluxes in U- and V- direction 99 99 DO ji = 1 , fs_jpim1 ! vector opt. 100 zflu(ji,jj) = pahu(ji,jj) * e2u(ji,jj) /e1u(ji,jj) * ( ptab(ji+1,jj) - ptab(ji,jj) )101 zflv(ji,jj) = pahv(ji,jj) * e1v(ji,jj) /e2v(ji,jj) * ( ptab(ji,jj+1) - ptab(ji,jj) )100 zflu(ji,jj) = pahu(ji,jj) * e2u(ji,jj) * r1_e1u(ji,jj) * ( ptab(ji+1,jj) - ptab(ji,jj) ) 101 zflv(ji,jj) = pahv(ji,jj) * e1v(ji,jj) * r1_e2v(ji,jj) * ( ptab(ji,jj+1) - ptab(ji,jj) ) 102 102 END DO 103 103 END DO … … 105 105 DO jj= 2, jpjm1 ! diffusive trend : divergence of the fluxes 106 106 DO ji = fs_2 , fs_jpim1 ! vector opt. 107 zdiv (ji,jj) = ( zflu(ji,jj) - zflu(ji-1,jj ) & 108 & + zflv(ji,jj) - zflv(ji ,jj-1) ) / ( e1t (ji,jj) * e2t (ji,jj) ) 107 zdiv(ji,jj) = ( zflu(ji,jj) - zflu(ji-1,jj) + zflv(ji,jj) - zflv(ji,jj-1) ) * r1_e12t(ji,jj) 109 108 END DO 110 109 END DO … … 116 115 zrlxint = ( ztab0(ji,jj) & 117 116 & + rdt_ice * ( zalfa * ( zdiv(ji,jj) + efact(ji,jj) * ptab(ji,jj) ) & 118 & + ( 1.0 - zalfa ) * zdiv0(ji,jj) ) )&119 & / ( 1.0 + zalfa * rdt_ice * efact(ji,jj) )117 & + ( 1.0 - zalfa ) * zdiv0(ji,jj) ) & 118 & ) / ( 1.0 + zalfa * rdt_ice * efact(ji,jj) ) 120 119 zrlx(ji,jj) = ptab(ji,jj) + zrelax * ( zrlxint - ptab(ji,jj) ) 121 120 END DO … … 139 138 DO jj = 1, jpjm1 ! diffusive fluxes in U- and V- direction 140 139 DO ji = 1 , fs_jpim1 ! vector opt. 141 zflu(ji,jj) = pahu(ji,jj) * e2u(ji,jj) /e1u(ji,jj) * ( ptab(ji+1,jj) - ptab(ji,jj) )142 zflv(ji,jj) = pahv(ji,jj) * e1v(ji,jj) /e2v(ji,jj) * ( ptab(ji,jj+1) - ptab(ji,jj) )140 zflu(ji,jj) = pahu(ji,jj) * e2u(ji,jj) * r1_e1u(ji,jj) * ( ptab(ji+1,jj) - ptab(ji,jj) ) 141 zflv(ji,jj) = pahv(ji,jj) * e1v(ji,jj) * r1_e2v(ji,jj) * ( ptab(ji,jj+1) - ptab(ji,jj) ) 143 142 END DO 144 143 END DO … … 146 145 DO jj= 2, jpjm1 ! diffusive trend : divergence of the fluxes 147 146 DO ji = fs_2 , fs_jpim1 ! vector opt. 148 zdiv (ji,jj) = ( zflu(ji,jj) - zflu(ji-1,jj ) & 149 & + zflv(ji,jj) - zflv(ji ,jj-1) ) / ( e1t (ji,jj) * e2t (ji,jj) ) 147 zdiv(ji,jj) = ( zflu(ji,jj) - zflu(ji-1,jj) + zflv(ji,jj) - zflv(ji,jj-1) ) * r1_e12t(ji,jj) 150 148 ptab(ji,jj) = ztab0(ji,jj) + 0.5 * ( zdiv(ji,jj) + zdiv0(ji,jj) ) 151 149 END DO
Note: See TracChangeset
for help on using the changeset viewer.