- Timestamp:
- 2013-01-23T15:33:04+01:00 (11 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
branches/2012/dev_MERGE_2012/NEMOGCM/NEMO/LIM_SRC_2/limwri_2.F90
r3625 r3764 105 105 !! 106 106 INTEGER :: ji, jj, jf ! dummy loop indices 107 CHARACTER(len = 40) :: clhstnam, clop107 CHARACTER(len = 80) :: clhstnam, clop 108 108 REAL(wp) :: zsto, zjulian, zout, & ! temporary scalars 109 109 & zindh, zinda, zindb, ztmu … … 161 161 zcmo(ji,jj,5) = sist (ji,jj) 162 162 zcmo(ji,jj,6) = fbif (ji,jj) 163 IF (lk_lim2_vp) THEN 163 164 zcmo(ji,jj,7) = zindb * ( u_ice(ji,jj ) * tmu(ji,jj ) + u_ice(ji+1,jj ) * tmu(ji+1,jj ) & 164 165 + u_ice(ji,jj+1) * tmu(ji,jj+1) + u_ice(ji+1,jj+1) * tmu(ji+1,jj+1) ) & … … 168 169 + v_ice(ji,jj+1) * tmu(ji,jj+1) + v_ice(ji+1,jj+1) * tmu(ji+1,jj+1) ) & 169 170 / ztmu 171 ELSE 172 173 zcmo(ji,jj,7) = zindb * ( u_ice(ji,jj ) * tmu(ji,jj) & 174 & + u_ice(ji-1,jj) * tmu(ji-1,jj) ) & 175 & / 2.0 176 zcmo(ji,jj,8) = zindb * ( v_ice(ji,jj ) * tmv(ji,jj) & 177 & + v_ice(ji,jj-1) * tmv(ji,jj-1) ) & 178 & / 2.0 179 180 ENDIF 170 181 zcmo(ji,jj,9) = sst_m(ji,jj) 171 182 zcmo(ji,jj,10) = sss_m(ji,jj) … … 187 198 niter = niter + 1 188 199 DO jf = 1 , noumef 189 DO jj = 1 , jpj 190 DO ji = 1 , jpi 191 zfield(ji,jj) = zcmo(ji,jj,jf) * cmulti(jf) + cadd(jf) 192 END DO 193 END DO 194 195 IF( jf == 7 .OR. jf == 8 .OR. jf == 15 .OR. jf == 16 ) THEN 200 zfield(:,:) = zcmo(:,:,jf) * cmulti(jf) + cadd(jf) * tmask(:,:,1) 201 SELECT CASE ( jf ) 202 CASE ( 7, 8, 15, 16, 20, 21 ) ! velocity or stress fields (vectors) 196 203 CALL lbc_lnk( zfield, 'T', -1. ) 197 ELSE204 CASE DEFAULT ! scalar fields 198 205 CALL lbc_lnk( zfield, 'T', 1. ) 199 END IF200 206 END SELECT 207 201 208 IF( nc(jf) == 1 ) CALL histwrite( nice, nam(jf), niter, zfield, ndim, ndex51 ) 202 209 203 210 END DO 204 211 205 212 IF( ( nn_fsbc * niter ) >= nitend ) CALL histclo( nice ) 206 213 … … 209 216 END SUBROUTINE lim_wri_2 210 217 211 # endif218 #endif 212 219 213 220 SUBROUTINE lim_wri_init_2
Note: See TracChangeset
for help on using the changeset viewer.