- Timestamp:
- 2015-12-16T10:25:22+01:00 (9 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
branches/2015/dev_merge_2015/NEMOGCM/NEMO/OPA_SRC/LDF/ldfc1d_c2d.F90
r5836 r6060 28 28 29 29 !! * Substitutions 30 # include "domzgr_substitute.h90"31 30 # include "vectopt_loop_substitute.h90" 32 31 !!---------------------------------------------------------------------- … … 72 71 CASE( 'DYN' ) ! T- and F-points 73 72 DO jk = 1, jpk ! pah1 at T-point 74 pah1(:,:,jk) = pahs1(:,:) * ( prat + zc * ( 1._wp + TANH( - ( fsdept(:,:,jk) - zh ) * zw) ) ) * tmask(:,:,jk)73 pah1(:,:,jk) = pahs1(:,:) * ( prat + zc * ( 1._wp + TANH( - ( gdept_n(:,:,jk) - zh ) * zw) ) ) * tmask(:,:,jk) 75 74 END DO 76 75 DO jk = 1, jpk ! pah2 at F-point (zdep2 is an approximation in zps-coord.) 77 76 DO jj = 1, jpjm1 78 77 DO ji = 1, fs_jpim1 79 zdep2 = ( fsdept(ji,jj+1,jk) + fsdept(ji+1,jj+1,jk) &80 & + fsdept(ji,jj ,jk) + fsdept(ji+1,jj ,jk) ) * 0.25_wp78 zdep2 = ( gdept_n(ji,jj+1,jk) + gdept_n(ji+1,jj+1,jk) & 79 & + gdept_n(ji,jj ,jk) + gdept_n(ji+1,jj ,jk) ) * 0.25_wp 81 80 pah2(ji,jj,jk) = pahs2(ji,jj) * ( prat + zc * ( 1._wp + TANH( - ( zdep2 - zh ) * zw) ) ) * fmask(ji,jj,jk) 82 81 END DO … … 89 88 DO jj = 1, jpjm1 90 89 DO ji = 1, fs_jpim1 91 zdep1 = ( fsdept(ji,jj,jk) + fsdept(ji+1,jj,jk) ) * 0.5_wp92 zdep2 = ( fsdept(ji,jj,jk) + fsdept(ji,jj+1,jk) ) * 0.5_wp90 zdep1 = ( gdept_n(ji,jj,jk) + gdept_n(ji+1,jj,jk) ) * 0.5_wp 91 zdep2 = ( gdept_n(ji,jj,jk) + gdept_n(ji,jj+1,jk) ) * 0.5_wp 93 92 pah1(ji,jj,jk) = pahs1(ji,jj) * ( prat + zc * ( 1._wp + TANH( - ( zdep1 - zh ) * zw) ) ) * umask(ji,jj,jk) 94 93 pah2(ji,jj,jk) = pahs2(ji,jj) * ( prat + zc * ( 1._wp + TANH( - ( zdep2 - zh ) * zw) ) ) * vmask(ji,jj,jk)
Note: See TracChangeset
for help on using the changeset viewer.