- Timestamp:
- 2020-09-15T12:56:56+02:00 (4 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
NEMO/branches/2020/temporary_r4_trunk/src/ICE/icedyn_adv_umx.F90
r13469 r13470 114 114 END WHERE 115 115 DO jl = 1, jpl 116 DO_2D _00_00116 DO_2D( 0, 0, 0, 0 ) 117 117 zhip_max(ji,jj,jl) = MAX( epsi20, ph_ip(ji,jj,jl), ph_ip(ji+1,jj ,jl), ph_ip(ji ,jj+1,jl), & 118 118 & ph_ip(ji-1,jj ,jl), ph_ip(ji ,jj-1,jl), & … … 147 147 END DO 148 148 DO jl = 1, jpl 149 DO_3D _00_00(1, nlay_i )149 DO_3D( 0, 0, 0, 0, 1, nlay_i ) 150 150 zei_max(ji,jj,jk,jl) = MAX( epsi20, ze_i(ji,jj,jk,jl), ze_i(ji+1,jj ,jk,jl), ze_i(ji ,jj+1,jk,jl), & 151 151 & ze_i(ji-1,jj ,jk,jl), ze_i(ji ,jj-1,jk,jl), & … … 155 155 END DO 156 156 DO jl = 1, jpl 157 DO_3D _00_00(1, nlay_s )157 DO_3D( 0, 0, 0, 0, 1, nlay_s ) 158 158 zes_max(ji,jj,jk,jl) = MAX( epsi20, ze_s(ji,jj,jk,jl), ze_s(ji+1,jj ,jk,jl), ze_s(ji ,jj+1,jk,jl), & 159 159 & ze_s(ji-1,jj ,jk,jl), ze_s(ji ,jj-1,jk,jl), & … … 191 191 ! 192 192 ! --- define velocity for advection: u*grad(H) --- ! 193 DO_2D _00_00193 DO_2D( 0, 0, 0, 0 ) 194 194 IF ( pu_ice(ji,jj) * pu_ice(ji-1,jj) <= 0._wp ) THEN ; zcu_box(ji,jj) = 0._wp 195 195 ELSEIF( pu_ice(ji,jj) > 0._wp ) THEN ; zcu_box(ji,jj) = pu_ice(ji-1,jj) … … 224 224 IF( .NOT. ALLOCATED(jmsk_small) ) ALLOCATE( jmsk_small(jpi,jpj,jpl) ) 225 225 DO jl = 1, jpl 226 DO_2D _10_10226 DO_2D( 1, 0, 1, 0 ) 227 227 zvi_cen = 0.5_wp * ( pv_i(ji+1,jj,jl) + pv_i(ji,jj,jl) ) 228 228 IF( zvi_cen < epsi06) THEN ; imsk_small(ji,jj,jl) = 0 … … 380 380 !== Open water area ==! 381 381 zati2(:,:) = SUM( pa_i(:,:,:), dim=3 ) 382 DO_2D _00_00382 DO_2D( 0, 0, 0, 0 ) 383 383 pato_i(ji,jj) = pato_i(ji,jj) - ( zati2(ji,jj) - zati1(ji,jj) ) & 384 384 & - ( zudy(ji,jj) - zudy(ji-1,jj) + zvdx(ji,jj) - zvdx(ji,jj-1) ) * r1_e1e2t(ji,jj) * zdt … … 490 490 IF( pamsk == 0._wp ) THEN 491 491 DO jl = 1, jpl 492 DO_2D _10_10492 DO_2D( 1, 0, 1, 0 ) 493 493 IF( ABS( pu(ji,jj) ) > epsi10 ) THEN 494 494 zfu_ho (ji,jj,jl) = zfu_ho (ji,jj,jl) * puc (ji,jj,jl) / pu(ji,jj) … … 512 512 ! thus we calculate the upstream solution and apply a limiter again 513 513 DO jl = 1, jpl 514 DO_2D _00_00514 DO_2D( 0, 0, 0, 0 ) 515 515 ztra = - ( zfu_ups(ji,jj,jl) - zfu_ups(ji-1,jj,jl) + zfv_ups(ji,jj,jl) - zfv_ups(ji,jj-1,jl) ) 516 516 ! … … 533 533 IF( PRESENT( pua_ho ) ) THEN 534 534 DO jl = 1, jpl 535 DO_2D _10_10535 DO_2D( 1, 0, 1, 0 ) 536 536 pua_ho (ji,jj,jl) = zfu_ho (ji,jj,jl) ; pva_ho (ji,jj,jl) = zfv_ho (ji,jj,jl) 537 537 pua_ups(ji,jj,jl) = zfu_ups(ji,jj,jl) ; pva_ups(ji,jj,jl) = zfv_ups(ji,jj,jl) … … 543 543 ! --------------------------------- 544 544 DO jl = 1, jpl 545 DO_2D _00_00545 DO_2D( 0, 0, 0, 0 ) 546 546 ztra = - ( zfu_ho(ji,jj,jl) - zfu_ho(ji-1,jj,jl) + zfv_ho(ji,jj,jl) - zfv_ho(ji,jj-1,jl) ) 547 547 ! … … 577 577 ! 578 578 DO jl = 1, jpl 579 DO_2D _10_10579 DO_2D( 1, 0, 1, 0 ) 580 580 pfu_ups(ji,jj,jl) = MAX( pu(ji,jj), 0._wp ) * pt(ji,jj,jl) + MIN( pu(ji,jj), 0._wp ) * pt(ji+1,jj,jl) 581 581 pfv_ups(ji,jj,jl) = MAX( pv(ji,jj), 0._wp ) * pt(ji,jj,jl) + MIN( pv(ji,jj), 0._wp ) * pt(ji,jj+1,jl) … … 588 588 ! 589 589 DO jl = 1, jpl !-- flux in x-direction 590 DO_2D _10_10590 DO_2D( 1, 0, 1, 0 ) 591 591 pfu_ups(ji,jj,jl) = MAX( pu(ji,jj), 0._wp ) * pt(ji,jj,jl) + MIN( pu(ji,jj), 0._wp ) * pt(ji+1,jj,jl) 592 592 END_2D … … 594 594 ! 595 595 DO jl = 1, jpl !-- first guess of tracer from u-flux 596 DO_2D _00_00596 DO_2D( 0, 0, 0, 0 ) 597 597 ztra = - ( pfu_ups(ji,jj,jl) - pfu_ups(ji-1,jj,jl) ) & 598 598 & + ( pu (ji,jj ) - pu (ji-1,jj ) ) * pt(ji,jj,jl) * (1.-pamsk) … … 604 604 ! 605 605 DO jl = 1, jpl !-- flux in y-direction 606 DO_2D _10_10606 DO_2D( 1, 0, 1, 0 ) 607 607 pfv_ups(ji,jj,jl) = MAX( pv(ji,jj), 0._wp ) * zpt(ji,jj,jl) + MIN( pv(ji,jj), 0._wp ) * zpt(ji,jj+1,jl) 608 608 END_2D … … 612 612 ! 613 613 DO jl = 1, jpl !-- flux in y-direction 614 DO_2D _10_10614 DO_2D( 1, 0, 1, 0 ) 615 615 pfv_ups(ji,jj,jl) = MAX( pv(ji,jj), 0._wp ) * pt(ji,jj,jl) + MIN( pv(ji,jj), 0._wp ) * pt(ji,jj+1,jl) 616 616 END_2D … … 618 618 ! 619 619 DO jl = 1, jpl !-- first guess of tracer from v-flux 620 DO_2D _00_00620 DO_2D( 0, 0, 0, 0 ) 621 621 ztra = - ( pfv_ups(ji,jj,jl) - pfv_ups(ji,jj-1,jl) ) & 622 622 & + ( pv (ji,jj ) - pv (ji,jj-1 ) ) * pt(ji,jj,jl) * (1.-pamsk) … … 628 628 ! 629 629 DO jl = 1, jpl !-- flux in x-direction 630 DO_2D _10_10630 DO_2D( 1, 0, 1, 0 ) 631 631 pfu_ups(ji,jj,jl) = MAX( pu(ji,jj), 0._wp ) * zpt(ji,jj,jl) + MIN( pu(ji,jj), 0._wp ) * zpt(ji+1,jj,jl) 632 632 END_2D … … 638 638 ! 639 639 DO jl = 1, jpl !-- after tracer with upstream scheme 640 DO_2D _00_00640 DO_2D( 0, 0, 0, 0 ) 641 641 ztra = - ( pfu_ups(ji,jj,jl) - pfu_ups(ji-1,jj ,jl) & 642 642 & + pfv_ups(ji,jj,jl) - pfv_ups(ji ,jj-1,jl) ) & … … 677 677 ! 678 678 DO jl = 1, jpl 679 DO_2D _10_10679 DO_2D( 1, 0, 1, 0 ) 680 680 pfu_ho(ji,jj,jl) = 0.5_wp * pu(ji,jj) * ( pt(ji,jj,jl) + pt(ji+1,jj ,jl) ) 681 681 pfv_ho(ji,jj,jl) = 0.5_wp * pv(ji,jj) * ( pt(ji,jj,jl) + pt(ji ,jj+1,jl) ) … … 695 695 ! 696 696 DO jl = 1, jpl !-- flux in x-direction 697 DO_2D _10_10697 DO_2D( 1, 0, 1, 0 ) 698 698 pfu_ho(ji,jj,jl) = 0.5_wp * pu(ji,jj) * ( pt(ji,jj,jl) + pt(ji+1,jj,jl) ) 699 699 END_2D … … 702 702 703 703 DO jl = 1, jpl !-- first guess of tracer from u-flux 704 DO_2D _00_00704 DO_2D( 0, 0, 0, 0 ) 705 705 ztra = - ( pfu_ho(ji,jj,jl) - pfu_ho(ji-1,jj,jl) ) & 706 706 & + ( pu (ji,jj ) - pu (ji-1,jj ) ) * pt(ji,jj,jl) * (1.-pamsk) … … 712 712 713 713 DO jl = 1, jpl !-- flux in y-direction 714 DO_2D _10_10714 DO_2D( 1, 0, 1, 0 ) 715 715 pfv_ho(ji,jj,jl) = 0.5_wp * pv(ji,jj) * ( zpt(ji,jj,jl) + zpt(ji,jj+1,jl) ) 716 716 END_2D … … 721 721 ! 722 722 DO jl = 1, jpl !-- flux in y-direction 723 DO_2D _10_10723 DO_2D( 1, 0, 1, 0 ) 724 724 pfv_ho(ji,jj,jl) = 0.5_wp * pv(ji,jj) * ( pt(ji,jj,jl) + pt(ji,jj+1,jl) ) 725 725 END_2D … … 728 728 ! 729 729 DO jl = 1, jpl !-- first guess of tracer from v-flux 730 DO_2D _00_00730 DO_2D( 0, 0, 0, 0 ) 731 731 ztra = - ( pfv_ho(ji,jj,jl) - pfv_ho(ji,jj-1,jl) ) & 732 732 & + ( pv (ji,jj ) - pv (ji,jj-1 ) ) * pt(ji,jj,jl) * (1.-pamsk) … … 738 738 ! 739 739 DO jl = 1, jpl !-- flux in x-direction 740 DO_2D _10_10740 DO_2D( 1, 0, 1, 0 ) 741 741 pfu_ho(ji,jj,jl) = 0.5_wp * pu(ji,jj) * ( zpt(ji,jj,jl) + zpt(ji+1,jj,jl) ) 742 742 END_2D … … 786 786 ! !-- advective form update in zpt --! 787 787 DO jl = 1, jpl 788 DO_2D _00_00788 DO_2D( 0, 0, 0, 0 ) 789 789 zpt(ji,jj,jl) = ( pt(ji,jj,jl) - ( pubox(ji,jj ) * ( zt_u(ji,jj,jl) - zt_u(ji-1,jj,jl) ) * r1_e1t (ji,jj) & 790 790 & + pt (ji,jj,jl) * ( pu (ji,jj ) - pu (ji-1,jj ) ) * r1_e1e2t(ji,jj) & … … 813 813 ! !-- advective form update in zpt --! 814 814 DO jl = 1, jpl 815 DO_2D _00_00815 DO_2D( 0, 0, 0, 0 ) 816 816 zpt(ji,jj,jl) = ( pt(ji,jj,jl) - ( pvbox(ji,jj ) * ( zt_v(ji,jj,jl) - zt_v(ji,jj-1,jl) ) * r1_e2t (ji,jj) & 817 817 & + pt (ji,jj,jl) * ( pv (ji,jj ) - pv (ji,jj-1 ) ) * r1_e1e2t(ji,jj) & … … 895 895 ! 896 896 DO jl = 1, jpl 897 DO_2D _10_10897 DO_2D( 1, 0, 1, 0 ) 898 898 pt_u(ji,jj,jl) = 0.5_wp * umask(ji,jj,1) * ( pt(ji+1,jj,jl) + pt(ji,jj,jl) & 899 899 & - SIGN( 1._wp, pu(ji,jj) ) * ( pt(ji+1,jj,jl) - pt(ji,jj,jl) ) ) … … 904 904 ! 905 905 DO jl = 1, jpl 906 DO_2D _10_10906 DO_2D( 1, 0, 1, 0 ) 907 907 zcu = pu(ji,jj) * r1_e2u(ji,jj) * pdt * r1_e1u(ji,jj) 908 908 pt_u(ji,jj,jl) = 0.5_wp * umask(ji,jj,1) * ( pt(ji+1,jj,jl) + pt(ji,jj,jl) & … … 914 914 ! 915 915 DO jl = 1, jpl 916 DO_2D _10_10916 DO_2D( 1, 0, 1, 0 ) 917 917 zcu = pu(ji,jj) * r1_e2u(ji,jj) * pdt * r1_e1u(ji,jj) 918 918 zdx2 = e1u(ji,jj) * e1u(ji,jj) … … 928 928 ! 929 929 DO jl = 1, jpl 930 DO_2D _10_10930 DO_2D( 1, 0, 1, 0 ) 931 931 zcu = pu(ji,jj) * r1_e2u(ji,jj) * pdt * r1_e1u(ji,jj) 932 932 zdx2 = e1u(ji,jj) * e1u(ji,jj) … … 942 942 ! 943 943 DO jl = 1, jpl 944 DO_2D _10_10944 DO_2D( 1, 0, 1, 0 ) 945 945 zcu = pu(ji,jj) * r1_e2u(ji,jj) * pdt * r1_e1u(ji,jj) 946 946 zdx2 = e1u(ji,jj) * e1u(ji,jj) … … 963 963 IF( ll_neg ) THEN 964 964 DO jl = 1, jpl 965 DO_2D _10_10965 DO_2D( 1, 0, 1, 0 ) 966 966 IF( pt_u(ji,jj,jl) < 0._wp .OR. ( imsk_small(ji,jj,jl) == 0 .AND. pamsk == 0. ) ) THEN 967 967 pt_u(ji,jj,jl) = 0.5_wp * umask(ji,jj,1) * ( pt(ji+1,jj,jl) + pt(ji,jj,jl) & … … 973 973 ! !-- High order flux in i-direction --! 974 974 DO jl = 1, jpl 975 DO_2D _10_10975 DO_2D( 1, 0, 1, 0 ) 976 976 pfu_ho(ji,jj,jl) = pu(ji,jj) * pt_u(ji,jj,jl) 977 977 END_2D … … 1006 1006 ! !-- Laplacian in j-direction --! 1007 1007 DO jl = 1, jpl 1008 DO_2D _10_001008 DO_2D( 1, 0, 0, 0 ) 1009 1009 ztv1(ji,jj,jl) = ( pt(ji,jj+1,jl) - pt(ji,jj,jl) ) * r1_e2v(ji,jj) * vmask(ji,jj,1) 1010 1010 END_2D 1011 DO_2D _00_001011 DO_2D( 0, 0, 0, 0 ) 1012 1012 ztv2(ji,jj,jl) = ( ztv1(ji,jj,jl) - ztv1(ji,jj-1,jl) ) * r1_e2t(ji,jj) 1013 1013 END_2D … … 1017 1017 ! !-- BiLaplacian in j-direction --! 1018 1018 DO jl = 1, jpl 1019 DO_2D _10_001019 DO_2D( 1, 0, 0, 0 ) 1020 1020 ztv3(ji,jj,jl) = ( ztv2(ji,jj+1,jl) - ztv2(ji,jj,jl) ) * r1_e2v(ji,jj) * vmask(ji,jj,1) 1021 1021 END_2D 1022 DO_2D _00_001022 DO_2D( 0, 0, 0, 0 ) 1023 1023 ztv4(ji,jj,jl) = ( ztv3(ji,jj,jl) - ztv3(ji,jj-1,jl) ) * r1_e2t(ji,jj) 1024 1024 END_2D … … 1031 1031 CASE( 1 ) !== 1st order central TIM ==! (Eq. 21) 1032 1032 DO jl = 1, jpl 1033 DO_2D _10_101033 DO_2D( 1, 0, 1, 0 ) 1034 1034 pt_v(ji,jj,jl) = 0.5_wp * vmask(ji,jj,1) * ( pt(ji,jj+1,jl) + pt(ji,jj,jl) & 1035 1035 & - SIGN( 1._wp, pv(ji,jj) ) * ( pt(ji,jj+1,jl) - pt(ji,jj,jl) ) ) … … 1039 1039 CASE( 2 ) !== 2nd order central TIM ==! (Eq. 23) 1040 1040 DO jl = 1, jpl 1041 DO_2D _10_101041 DO_2D( 1, 0, 1, 0 ) 1042 1042 zcv = pv(ji,jj) * r1_e1v(ji,jj) * pdt * r1_e2v(ji,jj) 1043 1043 pt_v(ji,jj,jl) = 0.5_wp * vmask(ji,jj,1) * ( pt(ji,jj+1,jl) + pt(ji,jj,jl) & … … 1048 1048 CASE( 3 ) !== 3rd order central TIM ==! (Eq. 24) 1049 1049 DO jl = 1, jpl 1050 DO_2D _10_101050 DO_2D( 1, 0, 1, 0 ) 1051 1051 zcv = pv(ji,jj) * r1_e1v(ji,jj) * pdt * r1_e2v(ji,jj) 1052 1052 zdy2 = e2v(ji,jj) * e2v(ji,jj) … … 1061 1061 CASE( 4 ) !== 4th order central TIM ==! (Eq. 27) 1062 1062 DO jl = 1, jpl 1063 DO_2D _10_101063 DO_2D( 1, 0, 1, 0 ) 1064 1064 zcv = pv(ji,jj) * r1_e1v(ji,jj) * pdt * r1_e2v(ji,jj) 1065 1065 zdy2 = e2v(ji,jj) * e2v(ji,jj) … … 1074 1074 CASE( 5 ) !== 5th order central TIM ==! (Eq. 29) 1075 1075 DO jl = 1, jpl 1076 DO_2D _10_101076 DO_2D( 1, 0, 1, 0 ) 1077 1077 zcv = pv(ji,jj) * r1_e1v(ji,jj) * pdt * r1_e2v(ji,jj) 1078 1078 zdy2 = e2v(ji,jj) * e2v(ji,jj) … … 1095 1095 IF( ll_neg ) THEN 1096 1096 DO jl = 1, jpl 1097 DO_2D _10_101097 DO_2D( 1, 0, 1, 0 ) 1098 1098 IF( pt_v(ji,jj,jl) < 0._wp .OR. ( jmsk_small(ji,jj,jl) == 0 .AND. pamsk == 0. ) ) THEN 1099 1099 pt_v(ji,jj,jl) = 0.5_wp * vmask(ji,jj,1) * ( ( pt(ji,jj+1,jl) + pt(ji,jj,jl) ) & … … 1105 1105 ! !-- High order flux in j-direction --! 1106 1106 DO jl = 1, jpl 1107 DO_2D _10_101107 DO_2D( 1, 0, 1, 0 ) 1108 1108 pfv_ho(ji,jj,jl) = pv(ji,jj) * pt_v(ji,jj,jl) 1109 1109 END_2D … … 1141 1141 ! -------------------------------------------------- 1142 1142 DO jl = 1, jpl 1143 DO_2D _10_101143 DO_2D( 1, 0, 1, 0 ) 1144 1144 pfu_ho(ji,jj,jl) = pfu_ho(ji,jj,jl) - pfu_ups(ji,jj,jl) 1145 1145 pfv_ho(ji,jj,jl) = pfv_ho(ji,jj,jl) - pfv_ups(ji,jj,jl) … … 1158 1158 1159 1159 DO jl = 1, jpl 1160 DO_2D _00_001160 DO_2D( 0, 0, 0, 0 ) 1161 1161 zti_ups(ji,jj,jl)= pt_ups(ji+1,jj ,jl) 1162 1162 ztj_ups(ji,jj,jl)= pt_ups(ji ,jj+1,jl) … … 1166 1166 1167 1167 DO jl = 1, jpl 1168 DO_2D _00_001168 DO_2D( 0, 0, 0, 0 ) 1169 1169 IF ( pfu_ho(ji,jj,jl) * ( pt_ups(ji+1,jj ,jl) - pt_ups(ji,jj,jl) ) <= 0._wp .AND. & 1170 1170 & pfv_ho(ji,jj,jl) * ( pt_ups(ji ,jj+1,jl) - pt_ups(ji,jj,jl) ) <= 0._wp ) THEN … … 1195 1195 DO jl = 1, jpl 1196 1196 1197 DO_2D _11_111197 DO_2D( 1, 1, 1, 1 ) 1198 1198 IF ( pt(ji,jj,jl) <= 0._wp .AND. pt_ups(ji,jj,jl) <= 0._wp ) THEN 1199 1199 zbup(ji,jj) = -zbig … … 1211 1211 END_2D 1212 1212 1213 DO_2D _00_001213 DO_2D( 0, 0, 0, 0 ) 1214 1214 ! 1215 1215 zup = MAX( zbup(ji,jj), zbup(ji-1,jj), zbup(ji+1,jj), zbup(ji,jj-1), zbup(ji,jj+1) ) ! search max/min in neighbourhood … … 1248 1248 ! --------------------------------- 1249 1249 DO jl = 1, jpl 1250 DO_2D _10_101250 DO_2D( 1, 0, 1, 0 ) 1251 1251 zau = MIN( 1._wp , zbetdo(ji,jj,jl) , zbetup(ji+1,jj,jl) ) 1252 1252 zbu = MIN( 1._wp , zbetup(ji,jj,jl) , zbetdo(ji+1,jj,jl) ) … … 1259 1259 END_2D 1260 1260 1261 DO_2D _10_101261 DO_2D( 1, 0, 1, 0 ) 1262 1262 zav = MIN( 1._wp , zbetdo(ji,jj,jl) , zbetup(ji,jj+1,jl) ) 1263 1263 zbv = MIN( 1._wp , zbetup(ji,jj,jl) , zbetdo(ji,jj+1,jl) ) … … 1293 1293 ! 1294 1294 DO jl = 1, jpl 1295 DO_2D _00_001295 DO_2D( 0, 0, 0, 0 ) 1296 1296 zslpx(ji,jj,jl) = ( pt(ji+1,jj,jl) - pt(ji,jj,jl) ) * umask(ji,jj,1) 1297 1297 END_2D … … 1300 1300 1301 1301 DO jl = 1, jpl 1302 DO_2D _00_001302 DO_2D( 0, 0, 0, 0 ) 1303 1303 uCFL = pdt * ABS( pu(ji,jj) ) * r1_e1e2t(ji,jj) 1304 1304 … … 1384 1384 ! 1385 1385 DO jl = 1, jpl 1386 DO_2D _00_001386 DO_2D( 0, 0, 0, 0 ) 1387 1387 zslpy(ji,jj,jl) = ( pt(ji,jj+1,jl) - pt(ji,jj,jl) ) * vmask(ji,jj,1) 1388 1388 END_2D … … 1391 1391 1392 1392 DO jl = 1, jpl 1393 DO_2D _00_001393 DO_2D( 0, 0, 0, 0 ) 1394 1394 vCFL = pdt * ABS( pv(ji,jj) ) * r1_e1e2t(ji,jj) 1395 1395 … … 1488 1488 ! 1489 1489 DO jl = 1, jpl 1490 DO_2D _11_111490 DO_2D( 1, 1, 1, 1 ) 1491 1491 IF ( pv_i(ji,jj,jl) > 0._wp ) THEN 1492 1492 ! … … 1535 1535 ! ! -- check e_i/v_i -- ! 1536 1536 DO jl = 1, jpl 1537 DO_3D _11_11(1, nlay_i )1537 DO_3D( 1, 1, 1, 1, 1, nlay_i ) 1538 1538 IF ( pv_i(ji,jj,jl) > 0._wp ) THEN 1539 1539 ! if e_i/v_i is larger than the surrounding 9 pts => put the heat excess in the ocean … … 1549 1549 ! ! -- check e_s/v_s -- ! 1550 1550 DO jl = 1, jpl 1551 DO_3D _11_11(1, nlay_s )1551 DO_3D( 1, 1, 1, 1, 1, nlay_s ) 1552 1552 IF ( pv_s(ji,jj,jl) > 0._wp ) THEN 1553 1553 ! if e_s/v_s is larger than the surrounding 9 pts => put the heat excess in the ocean … … 1592 1592 ! -- check snow load -- ! 1593 1593 DO jl = 1, jpl 1594 DO_2D _11_111594 DO_2D( 1, 1, 1, 1 ) 1595 1595 IF ( pv_i(ji,jj,jl) > 0._wp ) THEN 1596 1596 !
Note: See TracChangeset
for help on using the changeset viewer.