- Timestamp:
- 2013-11-20T11:17:17+01:00 (10 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
branches/2013/dev_r3858_NOC_ZTC/NEMOGCM/NEMO/TOP_SRC/trcsub.F90
r4258 r4280 79 79 ! 80 80 sshb_hold (:,:) = sshn (:,:) 81 82 81 !!Z~ sshu_b_hold(:,:) = sshu_n(:,:) 82 !!Z~ sshv_b_hold(:,:) = sshv_n(:,:) 83 83 emp_b_hold (:,:) = emp_b (:,:) 84 84 ! … … 117 117 ! 118 118 sshn_tm (:,:) = sshn_tm (:,:) + sshn (:,:) 119 120 119 !!Z~ sshu_n_tm(:,:) = sshu_n_tm(:,:) + sshu_n(:,:) 120 !!Z~ sshv_n_tm(:,:) = sshv_n_tm(:,:) + sshv_n(:,:) 121 121 rnf_tm (:,:) = rnf_tm (:,:) + rnf (:,:) 122 122 h_rnf_tm (:,:) = h_rnf_tm (:,:) + h_rnf (:,:) … … 197 197 # endif 198 198 sshn_temp (:,:) = sshn (:,:) 199 200 201 199 !!Z~ sshu_n_temp(:,:) = sshu_n(:,:) 200 !!Z~ sshv_n_temp(:,:) = sshv_n(:,:) 201 !!Z~ sshf_n_temp(:,:) = sshf_n(:,:) 202 202 sshb_temp (:,:) = sshb (:,:) 203 204 203 !!Z~ sshu_b_temp(:,:) = sshu_b(:,:) 204 !!Z~ sshv_b_temp(:,:) = sshv_b(:,:) 205 205 ssha_temp (:,:) = ssha (:,:) 206 207 206 !!Z~ sshu_a_temp(:,:) = sshu_a(:,:) 207 !!Z~ sshv_a_temp(:,:) = sshv_a(:,:) 208 208 rnf_temp (:,:) = rnf (:,:) 209 209 h_rnf_temp (:,:) = h_rnf (:,:) … … 309 309 # endif 310 310 sshn_tm (:,:) = sshn_tm (:,:) + sshn (:,:) 311 312 311 !!Z~ sshu_n_tm(:,:) = sshu_n_tm (:,:) + sshu_n(:,:) 312 !!Z~ sshv_n_tm(:,:) = sshv_n_tm (:,:) + sshv_n(:,:) 313 313 rnf_tm (:,:) = rnf_tm (:,:) + rnf (:,:) 314 314 h_rnf_tm (:,:) = h_rnf_tm (:,:) + h_rnf (:,:) … … 321 321 ! 322 322 sshn (:,:) = sshn_tm (:,:) * r1_ndttrcp1 323 324 323 !!Z~ sshu_n (:,:) = sshu_n_tm (:,:) * r1_ndttrcp1 324 !!Z~ sshv_n (:,:) = sshv_n_tm (:,:) * r1_ndttrcp1 325 325 sshb (:,:) = sshb_hold (:,:) 326 327 326 !!Z~ sshu_b (:,:) = sshu_b_hold(:,:) 327 !!Z~ sshv_b (:,:) = sshv_b_hold(:,:) 328 328 rnf (:,:) = rnf_tm (:,:) * r1_ndttrcp1 329 329 h_rnf (:,:) = h_rnf_tm (:,:) * r1_ndttrcp1 … … 486 486 #endif 487 487 CALL lbc_lnk( sshn (:,:) , 'T', 1. ) 488 489 490 488 !!Z~ CALL lbc_lnk( sshu_n(:,:) , 'U', 1. ) 489 !!Z~ CALL lbc_lnk( sshv_n(:,:) , 'V', 1. ) 490 !!Z~ CALL lbc_lnk( sshf_n(:,:) , 'F', 1. ) 491 491 CALL lbc_lnk( sshb (:,:) , 'T', 1. ) 492 493 492 !!Z~ CALL lbc_lnk( sshu_b(:,:) , 'U', 1. ) 493 !!Z~ CALL lbc_lnk( sshv_b(:,:) , 'V', 1. ) 494 494 CALL lbc_lnk( ssha (:,:) , 'T', 1. ) 495 496 495 !!Z~ CALL lbc_lnk( sshu_a(:,:) , 'U', 1. ) 496 !!Z~ CALL lbc_lnk( sshv_a(:,:) , 'V', 1. ) 497 497 CALL lbc_lnk( rnf (:,:) , 'T', 1. ) 498 498 CALL lbc_lnk( h_rnf (:,:) , 'T', 1. ) … … 592 592 #endif 593 593 sshn_tm (:,:) = sshn (:,:) 594 595 594 !!Z~ sshu_n_tm(:,:) = sshu_n(:,:) 595 !!Z~ sshv_n_tm(:,:) = sshv_n(:,:) 596 596 rnf_tm (:,:) = rnf (:,:) 597 597 h_rnf_tm (:,:) = h_rnf (:,:) … … 695 695 sshb (:,:) = sshb_temp (:,:) 696 696 ssha (:,:) = ssha_temp (:,:) 697 698 699 700 701 702 703 697 !!Z~ sshu_n(:,:) = sshu_n_temp(:,:) 698 !!Z~ sshu_b(:,:) = sshu_b_temp(:,:) 699 !!Z~ sshu_a(:,:) = sshu_a_temp(:,:) 700 !!Z~ sshv_n(:,:) = sshv_n_temp(:,:) 701 !!Z~ sshv_b(:,:) = sshv_b_temp(:,:) 702 !!Z~ sshv_a(:,:) = sshv_a_temp(:,:) 703 !!Z~ sshf_n(:,:) = sshf_n_temp(:,:) 704 704 rnf (:,:) = rnf_temp (:,:) 705 705 h_rnf (:,:) = h_rnf_temp (:,:) … … 816 816 ! 817 817 sshb_hold (:,:) = sshn (:,:) 818 819 818 !!Z~ sshu_b_hold(:,:) = sshu_n(:,:) 819 !!Z~ sshv_b_hold(:,:) = sshv_n(:,:) 820 820 emp_b_hold (:,:) = emp (:,:) 821 821 sshn_tm (:,:) = sshn (:,:) 822 823 822 !!Z~ sshu_n_tm (:,:) = sshu_n(:,:) 823 !!Z~ sshv_n_tm (:,:) = sshv_n(:,:) 824 824 rnf_tm (:,:) = rnf (:,:) 825 825 h_rnf_tm (:,:) = h_rnf (:,:) … … 951 951 END DO 952 952 ! 953 954 953 !!Z~ hu(:,:) = hu_0(:,:) + sshu_n(:,:) ! now ocean depth (at u- and v-points) 954 !!Z~ hv(:,:) = hv_0(:,:) + sshv_n(:,:) 955 955 ! ! now masked inverse of the ocean depth (at u- and v-points) 956 956 hur(:,:) = umask(:,:,1) / ( hu(:,:) + 1._wp - umask(:,:,1) ) … … 992 992 993 993 ! ! Sea Surface Height at u-,v- and f-points (vvl case only) 994 995 996 997 998 999 1000 1001 1002 1003 1004 1005 1006 994 !!Z~ IF( lk_vvl ) THEN ! (required only in key_vvl case) 995 !!Z~ DO jj = 1, jpjm1 996 !!Z~ DO ji = 1, jpim1 ! NO Vector Opt. 997 !!Z~ sshu_a(ji,jj) = 0.5 * umask(ji,jj,1) / ( e1u(ji ,jj) * e2u(ji ,jj) ) & 998 !!Z~ & * ( e1t(ji ,jj) * e2t(ji ,jj) * ssha(ji ,jj) & 999 !!Z~ & + e1t(ji+1,jj) * e2t(ji+1,jj) * ssha(ji+1,jj) ) 1000 !!Z~ sshv_a(ji,jj) = 0.5 * vmask(ji,jj,1) / ( e1v(ji,jj ) * e2v(ji,jj ) ) & 1001 !!Z~ & * ( e1t(ji,jj ) * e2t(ji,jj ) * ssha(ji,jj ) & 1002 !!Z~ & + e1t(ji,jj+1) * e2t(ji,jj+1) * ssha(ji,jj+1) ) 1003 !!Z~ END DO 1004 !!Z~ END DO 1005 !!Z~ CALL lbc_lnk( sshu_a, 'U', 1. ) ; CALL lbc_lnk( sshv_a, 'V', 1. ) ! Boundaries conditions 1006 !!Z~ ENDIF 1007 1007 1008 1008
Note: See TracChangeset
for help on using the changeset viewer.