- Timestamp:
- Feb 7, 2005, 4:47:11 PM (19 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
LMDZ4/branches/IPSL-CM4_IPCC_patches/libf/phylmd/oasis.true
r524 r589 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 c REAL zu10m(klon), zv10m(klon) 33 REAL zwindsp(klon) 34 c 35 c -- LOOP 25 36 c 26 37 #include "param_cou.h" … … 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 … … 309 326 310 327 c $Id$ 328 c -- LOOP 311 329 SUBROUTINE intocpl(kt, imjm, fsolice, fsolwat, fnsolice, fnsolwat, 312 330 $ fnsicedt, evice, evwat, lpre, spre, dirunoff, rivrunoff, 313 $ calving, tauxx_u, tauyy_u, tauzz_u, tauxx_v, tauyy_v, tauzz_v 314 $ , last) 331 $ calving, tauxx_u, tauyy_u, tauzz_u, tauxx_v, tauyy_v, tauzz_v, 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.
