Changeset 590 for LMDZ4/trunk/libf
- Timestamp:
- Feb 7, 2005, 5:41:35 PM (20 years ago)
- Location:
- LMDZ4/trunk/libf/phylmd
- Files:
-
- 6 edited
-
clmain.F (modified) (7 diffs)
-
ini_histhf.h (modified) (1 diff)
-
interface_surf.F90 (modified) (27 diffs)
-
oasis.F (modified) (1 diff)
-
oasis.dummy (modified) (9 diffs)
-
oasis.true (modified) (11 diffs)
Legend:
- Unmodified
- Added
- Removed
-
LMDZ4/trunk/libf/phylmd/clmain.F
r541 r590 241 241 c 242 242 REAL yt2m(klon), yq2m(klon), yu10m(klon) 243 c -- LOOP 244 REAL yu10mx(klon) 245 REAL yu10my(klon) 246 REAL ywindsp(klon) 247 c -- LOOP 243 248 c 244 249 REAL uzon(klon), vmer(klon) … … 351 356 ytsoil = 999999. 352 357 yrugoro = 0. 353 358 c -- LOOP 359 yu10mx = 0.0 360 yu10my = 0.0 361 ywindsp = 0.0 362 c -- LOOP 354 363 DO nsrf = 1, nbsrf 355 364 DO i = 1, klon … … 458 467 ypaprs(j,klev+1) = paprs(i,klev+1) 459 468 y_run_off_lic_0(j) = run_off_lic_0(i) 469 c -- LOOP 470 yu10mx(j) = u10m(i,nsrf) 471 yu10my(j) = v10m(i,nsrf) 472 ywindsp(j) = SQRT(yu10mx(j)*yu10mx(j) + yu10my(j)*yu10my(j) ) 473 c -- LOOP 460 474 END DO 461 475 C … … 667 681 e yt,yq,yts,ypaprs,ypplay, 668 682 e ydelp,yrads,yalb, yalblw, ysnow, yqsurf, 669 e yrain_f, ysnow_f, yfder, ytaux, ytauy, 683 e yrain_f, ysnow_f, yfder, ytaux, ytauy, 684 c -- LOOP 685 e ywindsp, 686 c -- LOOP 670 687 c$$$ e ysollw, ysolsw, 671 688 e ysollw, ysollwdown, ysolsw,yfluxlat, … … 900 917 e delp,radsol,albedo,alblw,snow,qsurf, 901 918 e precip_rain, precip_snow, fder, taux, tauy, 919 c -- LOOP 920 e ywindsp, 921 c -- LOOP 902 922 $ sollw, sollwdown, swnet,fluxlat, 903 923 s pctsrf_new, agesno, … … 958 978 character*6 ocean 959 979 integer npas, nexca 980 c -- LOOP 981 REAL yu10mx(klon) 982 REAL yu10my(klon) 983 REAL ywindsp(klon) 984 c -- LOOP 985 960 986 961 987 c … … 1184 1210 e tq_cdrag, petAcoef, peqAcoef, petBcoef, peqBcoef, 1185 1211 e precip_rain, precip_snow, sollw, sollwdown, swnet, swdown, 1186 e fder, taux, tauy, rugos, rugoro, 1212 e fder, taux, tauy, 1213 c -- LOOP 1214 e ywindsp, 1215 c -- LOOP 1216 e rugos, rugoro, 1187 1217 e albedo, snow, qsurf, 1188 1218 e ts, p1lay, psref, radsol, -
LMDZ4/trunk/libf/phylmd/ini_histhf.h
r524 r590 113 113 . "ave(X)", zsto1,zout) 114 114 c 115 c 115 c -- LOOP 116 116 CALL histdef(nid_hf, "SWdownOR", 117 . "Sfce incident SW radiation OR", "W/m^2", 118 . iim,jjmp1,nhori, 1,1,1, -99, 32, 119 . "ave(X)", zsto1,zout) 117 . "Sfce incident SW down radiation OR", "W/m^2", 118 . iim,jjmp1,nhori, 1,1,1, -99, 32, 119 . "ave(X)", zsto1,zout) 120 c 121 CALL histdef(nid_hf, "SWnetOR", 122 . "Sfce incident SW net radiation OR", "W/m^2", 123 . iim,jjmp1,nhori, 1,1,1, -99, 32, 124 . "ave(X)", zsto1,zout) 125 c -- LOOP 120 126 c 121 127 CALL histdef(nid_hf, "LWdownOR", -
LMDZ4/trunk/libf/phylmd/interface_surf.F90
r560 r590 64 64 & tq_cdrag, petAcoef, peqAcoef, petBcoef, peqBcoef, & 65 65 & precip_rain, precip_snow, sollw, sollwdown, swnet, swdown, & 66 & fder, taux, tauy, rugos, rugoro, & 66 & fder, taux, tauy, & 67 ! -- LOOP 68 & windsp, & 69 ! -- LOOP 70 & rugos, rugoro, & 67 71 & albedo, snow, qsurf, & 68 72 & tsurf, p1lay, ps, radsol, & … … 127 131 ! fder derivee des flux (pour le couplage) 128 132 ! taux, tauy tension de vents 133 ! -- LOOP 134 ! windsp module du vent a 10m 135 ! -- LOOP 129 136 ! rugos rugosite 130 137 ! zmasq masque terre/ocean … … 175 182 real, dimension(klon), intent(IN) :: zmasq 176 183 real, dimension(klon), intent(IN) :: taux, tauy, rugos, rugoro 184 ! -- LOOP 185 real, dimension(klon), intent(IN) :: windsp 186 ! -- LOOP 177 187 character (len = 6) :: ocean 178 188 integer :: npas, nexca ! nombre et pas de temps couplage … … 444 454 & ocean, npas, nexca, debut, lafin, & 445 455 & swdown, sollw, precip_rain, precip_snow, evap, tsurf, & 446 & fluxlat, fluxsens, fder, albedo, taux, tauy, zmasq, & 447 & tsurf_new, alb_new, pctsrf_new) 456 & fluxlat, fluxsens, fder, albedo, taux, tauy, & 457 ! -- LOOP 458 & windsp, & 459 ! -- LOOP 460 & zmasq, & 461 & tsurf_new, alb_new, & 462 & pctsrf_new) 448 463 449 464 ! else if (ocean == 'slab ') then … … 498 513 & ocean, npas, nexca, debut, lafin, & 499 514 & swdown, sollw, precip_rain, precip_snow, evap, tsurf, & 500 & fluxlat, fluxsens, fder, albedo, taux, tauy, zmasq, & 501 & tsurf_new, alb_new, pctsrf_new) 515 & fluxlat, fluxsens, fder, albedo, taux, tauy, & 516 ! -- LOOP 517 & windsp, & 518 ! -- LOOP 519 & zmasq, & 520 & tsurf_new, alb_new, & 521 & pctsrf_new) 502 522 503 523 ! else if (ocean == 'slab ') then … … 549 569 & ocean, npas, nexca, debut, lafin, & 550 570 & swdown, sollw, precip_rain, precip_snow, evap, tsurf, & 551 & fluxlat, fluxsens, fder, albedo, taux, tauy, zmasq, & 552 & tsurf_new, alb_new, pctsrf_new) 571 & fluxlat, fluxsens, fder, albedo, taux, tauy, & 572 ! -- LOOP 573 & windsp, & 574 ! -- LOOP 575 & zmasq, & 576 & tsurf_new, alb_new, & 577 & pctsrf_new) 553 578 554 579 tsurf_temp = tsurf_new … … 644 669 & ocean, npas, nexca, debut, lafin, & 645 670 & swdown, sollw, precip_rain, precip_snow, evap, tsurf, & 646 & fluxlat, fluxsens, fder, albedo, taux, tauy, zmasq, & 647 & tsurf_new, alb_new, pctsrf_new) 671 & fluxlat, fluxsens, fder, albedo, taux, tauy, & 672 ! -- LOOP 673 & windsp, & 674 ! -- LOOP 675 & zmasq, & 676 & tsurf_new, alb_new, & 677 & pctsrf_new) 648 678 649 679 ! else if (ocean == 'slab ') then … … 1198 1228 & ocean, npas, nexca, debut, lafin, & 1199 1229 & swdown, lwdown, precip_rain, precip_snow, evap, tsurf, & 1200 & fluxlat, fluxsens, fder, albsol, taux, tauy, zmasq, & 1201 & tsurf_new, alb_new, pctsrf_new) 1230 & fluxlat, fluxsens, fder, albsol, taux, tauy, & 1231 ! -- LOOP 1232 & windsp, & 1233 ! -- LOOP 1234 & zmasq, & 1235 & tsurf_new, alb_new, & 1236 & pctsrf_new) 1202 1237 1203 1238 ! Cette routine sert d'interface entre le modele atmospherique et un … … 1242 1277 ! taux tension de vent en x 1243 1278 ! tauy tension de vent en y 1279 ! -- LOOP 1280 ! windsp module du vent a 10m 1281 ! -- LOOP 1244 1282 ! nexca frequence de couplage 1245 1283 ! zmasq masque terre/ocean … … 1269 1307 real, dimension(klon), intent(IN) :: precip_rain, precip_snow 1270 1308 real, dimension(klon), intent(IN) :: tsurf, fder, albsol, taux, tauy 1309 ! -- LOOP 1310 real, dimension(klon), intent(IN) :: windsp 1311 ! -- LOOP 1271 1312 INTEGER :: nexca, npas, kstep 1272 1313 real, dimension(klon), intent(IN) :: zmasq … … 1281 1322 ! Variables locales 1282 1323 integer :: j, error, sum_error, ig, cpl_index,i 1324 ! -- LOOP 1325 INTEGER :: nsrf 1326 ! -- LOOP 1283 1327 character (len = 20) :: modname = 'interfoce_cpl' 1284 1328 character (len = 80) :: abort_message … … 1288 1332 real, allocatable, dimension(:,:),save :: cpl_snow, cpl_evap, cpl_tsol 1289 1333 real, allocatable, dimension(:,:),save :: cpl_fder, cpl_albe, cpl_taux 1334 ! -- LOOP 1335 real, allocatable, dimension(:,:),save :: cpl_windsp 1336 ! -- LOOP 1290 1337 real, allocatable, dimension(:,:),save :: cpl_tauy 1291 1338 REAL, ALLOCATABLE, DIMENSION(:,:),SAVE :: cpl_rriv, cpl_rcoa, cpl_rlic … … 1295 1342 real, allocatable, dimension(:,:,:),save :: tmp_snow, tmp_evap, tmp_tsol 1296 1343 real, allocatable, dimension(:,:,:),save :: tmp_fder, tmp_albe, tmp_taux 1344 ! -- LOOP 1345 real, allocatable, dimension(:,:,:),save :: tmp_windsp 1346 ! -- LOOP 1297 1347 !!$ real, allocatable, dimension(:,:,:),save :: tmp_tauy, tmp_rriv, tmp_rcoa 1298 1348 REAL, ALLOCATABLE, DIMENSION(:,:,:),SAVE :: tmp_tauy … … 1302 1352 REAL, DIMENSION(iim, jjm+1) :: wri_evap_sea, wri_rcoa, wri_rriv 1303 1353 REAL, DIMENSION(iim, jjm+1) :: wri_rain, wri_snow, wri_taux, wri_tauy 1354 ! -- LOOP 1355 REAL, DIMENSION(iim, jjm+1) :: wri_windsp 1356 ! -- LOOP 1304 1357 REAL, DIMENSION(iim, jjm+1) :: wri_calv 1305 1358 REAL, DIMENSION(iim, jjm+1) :: wri_tauxx, wri_tauyy, wri_tauzz … … 1334 1387 !med integer :: itau_w 1335 1388 integer,save :: itau_w 1389 ! -- LOOP 1390 integer :: nb_interf_cpl 1391 ! -- LOOP 1336 1392 #include "param_cou.h" 1337 1393 #include "inc_cpl.h" … … 1369 1425 allocate(cpl_albe(klon,2), stat = error); sum_error = sum_error + error 1370 1426 allocate(cpl_taux(klon,2), stat = error); sum_error = sum_error + error 1427 ! -- LOOP 1428 allocate(cpl_windsp(klon,2), stat = error); sum_error = sum_error + error 1429 ! -- LOOP 1371 1430 allocate(cpl_tauy(klon,2), stat = error); sum_error = sum_error + error 1372 1431 ALLOCATE(cpl_rriv(iim,jjm+1), stat=error); sum_error = sum_error + error … … 1386 1445 cpl_evap = 0.; cpl_tsol = 0.; cpl_fder = 0.; cpl_albe = 0. 1387 1446 cpl_taux = 0.; cpl_tauy = 0.; cpl_rriv = 0.; cpl_rcoa = 0.; cpl_rlic = 0. 1447 ! -- LOOP 1448 cpl_windsp = 0. 1449 ! -- LOOP 1388 1450 1389 1451 sum_error = 0 … … 1465 1527 1466 1528 ! calcul des fluxs a passer 1467 1529 ! -- LOOP 1530 nb_interf_cpl = nb_interf_cpl + 1 1531 if (check) write(lunout,*)'passage dans interface_surf.F90 : ',nb_interf_cpl 1532 ! -- LOOP 1468 1533 cpl_index = 1 1469 1534 if (nisurf == is_sic) cpl_index = 2 1470 1535 if (cumul) then 1536 ! -- LOOP 1537 if (check) write(lunout,*)'passage dans cumul ' 1538 if (check) write(lunout,*)'valeur de cpl_index ', cpl_index 1539 ! -- LOOP 1471 1540 if (check) write(*,*) modname, 'cumul des champs' 1472 1541 do ig = 1, knon … … 1492 1561 cpl_tauy(ig,cpl_index) = cpl_tauy(ig,cpl_index) & 1493 1562 & + tauy(ig) / FLOAT(nexca) 1563 ! -- LOOP 1564 IF (cpl_index .EQ. 1) THEN 1565 cpl_windsp(ig,cpl_index) = cpl_windsp(ig,cpl_index) & 1566 & + windsp(ig) / FLOAT(nexca) 1567 ENDIF 1568 ! -- LOOP 1494 1569 enddo 1495 1570 IF (cpl_index .EQ. 1) THEN … … 1603 1678 allocate(tmp_taux(iim,jjm+1,2), stat=error); sum_error = sum_error + error 1604 1679 allocate(tmp_tauy(iim,jjm+1,2), stat=error); sum_error = sum_error + error 1680 ! -- LOOP 1681 allocate(tmp_windsp(iim,jjm+1,2), stat=error); sum_error = sum_error + error 1682 ! -- LOOP 1605 1683 !!$ allocate(tmp_rriv(iim,jjm+1,2), stat=error); sum_error = sum_error + error 1606 1684 !!$ allocate(tmp_rcoa(iim,jjm+1,2), stat=error); sum_error = sum_error + error … … 1625 1703 call gath2cpl(cpl_albe(1,cpl_index), tmp_albe(1,1,cpl_index), klon, knon,iim,jjm, knindex) 1626 1704 call gath2cpl(cpl_taux(1,cpl_index), tmp_taux(1,1,cpl_index), klon, knon,iim,jjm, knindex) 1705 ! -- LOOP 1706 call gath2cpl(cpl_windsp(1,cpl_index), tmp_windsp(1,1,cpl_index), klon, knon,iim,jjm, knindex) 1707 ! -- LOOP 1627 1708 call gath2cpl(cpl_tauy(1,cpl_index), tmp_tauy(1,1,cpl_index), klon, knon,iim,jjm, knindex) 1628 1709 … … 1633 1714 wri_rain = 0.; wri_snow = 0.; wri_rcoa = 0.; wri_rriv = 0. 1634 1715 wri_taux = 0.; wri_tauy = 0. 1716 ! -- LOOP 1717 wri_windsp = 0. 1718 ! -- LOOP 1635 1719 call gath2cpl(pctsrf(1,is_oce), tamp_srf(1,1,1), klon, klon, iim, jjm, tamp_ind) 1636 1720 call gath2cpl(pctsrf(1,is_sic), tamp_srf(1,1,2), klon, klon, iim, jjm, tamp_ind) … … 1643 1727 wri_evap_ice = tmp_evap(:,:,2) 1644 1728 wri_evap_sea = tmp_evap(:,:,1) 1729 ! -- LOOP 1730 wri_windsp = tmp_windsp(:,:,1) 1731 ! -- LOOP 1732 1645 1733 !!$PB 1646 1734 wri_rriv = cpl_rriv(:,:) … … 1696 1784 ! envoi au coupleur 1697 1785 ! 1698 CALL histwrite(nidct,cl_writ(1),itau_w,wri_sol_ice,iim*(jjm+1),ndexct) 1699 CALL histwrite(nidct,cl_writ(2),itau_w,wri_sol_sea,iim*(jjm+1),ndexct) 1700 CALL histwrite(nidct,cl_writ(3),itau_w,wri_nsol_ice,iim*(jjm+1),ndexct) 1701 CALL histwrite(nidct,cl_writ(4),itau_w,wri_nsol_sea,iim*(jjm+1),ndexct) 1702 CALL histwrite(nidct,cl_writ(5),itau_w,wri_fder_ice,iim*(jjm+1),ndexct) 1703 CALL histwrite(nidct,cl_writ(6),itau_w,wri_evap_ice,iim*(jjm+1),ndexct) 1704 CALL histwrite(nidct,cl_writ(7),itau_w,wri_evap_sea,iim*(jjm+1),ndexct) 1705 CALL histwrite(nidct,cl_writ(8),itau_w,wri_rain,iim*(jjm+1),ndexct) 1706 CALL histwrite(nidct,cl_writ(9),itau_w,wri_snow,iim*(jjm+1),ndexct) 1707 CALL histwrite(nidct,cl_writ(10),itau_w,wri_rcoa,iim*(jjm+1),ndexct) 1708 CALL histwrite(nidct,cl_writ(11),itau_w,wri_rriv,iim*(jjm+1),ndexct) 1709 CALL histwrite(nidct,cl_writ(12),itau_w,wri_calv,iim*(jjm+1),ndexct) 1710 CALL histwrite(nidct,cl_writ(13),itau_w,wri_tauxx,iim*(jjm+1),ndexct) 1711 CALL histwrite(nidct,cl_writ(14),itau_w,wri_tauyy,iim*(jjm+1),ndexct) 1712 CALL histwrite(nidct,cl_writ(15),itau_w,wri_tauzz,iim*(jjm+1),ndexct) 1713 CALL histwrite(nidct,cl_writ(16),itau_w,wri_tauxx,iim*(jjm+1),ndexct) 1714 CALL histwrite(nidct,cl_writ(17),itau_w,wri_tauyy,iim*(jjm+1),ndexct) 1715 CALL histwrite(nidct,cl_writ(18),itau_w,wri_tauzz,iim*(jjm+1),ndexct) 1786 CALL histwrite(nidct,cl_writ(8),itau_w,wri_sol_ice,iim*(jjm+1),ndexct) 1787 CALL histwrite(nidct,cl_writ(9),itau_w,wri_sol_sea,iim*(jjm+1),ndexct) 1788 CALL histwrite(nidct,cl_writ(10),itau_w,wri_nsol_ice,iim*(jjm+1),ndexct) 1789 CALL histwrite(nidct,cl_writ(11),itau_w,wri_nsol_sea,iim*(jjm+1),ndexct) 1790 CALL histwrite(nidct,cl_writ(12),itau_w,wri_fder_ice,iim*(jjm+1),ndexct) 1791 CALL histwrite(nidct,cl_writ(13),itau_w,wri_evap_ice,iim*(jjm+1),ndexct) 1792 CALL histwrite(nidct,cl_writ(14),itau_w,wri_evap_sea,iim*(jjm+1),ndexct) 1793 CALL histwrite(nidct,cl_writ(15),itau_w,wri_rain,iim*(jjm+1),ndexct) 1794 CALL histwrite(nidct,cl_writ(16),itau_w,wri_snow,iim*(jjm+1),ndexct) 1795 CALL histwrite(nidct,cl_writ(17),itau_w,wri_rcoa,iim*(jjm+1),ndexct) 1796 CALL histwrite(nidct,cl_writ(18),itau_w,wri_rriv,iim*(jjm+1),ndexct) 1797 CALL histwrite(nidct,cl_writ(19),itau_w,wri_calv,iim*(jjm+1),ndexct) 1798 CALL histwrite(nidct,cl_writ(1),itau_w,wri_tauxx,iim*(jjm+1),ndexct) 1799 CALL histwrite(nidct,cl_writ(2),itau_w,wri_tauyy,iim*(jjm+1),ndexct) 1800 CALL histwrite(nidct,cl_writ(3),itau_w,wri_tauzz,iim*(jjm+1),ndexct) 1801 CALL histwrite(nidct,cl_writ(4),itau_w,wri_tauxx,iim*(jjm+1),ndexct) 1802 CALL histwrite(nidct,cl_writ(5),itau_w,wri_tauyy,iim*(jjm+1),ndexct) 1803 CALL histwrite(nidct,cl_writ(6),itau_w,wri_tauzz,iim*(jjm+1),ndexct) 1804 ! -- LOOP 1805 CALL histwrite(nidct,cl_writ(7),itau_w,wri_windsp,iim*(jjm+1),ndexct) 1806 ! -- LOOP 1716 1807 CALL histsync(nidct) 1717 1808 ! pas utile IF (lafin) CALL histclo(nidct) … … 1723 1814 & wri_nsol_sea, wri_fder_ice, wri_evap_ice, wri_evap_sea, wri_rain, & 1724 1815 & wri_snow, wri_rcoa, wri_rriv, wri_calv, wri_tauxx, wri_tauyy, & 1725 & wri_tauzz, wri_tauxx, wri_tauyy, wri_tauzz,lafin ) 1816 & wri_tauzz, wri_tauxx, wri_tauyy, wri_tauzz, & 1817 ! -- LOOP 1818 & wri_windsp,lafin) 1819 ! -- LOOP 1726 1820 #else 1727 1821 call intocpl(itime, (jjm+1)*iim, wri_sol_ice, wri_sol_sea, wri_nsol_ice,& 1728 1822 & wri_nsol_sea, wri_fder_ice, wri_evap_ice, wri_evap_sea, wri_rain, & 1729 1823 & wri_snow, wri_rcoa, wri_rriv, wri_calv, wri_tauxx, wri_tauyy, & 1730 & wri_tauzz, wri_tauxx, wri_tauyy, wri_tauzz,lafin ) 1824 & wri_tauzz, wri_tauxx, wri_tauyy, wri_tauzz, & 1825 ! -- LOOP 1826 & wri_windsp,lafin) 1827 ! -- LOOP 1731 1828 #endif 1732 1829 #endif … … 1735 1832 cpl_evap = 0.; cpl_tsol = 0.; cpl_fder = 0.; cpl_albe = 0. 1736 1833 cpl_taux = 0.; cpl_tauy = 0.; cpl_rriv = 0.; cpl_rcoa = 0.; cpl_rlic = 0. 1834 ! -- LOOP 1835 cpl_windsp = 0. 1836 ! -- LOOP 1737 1837 ! 1738 1838 ! deallocation memoire variables temporaires … … 1749 1849 deallocate(tmp_taux, stat=error); sum_error = sum_error + error 1750 1850 deallocate(tmp_tauy, stat=error); sum_error = sum_error + error 1851 ! -- LOOP 1852 deallocate(tmp_windsp, stat=error); sum_error = sum_error + error 1853 ! -- LOOP 1751 1854 !!$PB 1752 1855 !!$ deallocate(tmp_rriv, stat=error); sum_error = sum_error + error -
LMDZ4/trunk/libf/phylmd/oasis.F
r541 r590 2 2 ! $Header$ 3 3 ! 4 C $Id$ 4 5 C**** 5 6 C -
LMDZ4/trunk/libf/phylmd/oasis.dummy
r524 r590 20 20 c INCLUDE "param.h" 21 21 c 22 INTEGER kastp, kexch, kstep,imjm 22 INTEGER kastp, kexch, kstep,imjm,klon 23 23 INTEGER iparal(3) 24 24 INTEGER ifcpl, idt, info, imxtag, istep, jf 25 26 c -- LOOP 27 c 28 #include "dimensions.h" 29 INTEGER jjmp1 30 PARAMETER (jjmp1=jjm+1-1/jjm) 31 #include "dimphy.h" 32 REAL zwindsp(klon) 33 c 34 c -- LOOP 25 35 c 26 36 #include "param_cou.h" … … 72 82 c must be the same as (1) of the field definition in namcouple: 73 83 c 74 cl_writ( 1)='COSHFICE'75 cl_writ( 2)='COSHFOCE'76 cl_writ( 3)='CONSFICE'77 cl_writ( 4)='CONSFOCE'78 cl_writ( 5)='CODFLXDT'84 cl_writ(8)='COSHFICE' 85 cl_writ(9)='COSHFOCE' 86 cl_writ(10)='CONSFICE' 87 cl_writ(11)='CONSFOCE' 88 cl_writ(12)='CODFLXDT' 79 89 c cl_writ(6)='COICTEMP' 80 cl_writ( 6)='COTFSICE'81 cl_writ( 7)='COTFSOCE'82 cl_writ( 8)='COTOLPSU'83 cl_writ( 9)='COTOSPSU'84 cl_writ(1 0)='CORUNCOA'85 cl_writ(1 1)='CORIVFLU'86 cl_writ(1 2)='COCALVIN'90 cl_writ(13)='COTFSICE' 91 cl_writ(14)='COTFSOCE' 92 cl_writ(15)='COTOLPSU' 93 cl_writ(16)='COTOSPSU' 94 cl_writ(17)='CORUNCOA' 95 cl_writ(18)='CORIVFLU' 96 cl_writ(19)='COCALVIN' 87 97 c$$$ cl_writ(13)='COZOTAUX' 88 98 c$$$ cl_writ(14)='COZOTAUV' 89 99 c$$$ cl_writ(15)='COMETAUY' 90 100 c$$$ cl_writ(16)='COMETAUU' 91 cl_writ(13)='COTAUXXU' 92 cl_writ(14)='COTAUYYU' 93 cl_writ(15)='COTAUZZU' 94 cl_writ(16)='COTAUXXV' 95 cl_writ(17)='COTAUYYV' 96 cl_writ(18)='COTAUZZV' 101 cl_writ(1)='COTAUXXU' 102 cl_writ(2)='COTAUYYU' 103 cl_writ(3)='COTAUZZU' 104 cl_writ(4)='COTAUXXV' 105 cl_writ(5)='COTAUYYV' 106 cl_writ(6)='COTAUZZV' 107 c -- LOOP 108 cl_writ(7)='COWINDSP' 109 c -- LOOP 97 110 c 98 111 c Define files name for fields exchanged from atmos to coupler, … … 117 130 cl_f_writ(17)='flxatmos' 118 131 cl_f_writ(18)='flxatmos' 132 c -- LOOP 133 cl_f_writ(19)='flxatmos' 134 c -- LOOP 119 135 120 136 c … … 310 326 311 327 c $Id$ 328 c -- LOOP 312 329 SUBROUTINE intocpl(kt, imjm, fsolice, fsolwat, fnsolice, fnsolwat, 313 330 $ fnsicedt, evice, evwat, lpre, spre, dirunoff, rivrunoff, 314 $ calving, tauxx_u, tauyy_u, tauzz_u, tauxx_v, tauyy_v, tauzz_v 315 $ , last) 331 $ calving, tauxx_u, tauyy_u, tauzz_u, tauxx_v, tauyy_v, tauzz_v, 332 $ windsp, last) 333 c -- LOOP 316 334 c ====================================================================== 317 335 c S. Valcke (02/99) adapted From L.Z.X Li: this subroutine provides the … … 321 339 c ====================================================================== 322 340 IMPLICIT NONE 341 c -- LOOP 342 c 343 #include "dimensions.h" 344 INTEGER jjmp1 345 PARAMETER (jjmp1=jjm+1-1/jjm) 346 #include "dimphy.h" 347 c REAL zu10m(klon), zv10m(klon) 348 REAL zwindsp(klon) 349 c 350 c -- LOOP 351 c 352 323 353 INTEGER kt, imjm 324 354 c … … 345 375 REAL tauzz_u(imjm) 346 376 REAL tauzz_v(imjm) 377 c -- LOOP 378 REAL windsp(imjm) 379 c -- LOOP 347 380 LOGICAL last 348 381 c … … 424 457 c WRITE fields to files 425 458 DO jf=1, jpflda2o1 + jpflda2o2 426 IF (jf.eq. 1)459 IF (jf.eq.8) 427 460 $ CALL locwrite(cl_writ(jf),fsolice, imjm, 428 461 $ file_unit_field(jf), ierror) 429 IF (jf.eq. 2)462 IF (jf.eq.9) 430 463 $ CALL locwrite(cl_writ(jf),fsolwat, imjm, 431 464 $ file_unit_field(jf), ierror) 432 IF (jf.eq. 3)465 IF (jf.eq.10) 433 466 $ CALL locwrite(cl_writ(jf),fnsolice, imjm, 434 467 $ file_unit_field(jf), ierror) 435 IF (jf.eq. 4)468 IF (jf.eq.11) 436 469 $ CALL locwrite(cl_writ(jf),fnsolwat, imjm, 437 470 $ file_unit_field(jf), ierror) 438 IF (jf.eq. 5)471 IF (jf.eq.12) 439 472 $ CALL locwrite(cl_writ(jf),fnsicedt, imjm, 440 473 $ file_unit_field(jf), ierror) 441 c IF (jf.eq. 6)474 c IF (jf.eq.13) 442 475 c $ CALL locwrite(cl_writ(jf),ictemp, imjm, 443 476 c $ file_unit_field(jf), ierror) 444 IF (jf.eq. 6)477 IF (jf.eq.13) 445 478 $ CALL locwrite(cl_writ(jf),evice, imjm, 446 479 $ file_unit_field(jf), ierror) 447 IF (jf.eq. 7)480 IF (jf.eq.14) 448 481 $ CALL locwrite(cl_writ(jf),evwat, imjm, 449 482 $ file_unit_field(jf), ierror) 450 IF (jf.eq. 8)483 IF (jf.eq.15) 451 484 $ CALL locwrite(cl_writ(jf),lpre, imjm, 452 485 $ file_unit_field(jf), ierror) 453 IF (jf.eq. 9)486 IF (jf.eq.16) 454 487 $ CALL locwrite(cl_writ(jf),spre, imjm, 455 488 $ file_unit_field(jf), ierror) 456 IF (jf.eq.1 0)489 IF (jf.eq.17) 457 490 $ CALL locwrite(cl_writ(jf),dirunoff, imjm, 458 491 $ file_unit_field(jf), ierror) 459 IF (jf.eq.1 1)492 IF (jf.eq.18) 460 493 $ CALL locwrite(cl_writ(jf),rivrunoff, imjm, 461 494 $ file_unit_field(jf), ierror) 462 IF (jf.eq.1 2)495 IF (jf.eq.19) 463 496 $ CALL locwrite(cl_writ(jf),calving, imjm, 464 497 $ file_unit_field(jf), ierror) … … 475 508 c$$$ $ CALL locwrite(cl_writ(jf),tauyu, imjm, 476 509 c$$$ $ file_unit_field(jf), ierror) 477 IF (jf.eq.1 3)510 IF (jf.eq.1) 478 511 $ CALL locwrite(cl_writ(jf),tauxx_u, imjm, 479 512 $ file_unit_field(jf),ierror) 480 IF (jf.eq. 14)513 IF (jf.eq.2) 481 514 $ CALL locwrite(cl_writ(jf),tauyy_u, imjm, 482 515 $ file_unit_field(jf),ierror) 483 IF (jf.eq. 15)516 IF (jf.eq.3) 484 517 $ CALL locwrite(cl_writ(jf),tauzz_u, imjm, 485 518 $ file_unit_field(jf),ierror) 486 IF (jf.eq. 16)519 IF (jf.eq.4) 487 520 $ CALL locwrite(cl_writ(jf),tauxx_v, imjm, 488 521 $ file_unit_field(jf),ierror) 489 IF (jf.eq. 17)522 IF (jf.eq.5) 490 523 $ CALL locwrite(cl_writ(jf),tauyy_v, imjm, 491 524 $ file_unit_field(jf),ierror) 492 IF (jf.eq. 18)525 IF (jf.eq.6) 493 526 $ CALL locwrite(cl_writ(jf),tauzz_v, imjm, 494 527 $ file_unit_field(jf),ierror) 528 c -- LOOP 529 IF (jf.eq.7) 530 CALL locwrite(cl_writ(jf),windsp, imjm, 531 $ file_unit_field(jf),ierror) 532 c -- LOOP 533 495 534 END DO 496 535 C … … 522 561 DO jn=1, jpflda2o1 + jpflda2o2 523 562 C 524 IF (jn.eq. 1) CALL CLIM_Export(cl_writ(jn), kt, fsolice, info)525 IF (jn.eq. 2) CALL CLIM_Export(cl_writ(jn), kt, fsolwat, info)526 IF (jn.eq.3) CALL CLIM_Export(cl_writ(jn), kt, fnsolice, info)527 IF (jn.eq.4) CALL CLIM_Export(cl_writ(jn), kt, fnsolwat, info)528 IF (jn.eq.5) CALL CLIM_Export(cl_writ(jn), kt, fnsicedt, info)563 IF (jn.eq.8) CALL CLIM_Export(cl_writ(jn), kt, fsolice, info) 564 IF (jn.eq.9) CALL CLIM_Export(cl_writ(jn), kt, fsolwat, info) 565 IF (jn.eq.10) CALL CLIM_Export(cl_writ(jn), kt, fnsolice, info) 566 IF (jn.eq.11) CALL CLIM_Export(cl_writ(jn), kt, fnsolwat, info) 567 IF (jn.eq.12) CALL CLIM_Export(cl_writ(jn), kt, fnsicedt, info) 529 568 c IF (jn.eq.6) CALL CLIM_Export(cl_writ(jn), kt, ictemp, info) 530 IF (jn.eq. 6) CALL CLIM_Export(cl_writ(jn), kt, evice, info)531 IF (jn.eq. 7) CALL CLIM_Export(cl_writ(jn), kt, evwat, info)532 IF (jn.eq. 8) CALL CLIM_Export(cl_writ(jn), kt, lpre, info)533 IF (jn.eq. 9) CALL CLIM_Export(cl_writ(jn), kt, spre, info)534 IF (jn.eq.1 0) CALL CLIM_Export(cl_writ(jn),kt,dirunoff, info)535 IF (jn.eq.1 1) CALL CLIM_Export(cl_writ(jn),kt,rivrunoff,info)536 IF (jn.eq.1 2) CALL CLIM_Export(cl_writ(jn),kt,calving,info)569 IF (jn.eq.13) CALL CLIM_Export(cl_writ(jn), kt, evice, info) 570 IF (jn.eq.14) CALL CLIM_Export(cl_writ(jn), kt, evwat, info) 571 IF (jn.eq.15) CALL CLIM_Export(cl_writ(jn), kt, lpre, info) 572 IF (jn.eq.16) CALL CLIM_Export(cl_writ(jn), kt, spre, info) 573 IF (jn.eq.17) CALL CLIM_Export(cl_writ(jn),kt,dirunoff, info) 574 IF (jn.eq.18) CALL CLIM_Export(cl_writ(jn),kt,rivrunoff,info) 575 IF (jn.eq.19) CALL CLIM_Export(cl_writ(jn),kt,calving,info) 537 576 c$$$ IF (jn.eq.13) CALL CLIM_Export(cl_writ(jn), kt, tauxu, info) 538 577 c$$$ IF (jn.eq.14) CALL CLIM_Export(cl_writ(jn), kt, tauxv, info) 539 578 c$$$ IF (jn.eq.15) CALL CLIM_Export(cl_writ(jn), kt, tauyv, info) 540 579 c$$$ IF (jn.eq.16) CALL CLIM_Export(cl_writ(jn), kt, tauyu, info) 541 IF (jn.eq.13) CALL CLIM_Export(cl_writ(jn), kt, tauxx_u, info) 542 IF (jn.eq.14) CALL CLIM_Export(cl_writ(jn), kt, tauyy_u, info) 543 IF (jn.eq.15) CALL CLIM_Export(cl_writ(jn), kt, tauzz_u, info) 544 IF (jn.eq.16) CALL CLIM_Export(cl_writ(jn), kt, tauxx_v, info) 545 IF (jn.eq.17) CALL CLIM_Export(cl_writ(jn), kt, tauyy_v, info) 546 IF (jn.eq.18) CALL CLIM_Export(cl_writ(jn), kt, tauzz_v, info) 547 580 IF (jn.eq.1) CALL CLIM_Export(cl_writ(jn), kt, tauxx_u, info) 581 IF (jn.eq.2) CALL CLIM_Export(cl_writ(jn), kt, tauyy_u, info) 582 IF (jn.eq.3) CALL CLIM_Export(cl_writ(jn), kt, tauzz_u, info) 583 IF (jn.eq.4) CALL CLIM_Export(cl_writ(jn), kt, tauxx_v, info) 584 IF (jn.eq.5) CALL CLIM_Export(cl_writ(jn), kt, tauyy_v, info) 585 IF (jn.eq.6) CALL CLIM_Export(cl_writ(jn), kt, tauzz_v, info) 586 c -- LOOP 587 IF (jn.eq.7) CALL CLIM_Export(cl_writ(jn), kt, windsp, info) 588 c -- LOOP 548 589 IF (info .NE. CLIM_Ok) THEN 549 590 WRITE (nuout,*) 'STEP : Pb giving ',cl_writ(jn), ':',jn -
LMDZ4/trunk/libf/phylmd/oasis.true
r524 r590 2 2 ! $Header$ 3 3 ! 4 C $Id$5 4 C**** 6 5 C … … 20 19 c INCLUDE "param.h" 21 20 c 22 INTEGER kastp, kexch, kstep,imjm 21 INTEGER kastp, kexch, kstep,imjm,klon 23 22 INTEGER iparal(3) 24 23 INTEGER ifcpl, idt, info, imxtag, istep, jf 24 25 c -- LOOP 26 c 27 #include "dimensions.h" 28 INTEGER jjmp1 29 PARAMETER (jjmp1=jjm+1-1/jjm) 30 #include "dimphy.h" 31 c REAL zu10m(klon), zv10m(klon) 32 REAL zwindsp(klon) 33 c 34 c -- LOOP 25 35 c 26 36 #include "param_cou.h" … … 43 53 #include "mpiclim.h" 44 54 c 45 #include "oasis.h" ! contains the name of communication technique. Here 55 #include "oasis.h" 56 ! contains the name of communication technique. Here 46 57 ! cchan=CLIM only is possible. 47 58 c ! ctype=MPI2 … … 71 82 c must be the same as (1) of the field definition in namcouple: 72 83 c 73 cl_writ( 1)='COSHFICE'74 cl_writ( 2)='COSHFOCE'75 cl_writ( 3)='CONSFICE'76 cl_writ( 4)='CONSFOCE'77 cl_writ( 5)='CODFLXDT'84 cl_writ(8)='COSHFICE' 85 cl_writ(9)='COSHFOCE' 86 cl_writ(10)='CONSFICE' 87 cl_writ(11)='CONSFOCE' 88 cl_writ(12)='CODFLXDT' 78 89 c cl_writ(6)='COICTEMP' 79 cl_writ( 6)='COTFSICE'80 cl_writ( 7)='COTFSOCE'81 cl_writ( 8)='COTOLPSU'82 cl_writ( 9)='COTOSPSU'83 cl_writ(1 0)='CORUNCOA'84 cl_writ(1 1)='CORIVFLU'85 cl_writ(1 2)='COCALVIN'90 cl_writ(13)='COTFSICE' 91 cl_writ(14)='COTFSOCE' 92 cl_writ(15)='COTOLPSU' 93 cl_writ(16)='COTOSPSU' 94 cl_writ(17)='CORUNCOA' 95 cl_writ(18)='CORIVFLU' 96 cl_writ(19)='COCALVIN' 86 97 c$$$ cl_writ(13)='COZOTAUX' 87 98 c$$$ cl_writ(14)='COZOTAUV' 88 99 c$$$ cl_writ(15)='COMETAUY' 89 100 c$$$ cl_writ(16)='COMETAUU' 90 cl_writ(13)='COTAUXXU' 91 cl_writ(14)='COTAUYYU' 92 cl_writ(15)='COTAUZZU' 93 cl_writ(16)='COTAUXXV' 94 cl_writ(17)='COTAUYYV' 95 cl_writ(18)='COTAUZZV' 101 cl_writ(1)='COTAUXXU' 102 cl_writ(2)='COTAUYYU' 103 cl_writ(3)='COTAUZZU' 104 cl_writ(4)='COTAUXXV' 105 cl_writ(5)='COTAUYYV' 106 cl_writ(6)='COTAUZZV' 107 c -- LOOP 108 cl_writ(7)='COWINDSP' 109 c -- LOOP 96 110 c 97 111 c Define files name for fields exchanged from atmos to coupler, … … 116 130 cl_f_writ(17)='flxatmos' 117 131 cl_f_writ(18)='flxatmos' 132 c -- LOOP 133 cl_f_writ(19)='flxatmos' 134 c -- LOOP 118 135 119 136 c … … 311 328 SUBROUTINE intocpl(kt, imjm, fsolice, fsolwat, fnsolice, fnsolwat, 312 329 $ fnsicedt, evice, evwat, lpre, spre, dirunoff, rivrunoff, 313 $ calving, tauxx_u, tauyy_u, tauzz_u, tauxx_v, tauyy_v, tauzz_v 314 $ , last) 330 $ calving, tauxx_u, tauyy_u, tauzz_u, tauxx_v, tauyy_v, tauzz_v, 331 c -- LOOP 332 $ windsp, last) 333 c -- LOOP 315 334 c ====================================================================== 316 335 c S. Valcke (02/99) adapted From L.Z.X Li: this subroutine provides the … … 320 339 c ====================================================================== 321 340 IMPLICIT NONE 341 c -- LOOP 342 c 343 #include "dimensions.h" 344 INTEGER jjmp1 345 PARAMETER (jjmp1=jjm+1-1/jjm) 346 #include "dimphy.h" 347 c REAL zu10m(klon), zv10m(klon) 348 REAL zwindsp(klon) 349 c 350 c -- LOOP 351 c 352 322 353 INTEGER kt, imjm 323 354 c … … 344 375 REAL tauzz_u(imjm) 345 376 REAL tauzz_v(imjm) 377 c -- LOOP 378 REAL windsp(imjm) 379 c -- LOOP 346 380 LOGICAL last 347 381 c … … 423 457 c WRITE fields to files 424 458 DO jf=1, jpflda2o1 + jpflda2o2 425 IF (jf.eq. 1)459 IF (jf.eq.8) 426 460 $ CALL locwrite(cl_writ(jf),fsolice, imjm, 427 461 $ file_unit_field(jf), ierror) 428 IF (jf.eq. 2)462 IF (jf.eq.9) 429 463 $ CALL locwrite(cl_writ(jf),fsolwat, imjm, 430 464 $ file_unit_field(jf), ierror) 431 IF (jf.eq. 3)465 IF (jf.eq.10) 432 466 $ CALL locwrite(cl_writ(jf),fnsolice, imjm, 433 467 $ file_unit_field(jf), ierror) 434 IF (jf.eq. 4)468 IF (jf.eq.11) 435 469 $ CALL locwrite(cl_writ(jf),fnsolwat, imjm, 436 470 $ file_unit_field(jf), ierror) 437 IF (jf.eq. 5)471 IF (jf.eq.12) 438 472 $ CALL locwrite(cl_writ(jf),fnsicedt, imjm, 439 473 $ file_unit_field(jf), ierror) 440 c IF (jf.eq. 6)474 c IF (jf.eq.13) 441 475 c $ CALL locwrite(cl_writ(jf),ictemp, imjm, 442 476 c $ file_unit_field(jf), ierror) 443 IF (jf.eq. 6)477 IF (jf.eq.13) 444 478 $ CALL locwrite(cl_writ(jf),evice, imjm, 445 479 $ file_unit_field(jf), ierror) 446 IF (jf.eq. 7)480 IF (jf.eq.14) 447 481 $ CALL locwrite(cl_writ(jf),evwat, imjm, 448 482 $ file_unit_field(jf), ierror) 449 IF (jf.eq. 8)483 IF (jf.eq.15) 450 484 $ CALL locwrite(cl_writ(jf),lpre, imjm, 451 485 $ file_unit_field(jf), ierror) 452 IF (jf.eq. 9)486 IF (jf.eq.16) 453 487 $ CALL locwrite(cl_writ(jf),spre, imjm, 454 488 $ file_unit_field(jf), ierror) 455 IF (jf.eq.1 0)489 IF (jf.eq.17) 456 490 $ CALL locwrite(cl_writ(jf),dirunoff, imjm, 457 491 $ file_unit_field(jf), ierror) 458 IF (jf.eq.1 1)492 IF (jf.eq.18) 459 493 $ CALL locwrite(cl_writ(jf),rivrunoff, imjm, 460 494 $ file_unit_field(jf), ierror) 461 IF (jf.eq.1 2)495 IF (jf.eq.19) 462 496 $ CALL locwrite(cl_writ(jf),calving, imjm, 463 497 $ file_unit_field(jf), ierror) … … 474 508 c$$$ $ CALL locwrite(cl_writ(jf),tauyu, imjm, 475 509 c$$$ $ file_unit_field(jf), ierror) 476 IF (jf.eq.1 3)510 IF (jf.eq.1) 477 511 $ CALL locwrite(cl_writ(jf),tauxx_u, imjm, 478 512 $ file_unit_field(jf),ierror) 479 IF (jf.eq. 14)513 IF (jf.eq.2) 480 514 $ CALL locwrite(cl_writ(jf),tauyy_u, imjm, 481 515 $ file_unit_field(jf),ierror) 482 IF (jf.eq. 15)516 IF (jf.eq.3) 483 517 $ CALL locwrite(cl_writ(jf),tauzz_u, imjm, 484 518 $ file_unit_field(jf),ierror) 485 IF (jf.eq. 16)519 IF (jf.eq.4) 486 520 $ CALL locwrite(cl_writ(jf),tauxx_v, imjm, 487 521 $ file_unit_field(jf),ierror) 488 IF (jf.eq. 17)522 IF (jf.eq.5) 489 523 $ CALL locwrite(cl_writ(jf),tauyy_v, imjm, 490 524 $ file_unit_field(jf),ierror) 491 IF (jf.eq. 18)525 IF (jf.eq.6) 492 526 $ CALL locwrite(cl_writ(jf),tauzz_v, imjm, 493 527 $ file_unit_field(jf),ierror) 528 c -- LOOP 529 IF (jf.eq.7) 530 $ CALL locwrite(cl_writ(jf),windsp, imjm, 531 $ file_unit_field(jf),ierror) 532 c -- LOOP 533 494 534 END DO 495 535 C … … 521 561 DO jn=1, jpflda2o1 + jpflda2o2 522 562 C 523 IF (jn.eq. 1) CALL CLIM_Export(cl_writ(jn), kt, fsolice, info)524 IF (jn.eq. 2) CALL CLIM_Export(cl_writ(jn), kt, fsolwat, info)525 IF (jn.eq.3) CALL CLIM_Export(cl_writ(jn), kt, fnsolice, info)526 IF (jn.eq.4) CALL CLIM_Export(cl_writ(jn), kt, fnsolwat, info)527 IF (jn.eq.5) CALL CLIM_Export(cl_writ(jn), kt, fnsicedt, info)563 IF (jn.eq.8) CALL CLIM_Export(cl_writ(jn), kt, fsolice, info) 564 IF (jn.eq.9) CALL CLIM_Export(cl_writ(jn), kt, fsolwat, info) 565 IF (jn.eq.10) CALL CLIM_Export(cl_writ(jn), kt, fnsolice, info) 566 IF (jn.eq.11) CALL CLIM_Export(cl_writ(jn), kt, fnsolwat, info) 567 IF (jn.eq.12) CALL CLIM_Export(cl_writ(jn), kt, fnsicedt, info) 528 568 c IF (jn.eq.6) CALL CLIM_Export(cl_writ(jn), kt, ictemp, info) 529 IF (jn.eq. 6) CALL CLIM_Export(cl_writ(jn), kt, evice, info)530 IF (jn.eq. 7) CALL CLIM_Export(cl_writ(jn), kt, evwat, info)531 IF (jn.eq. 8) CALL CLIM_Export(cl_writ(jn), kt, lpre, info)532 IF (jn.eq. 9) CALL CLIM_Export(cl_writ(jn), kt, spre, info)533 IF (jn.eq.1 0) CALL CLIM_Export(cl_writ(jn),kt,dirunoff, info)534 IF (jn.eq.1 1) CALL CLIM_Export(cl_writ(jn),kt,rivrunoff,info)535 IF (jn.eq.1 2) CALL CLIM_Export(cl_writ(jn),kt,calving,info)569 IF (jn.eq.13) CALL CLIM_Export(cl_writ(jn), kt, evice, info) 570 IF (jn.eq.14) CALL CLIM_Export(cl_writ(jn), kt, evwat, info) 571 IF (jn.eq.15) CALL CLIM_Export(cl_writ(jn), kt, lpre, info) 572 IF (jn.eq.16) CALL CLIM_Export(cl_writ(jn), kt, spre, info) 573 IF (jn.eq.17) CALL CLIM_Export(cl_writ(jn),kt,dirunoff, info) 574 IF (jn.eq.18) CALL CLIM_Export(cl_writ(jn),kt,rivrunoff,info) 575 IF (jn.eq.19) CALL CLIM_Export(cl_writ(jn),kt,calving,info) 536 576 c$$$ IF (jn.eq.13) CALL CLIM_Export(cl_writ(jn), kt, tauxu, info) 537 577 c$$$ IF (jn.eq.14) CALL CLIM_Export(cl_writ(jn), kt, tauxv, info) 538 578 c$$$ IF (jn.eq.15) CALL CLIM_Export(cl_writ(jn), kt, tauyv, info) 539 579 c$$$ IF (jn.eq.16) CALL CLIM_Export(cl_writ(jn), kt, tauyu, info) 540 IF (jn.eq.13) CALL CLIM_Export(cl_writ(jn), kt, tauxx_u, info) 541 IF (jn.eq.14) CALL CLIM_Export(cl_writ(jn), kt, tauyy_u, info) 542 IF (jn.eq.15) CALL CLIM_Export(cl_writ(jn), kt, tauzz_u, info) 543 IF (jn.eq.16) CALL CLIM_Export(cl_writ(jn), kt, tauxx_v, info) 544 IF (jn.eq.17) CALL CLIM_Export(cl_writ(jn), kt, tauyy_v, info) 545 IF (jn.eq.18) CALL CLIM_Export(cl_writ(jn), kt, tauzz_v, info) 546 580 IF (jn.eq.1) CALL CLIM_Export(cl_writ(jn), kt, tauxx_u, info) 581 IF (jn.eq.2) CALL CLIM_Export(cl_writ(jn), kt, tauyy_u, info) 582 IF (jn.eq.3) CALL CLIM_Export(cl_writ(jn), kt, tauzz_u, info) 583 IF (jn.eq.4) CALL CLIM_Export(cl_writ(jn), kt, tauxx_v, info) 584 IF (jn.eq.5) CALL CLIM_Export(cl_writ(jn), kt, tauyy_v, info) 585 IF (jn.eq.6) CALL CLIM_Export(cl_writ(jn), kt, tauzz_v, info) 586 c -- LOOP 587 IF (jn.eq.7) CALL CLIM_Export(cl_writ(jn), kt, windsp, info) 588 c -- LOOP 547 589 IF (info .NE. CLIM_Ok) THEN 548 590 WRITE (nuout,*) 'STEP : Pb giving ',cl_writ(jn), ':',jn
Note: See TracChangeset
for help on using the changeset viewer.
