Changeset 5964
- Timestamp:
- 2015-12-01T16:09:09+01:00 (8 years ago)
- Location:
- branches/NERC/dev_r5589_is_oce_cpl/NEMOGCM/NEMO/OPA_SRC
- Files:
-
- 1 deleted
- 2 edited
Legend:
- Unmodified
- Added
- Removed
-
branches/NERC/dev_r5589_is_oce_cpl/NEMOGCM/NEMO/OPA_SRC/DIA/diahsb.F90
r5955 r5964 276 276 IF(lwp) WRITE(numout,*) ' dia_hsb at initial state ' 277 277 IF(lwp) WRITE(numout,*) '~~~~~~~' 278 surf_ini(:,:) = e1 2t(:,:) * tmask_i(:,:) ! initial ocean surface278 surf_ini(:,:) = e1e2t(:,:) * tmask_i(:,:) ! initial ocean surface 279 279 ssh_ini(:,:) = sshn(:,:) ! initial ssh 280 280 DO jk = 1, jpk -
branches/NERC/dev_r5589_is_oce_cpl/NEMOGCM/NEMO/OPA_SRC/DOM/iscplhsb.F90
r5945 r5964 128 128 jip1=ji+1 ; jim1=ji-1 ; jjp1=jj+1 ; jjm1=jj-1 ; 129 129 130 zsum = e1 2t(ji ,jjp1) * tmask(ji ,jjp1,jk) + e12t(ji ,jjm1) * tmask(ji ,jjm1,jk) &131 & + e1 2t(jim1,jj ) * tmask(jim1,jj ,jk) + e12t(jip1,jj ) * tmask(jip1,jj ,jk)130 zsum = e1e2t(ji ,jjp1) * tmask(ji ,jjp1,jk) + e1e2t(ji ,jjm1) * tmask(ji ,jjm1,jk) & 131 & + e1e2t(jim1,jj ) * tmask(jim1,jj ,jk) + e1e2t(jip1,jj ) * tmask(jip1,jj ,jk) 132 132 133 133 IF ( zsum /= 0._wp ) THEN 134 zjip1_ratio = e1 2t(jip1,jj ) * tmask(jip1,jj ,jk) / zsum135 zjim1_ratio = e1 2t(jim1,jj ) * tmask(jim1,jj ,jk) / zsum136 zjjp1_ratio = e1 2t(ji ,jjp1) * tmask(ji ,jjp1,jk) / zsum137 zjjm1_ratio = e1 2t(ji ,jjm1) * tmask(ji ,jjm1,jk) / zsum134 zjip1_ratio = e1e2t(jip1,jj ) * tmask(jip1,jj ,jk) / zsum 135 zjim1_ratio = e1e2t(jim1,jj ) * tmask(jim1,jj ,jk) / zsum 136 zjjp1_ratio = e1e2t(ji ,jjp1) * tmask(ji ,jjp1,jk) / zsum 137 zjjm1_ratio = e1e2t(ji ,jjm1) * tmask(ji ,jjm1,jk) / zsum 138 138 139 139 pvol_flx(ji ,jjp1,jk ) = pvol_flx(ji ,jjp1,jk ) + pvol_flx(ji,jj,jk ) * zjjp1_ratio
Note: See TracChangeset
for help on using the changeset viewer.