Changeset 5951 for branches/2014/dev_r4650_UKMO14.4_OBS_GENERAL_VINTERP/NEMOGCM/NEMO/LIM_SRC_3/limtrp.F90
- Timestamp:
- 2015-11-30T12:48:01+01:00 (8 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
branches/2014/dev_r4650_UKMO14.4_OBS_GENERAL_VINTERP/NEMOGCM/NEMO/LIM_SRC_3/limtrp.F90
r5950 r5951 95 95 ENDIF 96 96 97 zsm(:,:) = e1 2t(:,:)97 zsm(:,:) = e1e2t(:,:) 98 98 99 99 ! !-------------------------------------! … … 162 162 ! transported fields 163 163 !------------------------- 164 z0opw(:,:,1) = ato_i(:,:) * e1 2t(:,:) ! Open water area165 DO jl = 1, jpl 166 z0snw (:,:,jl) = v_s (:,:, jl) * e12t(:,:)! Snow volume167 z0ice(:,:,jl) = v_i (:,:, jl) * e12t(:,:)! Ice volume168 z0ai (:,:,jl) = a_i (:,:, jl) * e12t(:,:)! Ice area169 z0smi (:,:,jl) = smv_i(:,:, jl) * e12t(:,:)! Salt content170 z0oi (:,:,jl) = oa_i (:,:, jl) * e12t(:,:)! Age content171 z0es (:,:,jl) = e_s (:,:,1,jl) * e1 2t(:,:) ! Snow heat content164 z0opw(:,:,1) = ato_i(:,:) * e1e2t(:,:) ! Open water area 165 DO jl = 1, jpl 166 z0snw (:,:,jl) = v_s (:,:, jl) * e1e2t(:,:) ! Snow volume 167 z0ice(:,:,jl) = v_i (:,:, jl) * e1e2t(:,:) ! Ice volume 168 z0ai (:,:,jl) = a_i (:,:, jl) * e1e2t(:,:) ! Ice area 169 z0smi (:,:,jl) = smv_i(:,:, jl) * e1e2t(:,:) ! Salt content 170 z0oi (:,:,jl) = oa_i (:,:, jl) * e1e2t(:,:) ! Age content 171 z0es (:,:,jl) = e_s (:,:,1,jl) * e1e2t(:,:) ! Snow heat content 172 172 DO jk = 1, nlay_i 173 z0ei (:,:,jk,jl) = e_i (:,:,jk,jl) * e1 2t(:,:) ! Ice heat content173 z0ei (:,:,jk,jl) = e_i (:,:,jk,jl) * e1e2t(:,:) ! Ice heat content 174 174 END DO 175 175 END DO … … 263 263 ! Recover the properties from their contents 264 264 !------------------------------------------- 265 ato_i(:,:) = z0opw(:,:,1) * r1_e1 2t(:,:)266 DO jl = 1, jpl 267 v_i (:,:, jl) = z0ice(:,:,jl) * r1_e12t(:,:)268 v_s (:,:, jl) = z0snw(:,:,jl) * r1_e12t(:,:)269 smv_i(:,:, jl) = z0smi(:,:,jl) * r1_e12t(:,:)270 oa_i (:,:, jl) = z0oi (:,:,jl) * r1_e12t(:,:)271 a_i (:,:, jl) = z0ai (:,:,jl) * r1_e12t(:,:)272 e_s (:,:,1,jl) = z0es (:,:,jl) * r1_e1 2t(:,:)265 ato_i(:,:) = z0opw(:,:,1) * r1_e1e2t(:,:) 266 DO jl = 1, jpl 267 v_i (:,:, jl) = z0ice(:,:,jl) * r1_e1e2t(:,:) 268 v_s (:,:, jl) = z0snw(:,:,jl) * r1_e1e2t(:,:) 269 smv_i(:,:, jl) = z0smi(:,:,jl) * r1_e1e2t(:,:) 270 oa_i (:,:, jl) = z0oi (:,:,jl) * r1_e1e2t(:,:) 271 a_i (:,:, jl) = z0ai (:,:,jl) * r1_e1e2t(:,:) 272 e_s (:,:,1,jl) = z0es (:,:,jl) * r1_e1e2t(:,:) 273 273 DO jk = 1, nlay_i 274 e_i(:,:,jk,jl) = z0ei(:,:,jk,jl) * r1_e1 2t(:,:)274 e_i(:,:,jk,jl) = z0ei(:,:,jk,jl) * r1_e1e2t(:,:) 275 275 END DO 276 276 END DO
Note: See TracChangeset
for help on using the changeset viewer.