- Timestamp:
- 2017-12-07T16:28:21+01:00 (7 years ago)
- Location:
- branches/2017/dev_CNRS_2017/NEMOGCM/NEMO/OPA_SRC/SBC
- Files:
-
- 2 edited
Legend:
- Unmodified
- Added
- Removed
-
branches/2017/dev_CNRS_2017/NEMOGCM/NEMO/OPA_SRC/SBC/sbc_oce.F90
r8933 r8934 65 65 ! 66 66 INTEGER , PUBLIC :: nn_lsm !: Number of iteration if seaoverland is applied 67 ! 68 ! !!* namsbc_cpl namelist * 69 INTEGER , PUBLIC :: nn_cats_cpl !: Number of sea ice categories over which the coupling is carried out 70 67 71 !!---------------------------------------------------------------------- 68 72 !! switch definition (improve readability) -
branches/2017/dev_CNRS_2017/NEMOGCM/NEMO/OPA_SRC/SBC/sbccpl.F90
r8933 r8934 116 116 INTEGER, PARAMETER :: jpr_ts_ice = 54 ! Sea ice surface temp 117 117 118 INTEGER, PARAMETER :: jprcv = 5 5! total number of fields received118 INTEGER, PARAMETER :: jprcv = 54 ! total number of fields received 119 119 120 120 INTEGER, PARAMETER :: jps_fice = 1 ! ice fraction sent to the atmosphere … … 159 159 INTEGER, PARAMETER :: jpsnd = 38 ! total number of fields sent 160 160 161 INTEGER :: nn_cats_cpl ! number of sea ice categories over which the coupling is carried out162 163 161 ! !!** namelist namsbc_cpl ** 164 162 TYPE :: FLD_C ! … … 171 169 ! ! Send to the atmosphere 172 170 TYPE(FLD_C) :: sn_snd_temp , sn_snd_alb , sn_snd_thick, sn_snd_crt , sn_snd_co2, & 173 & sn_snd_thick1, sn_snd_cond, sn_snd_mpnd , sn_snd_sstfrz, sn_snd_ttilyr 171 & sn_snd_thick1, sn_snd_cond, sn_snd_mpnd , sn_snd_sstfrz, sn_snd_ttilyr 174 172 ! ! Received from the atmosphere 175 173 TYPE(FLD_C) :: sn_rcv_w10m, sn_rcv_taumod, sn_rcv_tau, sn_rcv_dqnsdt, sn_rcv_qsr, & 176 174 & sn_rcv_qns , sn_rcv_emp , sn_rcv_rnf, sn_rcv_ts_ice 177 TYPE(FLD_C) :: sn_rcv_cal, sn_rcv_iceflx, sn_rcv_co2, sn_rcv_mslp, sn_rcv_icb, sn_rcv_isf 175 TYPE(FLD_C) :: sn_rcv_cal, sn_rcv_iceflx, sn_rcv_co2, sn_rcv_mslp, sn_rcv_icb, sn_rcv_isf 178 176 ! Send to waves 179 177 TYPE(FLD_C) :: sn_snd_ifrac, sn_snd_crtw, sn_snd_wlev 180 178 ! Received from waves 181 TYPE(FLD_C) :: sn_rcv_hsig, sn_rcv_phioc,sn_rcv_sdrfx,sn_rcv_sdrfy,sn_rcv_wper,sn_rcv_wnum,sn_rcv_wstrf,sn_rcv_wdrag179 TYPE(FLD_C) :: sn_rcv_hsig, sn_rcv_phioc, sn_rcv_sdrfx, sn_rcv_sdrfy, sn_rcv_wper, sn_rcv_wnum, sn_rcv_wstrf, sn_rcv_wdrag 182 180 ! ! Other namelist parameters 183 181 INTEGER :: nn_cplmodel ! Maximum number of models to/from which NEMO is potentialy sending/receiving data 184 182 LOGICAL :: ln_usecplmask ! use a coupling mask file to merge data received from several models 185 ! -> file cplmask.nc with the float variable called cplmask (jpi,jpj,nn_cplmodel)183 ! -> file cplmask.nc with the float variable called cplmask (jpi,jpj,nn_cplmodel) 186 184 TYPE :: DYNARR 187 185 REAL(wp), POINTER, DIMENSION(:,:,:) :: z3 188 186 END TYPE DYNARR 189 187 190 TYPE( DYNARR ), SAVE, DIMENSION(jprcv) :: frcv 188 TYPE( DYNARR ), SAVE, DIMENSION(jprcv) :: frcv ! all fields recieved from the atmosphere 191 189 192 190 REAL(wp), ALLOCATABLE, SAVE, DIMENSION(:,:) :: alb_oce_mix ! ocean albedo sent to atmosphere (mix clear/overcast sky) … … 261 259 ! 262 260 CALL wrk_alloc( jpi,jpj, zacs, zaos ) 263 264 IF( ln_meto_cpl ) THEN265 tsfc_ice(:,:,:) = 0.0266 a_ip (:,:,:) = 0.0267 v_ip (:,:,:) = 0.0268 t1_ice (:,:,:) = rt0269 cnd_ice (:,:,:) = 0.0270 sstfrz (:,:) = 0.0271 ENDIF272 261 273 262 ! ================================ ! … … 459 448 SELECT CASE( TRIM( sn_rcv_emp%cldes ) ) 460 449 CASE( 'none' ) ! nothing to do 461 CASE( 'oce only' ) ; srcv( jpr_oemp)%laction = .TRUE.450 CASE( 'oce only' ) ; srcv(jpr_oemp)%laction = .TRUE. 462 451 CASE( 'conservative' ) 463 452 srcv( (/jpr_rain, jpr_snow, jpr_ievp, jpr_tevp/) )%laction = .TRUE. … … 503 492 CASE default ; CALL ctl_stop( 'sbc_cpl_init: wrong definition of sn_rcv_qns%cldes' ) 504 493 END SELECT 505 IF( TRIM( sn_rcv_qns%cldes ) == 'mixed oce-ice' .AND. jpl > 1 ) &494 IF( TRIM( sn_rcv_qns%cldes ) == 'mixed oce-ice' .AND. nn_cats_cpl > 1 ) & 506 495 CALL ctl_stop( 'sbc_cpl_init: sn_rcv_qns%cldes not currently allowed to be mixed oce-ice for multi-category ice' ) 496 ! 507 497 ! ! ------------------------- ! 508 498 ! ! solar radiation ! Qsr … … 519 509 CASE default ; CALL ctl_stop( 'sbc_cpl_init: wrong definition of sn_rcv_qsr%cldes' ) 520 510 END SELECT 521 IF( TRIM( sn_rcv_qsr%cldes ) == 'mixed oce-ice' .AND. jpl > 1 ) &511 IF( TRIM( sn_rcv_qsr%cldes ) == 'mixed oce-ice' .AND. nn_cats_cpl > 1 ) & 522 512 CALL ctl_stop( 'sbc_cpl_init: sn_rcv_qsr%cldes not currently allowed to be mixed oce-ice for multi-category ice' ) 513 ! 523 514 ! ! ------------------------- ! 524 515 ! ! non solar sensitivity ! d(Qns)/d(T) … … 534 525 535 526 ! non solar sensitivity mandatory for mixed oce-ice solar radiation coupling technique 536 IF( TRIM( sn_rcv_dqnsdt%cldes ) == 'none' .AND. TRIM( sn_rcv_qns%cldes ) == 'mixed oce-ice' ) & 537 CALL ctl_stop( 'sbc_cpl_init: namsbc_cpl namelist mismatch between sn_rcv_qns%cldes and sn_rcv_dqnsdt%cldes' ) 527 IF( TRIM( sn_rcv_dqnsdt%cldes ) == 'none' .AND. TRIM( sn_rcv_qns%cldes ) == 'mixed oce-ice' ) & 528 & CALL ctl_stop( 'sbc_cpl_init: namsbc_cpl namelist mismatch between sn_rcv_qns%cldes and sn_rcv_dqnsdt%cldes' ) 529 ! 538 530 ! ! ------------------------- ! 539 531 ! ! 10m wind module ! … … 546 538 srcv(jpr_taum)%clname = 'O_TauMod' ; IF( TRIM(sn_rcv_taumod%cldes) == 'coupled' ) srcv(jpr_taum)%laction = .TRUE. 547 539 lhftau = srcv(jpr_taum)%laction 548 540 ! 549 541 ! ! ------------------------- ! 550 542 ! ! Atmospheric CO2 ! … … 558 550 IF(lwp) WRITE(numout,*) 559 551 ENDIF 560 552 ! 561 553 ! ! ------------------------- ! 562 554 ! ! Mean Sea Level Pressure ! 563 555 ! ! ------------------------- ! 564 556 srcv(jpr_mslp)%clname = 'O_MSLP' ; IF( TRIM(sn_rcv_mslp%cldes ) == 'coupled' ) srcv(jpr_mslp)%laction = .TRUE. 565 566 ! ! ------------------------- ! 567 ! ! topmelt and botmelt!557 ! 558 ! ! ------------------------- ! 559 ! ! ice topmelt and botmelt ! 568 560 ! ! ------------------------- ! 569 561 srcv(jpr_topm )%clname = 'OTopMlt' … … 571 563 IF( TRIM(sn_rcv_iceflx%cldes) == 'coupled' ) THEN 572 564 IF ( TRIM( sn_rcv_iceflx%clcat ) == 'yes' ) THEN 573 srcv(jpr_topm:jpr_botm)%nct = jpl565 srcv(jpr_topm:jpr_botm)%nct = nn_cats_cpl 574 566 ELSE 575 567 CALL ctl_stop( 'sbc_cpl_init: sn_rcv_iceflx%clcat should always be set to yes currently' ) … … 577 569 srcv(jpr_topm:jpr_botm)%laction = .TRUE. 578 570 ENDIF 579 ! ! ------------------------- ----!580 581 ! !!!! To get NEMO4-LIM working at Met Office582 srcv(jpr_ts_ice)%clname = 'OTsfIce' 571 ! ! ------------------------- ! 572 ! ! ice skin temperature ! 573 ! ! ------------------------- ! 574 srcv(jpr_ts_ice)%clname = 'OTsfIce' ! needed by Met Office 583 575 IF ( TRIM( sn_rcv_ts_ice%cldes ) == 'ice' ) srcv(jpr_ts_ice)%laction = .TRUE. 584 576 IF ( TRIM( sn_rcv_ts_ice%clcat ) == 'yes' ) srcv(jpr_ts_ice)%nct = nn_cats_cpl 585 577 IF ( TRIM( sn_rcv_emp%clcat ) == 'yes' ) srcv(jpr_ievp)%nct = nn_cats_cpl 586 !!!!!587 578 588 579 ! ! ------------------------- ! … … 629 620 cpl_wdrag = .TRUE. 630 621 ENDIF 631 !632 622 ! ! ------------------------------- ! 633 623 ! ! OPA-SAS coupling - rcv by opa ! … … 755 745 CASE( 'oce and ice' , 'weighted oce and ice' , 'oce and weighted ice' ) 756 746 ssnd( (/jps_toce, jps_tice/) )%laction = .TRUE. 757 IF ( TRIM( sn_snd_temp%clcat ) == 'yes' ) ssnd(jps_tice)%nct = jpl747 IF ( TRIM( sn_snd_temp%clcat ) == 'yes' ) ssnd(jps_tice)%nct = nn_cats_cpl 758 748 CASE( 'mixed oce-ice' ) ; ssnd( jps_tmix )%laction = .TRUE. 759 749 CASE default ; CALL ctl_stop( 'sbc_cpl_init: wrong definition of sn_snd_temp%cldes' ) … … 780 770 alb_oce_mix(:,:) = ( zacs(:,:) + zaos(:,:) ) * 0.5 781 771 ENDIF 782 783 772 ! ! ------------------------- ! 784 773 ! ! Ice fraction & Thickness ! … … 792 781 ssnd(jps_fice1)%clname = 'OIceFrd' 793 782 IF( k_ice /= 0 ) THEN 794 ssnd(jps_fice)%laction = .TRUE. ! if ice treated in the ocean (even in climato case) 795 ssnd(jps_fice1)%laction = .TRUE. ! First-order regridded ice concentration, to be used 796 ! producing atmos-to-ice fluxes 783 ssnd(jps_fice)%laction = .TRUE. ! if ice treated in the ocean (even in climato case) 784 ssnd(jps_fice1)%laction = .TRUE. ! First-order regridded ice concentration, to be used producing atmos-to-ice fluxes (Met Office requirement) 797 785 ! Currently no namelist entry to determine sending of multi-category ice fraction so use the thickness entry for now 798 IF ( TRIM( sn_snd_thick%clcat ) == 'yes' ) ssnd(jps_fice)%nct = jpl799 IF ( TRIM( sn_snd_thick1%clcat ) == 'yes' ) ssnd(jps_fice1)%nct = jpl786 IF ( TRIM( sn_snd_thick%clcat ) == 'yes' ) ssnd(jps_fice)%nct = nn_cats_cpl 787 IF ( TRIM( sn_snd_thick1%clcat ) == 'yes' ) ssnd(jps_fice1)%nct = nn_cats_cpl 800 788 ENDIF 801 789 … … 807 795 ssnd(jps_hice:jps_hsnw)%laction = .TRUE. 808 796 IF ( TRIM( sn_snd_thick%clcat ) == 'yes' ) THEN 809 ssnd(jps_hice:jps_hsnw)%nct = jpl797 ssnd(jps_hice:jps_hsnw)%nct = nn_cats_cpl 810 798 ENDIF 811 799 CASE ( 'weighted ice and snow' ) 812 800 ssnd(jps_hice:jps_hsnw)%laction = .TRUE. 813 IF ( TRIM( sn_snd_thick%clcat ) == 'yes' ) ssnd(jps_hice:jps_hsnw)%nct = jpl801 IF ( TRIM( sn_snd_thick%clcat ) == 'yes' ) ssnd(jps_hice:jps_hsnw)%nct = nn_cats_cpl 814 802 CASE default ; CALL ctl_stop( 'sbc_cpl_init: wrong definition of sn_snd_thick%cldes' ) 815 803 END SELECT 816 804 817 805 ! ! ------------------------- ! 818 ! ! Ice Meltponds!806 ! ! Ice Meltponds ! 819 807 ! ! ------------------------- ! 820 821 822 !!!!! Getting NEMO4-LIM to work at Met Office 823 ssnd(jps_a_p)%clname = 'OPndFrc' 808 ! Needed by Met Office 809 ssnd(jps_a_p)%clname = 'OPndFrc' 824 810 ssnd(jps_ht_p)%clname = 'OPndTck' 825 811 SELECT CASE ( TRIM( sn_snd_mpnd%cldes ) ) 826 812 CASE ( 'none' ) 827 ssnd(jps_a_p)%laction = .FALSE.813 ssnd(jps_a_p)%laction = .FALSE. 828 814 ssnd(jps_ht_p)%laction = .FALSE. 829 815 CASE ( 'ice only' ) 830 ssnd(jps_a_p)%laction = .TRUE.816 ssnd(jps_a_p)%laction = .TRUE. 831 817 ssnd(jps_ht_p)%laction = .TRUE. 832 818 IF ( TRIM( sn_snd_mpnd%clcat ) == 'yes' ) THEN 833 ssnd(jps_a_p)%nct = nn_cats_cpl819 ssnd(jps_a_p)%nct = nn_cats_cpl 834 820 ssnd(jps_ht_p)%nct = nn_cats_cpl 835 821 ELSE 836 IF ( jpl > 1 ) THEN822 IF ( nn_cats_cpl > 1 ) THEN 837 823 CALL ctl_stop( 'sbc_cpl_init: use weighted ice option for sn_snd_mpnd%cldes if not exchanging category fields' ) 838 824 ENDIF 839 825 ENDIF 840 826 CASE ( 'weighted ice' ) 841 ssnd(jps_a_p)%laction = .TRUE.827 ssnd(jps_a_p)%laction = .TRUE. 842 828 ssnd(jps_ht_p)%laction = .TRUE. 843 829 IF ( TRIM( sn_snd_mpnd%clcat ) == 'yes' ) THEN 844 ssnd(jps_a_p)%nct = nn_cats_cpl830 ssnd(jps_a_p)%nct = nn_cats_cpl 845 831 ssnd(jps_ht_p)%nct = nn_cats_cpl 846 832 ENDIF 847 833 CASE default ; CALL ctl_stop( 'sbc_cpl_init: wrong definition of sn_snd_mpnd%cldes; '//sn_snd_mpnd%cldes ) 848 834 END SELECT 849 !!!!!850 835 851 836 ! ! ------------------------- ! … … 899 884 ssnd(jps_co2)%clname = 'O_CO2FLX' ; IF( TRIM(sn_snd_co2%cldes) == 'coupled' ) ssnd(jps_co2 )%laction = .TRUE. 900 885 ! 901 902 !!!!! Getting NEMO4-LIM to work at the Met Office903 886 ! ! ------------------------- ! 904 887 ! ! Sea surface freezing temp ! 905 888 ! ! ------------------------- ! 889 ! needed by Met Office 906 890 ssnd(jps_sstfrz)%clname = 'O_SSTFrz' ; IF( TRIM(sn_snd_sstfrz%cldes) == 'coupled' ) ssnd(jps_sstfrz)%laction = .TRUE. 907 !!!!!908 909 891 ! 910 892 ! ! ------------------------- ! 911 893 ! ! Ice conductivity ! 912 894 ! ! ------------------------- ! 895 ! needed by Met Office 913 896 ! Note that ultimately we will move to passing an ocean effective conductivity as well so there 914 897 ! will be some changes to the parts of the code which currently relate only to ice conductivity 915 916 898 ssnd(jps_ttilyr )%clname = 'O_TtiLyr' 917 899 SELECT CASE ( TRIM( sn_snd_ttilyr%cldes ) ) … … 952 934 END SELECT 953 935 ! 954 955 936 ! ! ------------------------- ! 956 937 ! ! Sea surface height ! … … 1119 1100 REAL(wp) :: zcdrag = 1.5e-3 ! drag coefficient 1120 1101 REAL(wp) :: zzx, zzy ! temporary variables 1121 REAL(wp), POINTER, DIMENSION(:,:) :: ztx, zty, zmsk, zemp, zqns, zqsr1102 REAL(wp), DIMENSION(jpi,jpj) :: ztx, zty, zmsk, zemp, zqns, zqsr 1122 1103 !!---------------------------------------------------------------------- 1123 1104 ! 1124 1105 IF( nn_timing == 1 ) CALL timing_start('sbc_cpl_rcv') 1125 !1126 CALL wrk_alloc( jpi,jpj, ztx, zty, zmsk, zemp, zqns, zqsr )1127 1106 ! 1128 1107 IF( ln_mixcpl ) zmsk(:,:) = 1. - xcplmask(:,:,0) … … 1196 1175 ! ! wind stress module ! (taum) 1197 1176 ! ! ========================= ! 1198 !1199 1177 IF( .NOT. srcv(jpr_taum)%laction ) THEN ! compute wind stress module from its components if not received 1200 1178 ! => need to be done only when otx1 was changed … … 1223 1201 ! ! 10 m wind speed ! (wndm) 1224 1202 ! ! ========================= ! 1225 !1226 1203 IF( .NOT. srcv(jpr_w10m)%laction ) THEN ! compute wind spreed from wind stress module if not received 1227 1204 ! => need to be done only when taumod was changed … … 1259 1236 ! ! ================== ! 1260 1237 IF( srcv(jpr_co2)%laction ) atm_co2(:,:) = frcv(jpr_co2)%z3(:,:,1) 1261 ! 1262 1263 !!!!! Getting NEMO4-LIM to work at the Met Office 1264 ! ! Sea ice surface skin temp: 1238 ! 1239 ! ! ================== ! 1240 ! ! ice skin temp. ! 1241 ! ! ================== ! 1242 #if defined key_lim3 1243 ! needed by Met Office 1265 1244 IF( srcv(jpr_ts_ice)%laction ) THEN 1266 1245 DO jn = 1, jpl … … 1278 1257 END DO 1279 1258 ENDIF 1280 !!!!! 1281 1282 1259 #endif 1283 1260 ! ! ========================= ! 1284 1261 ! ! Mean Sea Level Pressure ! (taum) 1285 1262 ! ! ========================= ! 1286 !1287 1263 IF( srcv(jpr_mslp)%laction ) THEN ! UKMO SHELF effect of atmospheric pressure on SSH 1288 1264 IF( kt /= nit000 ) ssh_ibb(:,:) = ssh_ib(:,:) !* Swap of ssh_ib fields … … 1460 1436 ! 1461 1437 ENDIF 1462 !1463 CALL wrk_dealloc( jpi,jpj, ztx, zty, zmsk, zemp, zqns, zqsr )1464 1438 ! 1465 1439 IF( nn_timing == 1 ) CALL timing_stop('sbc_cpl_rcv') … … 1506 1480 INTEGER :: ji, jj ! dummy loop indices 1507 1481 INTEGER :: itx ! index of taux over ice 1508 REAL(wp), POINTER, DIMENSION(:,:) :: ztx, zty1482 REAL(wp), DIMENSION(jpi,jpj) :: ztx, zty 1509 1483 !!---------------------------------------------------------------------- 1510 1484 ! 1511 1485 IF( nn_timing == 1 ) CALL timing_start('sbc_cpl_ice_tau') 1512 1486 ! 1513 CALL wrk_alloc( jpi,jpj, ztx, zty )1514 1515 1487 IF( srcv(jpr_itx1)%laction ) THEN ; itx = jpr_itx1 1516 1488 ELSE ; itx = jpr_otx1 … … 1671 1643 ENDIF 1672 1644 ! 1673 CALL wrk_dealloc( jpi,jpj, ztx, zty )1674 !1675 1645 IF( nn_timing == 1 ) CALL timing_stop('sbc_cpl_ice_tau') 1676 1646 ! … … 2093 2063 dqns_ice(:,:,:) = zdqns_ice(:,:,:) 2094 2064 ENDIF 2095 2065 2066 #if defined key_lim3 2096 2067 IF( ln_meto_cpl ) THEN 2097 2068 ! ! ========================= ! … … 2104 2075 ENDIF 2105 2076 2106 2107 #if defined key_lim32108 2077 ! ! ========================= ! 2109 2078 ! ! Transmitted Qsr ! [W/m2] … … 2147 2116 INTEGER :: isec, info ! local integer 2148 2117 REAL(wp) :: zumax, zvmax 2149 REAL(wp), POINTER, DIMENSION(:,:):: zfr_l, ztmp1, ztmp2, zotx1, zoty1, zotz1, zitx1, zity1, zitz12150 REAL(wp), POINTER, DIMENSION(:,:,:) :: ztmp3, ztmp42118 REAL(wp), DIMENSION(jpi,jpj) :: zfr_l, ztmp1, ztmp2, zotx1, zoty1, zotz1, zitx1, zity1, zitz1 2119 REAL(wp), DIMENSION(jpi,jpj,jpl) :: ztmp3, ztmp4 2151 2120 !!---------------------------------------------------------------------- 2152 2121 ! 2153 2122 IF( nn_timing == 1 ) CALL timing_start('sbc_cpl_snd') 2154 2123 ! 2155 CALL wrk_alloc( jpi,jpj, zfr_l, ztmp1, ztmp2, zotx1, zoty1, zotz1, zitx1, zity1, zitz1 )2156 CALL wrk_alloc( jpi,jpj,jpl, ztmp3, ztmp4 )2157 2158 2124 isec = ( kt - nit000 ) * NINT( rdt ) ! date of exchanges 2159 2125 … … 2221 2187 ENDIF 2222 2188 2189 #if defined key_lim3 2223 2190 !!!!! Getting NEMO4-LIM working at Met Office 2224 2191 ! Top layer ice temperature … … 2232 2199 ENDIF 2233 2200 !!!!! 2234 2201 #endif 2235 2202 2236 2203 ! ! ------------------------- ! … … 2347 2314 ENDIF 2348 2315 2316 #if defined key_lim3 2349 2317 ! NEMO4 - Jules coupling - Met Office 2350 2318 ! Send meltpond fields … … 2390 2358 ! 2391 2359 !!!!! 2392 2360 #endif 2393 2361 2394 2362 ! ! ------------------------- ! … … 2713 2681 IF( ssnd(jps_taum )%laction ) CALL cpl_snd( jps_taum , isec, RESHAPE ( taum, (/jpi,jpj,1/) ), info ) 2714 2682 2683 #if defined key_lim3 2715 2684 ! NEMO4 - Jules coupling - Met Office 2716 2685 CALL eos_fzp(tsn(:,:,1,jp_sal), sstfrz) 2717 2686 ztmp1(:,:) = sstfrz(:,:) + rt0 2718 2687 IF( ssnd(jps_sstfrz)%laction ) CALL cpl_snd( jps_sstfrz, isec, RESHAPE ( ztmp1, (/jpi,jpj,1/) ), info) 2719 2720 CALL wrk_dealloc( jpi,jpj, zfr_l, ztmp1, ztmp2, zotx1, zoty1, zotz1, zitx1, zity1, zitz1 ) 2721 CALL wrk_dealloc( jpi,jpj,jpl, ztmp3, ztmp4 ) 2688 #endif 2722 2689 ! 2723 2690 IF( nn_timing == 1 ) CALL timing_stop('sbc_cpl_snd')
Note: See TracChangeset
for help on using the changeset viewer.