- 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.dummy
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 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 … … 312 328 SUBROUTINE intocpl(kt, imjm, fsolice, fsolwat, fnsolice, fnsolwat, 313 329 $ fnsicedt, evice, evwat, lpre, spre, dirunoff, rivrunoff, 314 $ calving, tauxx_u, tauyy_u, tauzz_u, tauxx_v, tauyy_v, tauzz_v 315 $ , last) 330 $ calving, tauxx_u, tauyy_u, tauzz_u, tauxx_v, tauyy_v, tauzz_v, 331 $ windsp, last) 332 c -- LOOP 316 333 c ====================================================================== 317 334 c S. Valcke (02/99) adapted From L.Z.X Li: this subroutine provides the … … 321 338 c ====================================================================== 322 339 IMPLICIT NONE 340 c -- LOOP 341 c 342 #include "dimensions.h" 343 INTEGER jjmp1 344 PARAMETER (jjmp1=jjm+1-1/jjm) 345 #include "dimphy.h" 346 c REAL zu10m(klon), zv10m(klon) 347 REAL zwindsp(klon) 348 c 349 c -- LOOP 350 c 351 323 352 INTEGER kt, imjm 324 353 c … … 345 374 REAL tauzz_u(imjm) 346 375 REAL tauzz_v(imjm) 376 c -- LOOP 377 REAL windsp(imjm) 378 c -- LOOP 347 379 LOGICAL last 348 380 c … … 424 456 c WRITE fields to files 425 457 DO jf=1, jpflda2o1 + jpflda2o2 426 IF (jf.eq. 1)458 IF (jf.eq.8) 427 459 $ CALL locwrite(cl_writ(jf),fsolice, imjm, 428 460 $ file_unit_field(jf), ierror) 429 IF (jf.eq. 2)461 IF (jf.eq.9) 430 462 $ CALL locwrite(cl_writ(jf),fsolwat, imjm, 431 463 $ file_unit_field(jf), ierror) 432 IF (jf.eq. 3)464 IF (jf.eq.10) 433 465 $ CALL locwrite(cl_writ(jf),fnsolice, imjm, 434 466 $ file_unit_field(jf), ierror) 435 IF (jf.eq. 4)467 IF (jf.eq.11) 436 468 $ CALL locwrite(cl_writ(jf),fnsolwat, imjm, 437 469 $ file_unit_field(jf), ierror) 438 IF (jf.eq. 5)470 IF (jf.eq.12) 439 471 $ CALL locwrite(cl_writ(jf),fnsicedt, imjm, 440 472 $ file_unit_field(jf), ierror) 441 c IF (jf.eq. 6)473 c IF (jf.eq.13) 442 474 c $ CALL locwrite(cl_writ(jf),ictemp, imjm, 443 475 c $ file_unit_field(jf), ierror) 444 IF (jf.eq. 6)476 IF (jf.eq.13) 445 477 $ CALL locwrite(cl_writ(jf),evice, imjm, 446 478 $ file_unit_field(jf), ierror) 447 IF (jf.eq. 7)479 IF (jf.eq.14) 448 480 $ CALL locwrite(cl_writ(jf),evwat, imjm, 449 481 $ file_unit_field(jf), ierror) 450 IF (jf.eq. 8)482 IF (jf.eq.15) 451 483 $ CALL locwrite(cl_writ(jf),lpre, imjm, 452 484 $ file_unit_field(jf), ierror) 453 IF (jf.eq. 9)485 IF (jf.eq.16) 454 486 $ CALL locwrite(cl_writ(jf),spre, imjm, 455 487 $ file_unit_field(jf), ierror) 456 IF (jf.eq.1 0)488 IF (jf.eq.17) 457 489 $ CALL locwrite(cl_writ(jf),dirunoff, imjm, 458 490 $ file_unit_field(jf), ierror) 459 IF (jf.eq.1 1)491 IF (jf.eq.18) 460 492 $ CALL locwrite(cl_writ(jf),rivrunoff, imjm, 461 493 $ file_unit_field(jf), ierror) 462 IF (jf.eq.1 2)494 IF (jf.eq.19) 463 495 $ CALL locwrite(cl_writ(jf),calving, imjm, 464 496 $ file_unit_field(jf), ierror) … … 475 507 c$$$ $ CALL locwrite(cl_writ(jf),tauyu, imjm, 476 508 c$$$ $ file_unit_field(jf), ierror) 477 IF (jf.eq.1 3)509 IF (jf.eq.1) 478 510 $ CALL locwrite(cl_writ(jf),tauxx_u, imjm, 479 511 $ file_unit_field(jf),ierror) 480 IF (jf.eq. 14)512 IF (jf.eq.2) 481 513 $ CALL locwrite(cl_writ(jf),tauyy_u, imjm, 482 514 $ file_unit_field(jf),ierror) 483 IF (jf.eq. 15)515 IF (jf.eq.3) 484 516 $ CALL locwrite(cl_writ(jf),tauzz_u, imjm, 485 517 $ file_unit_field(jf),ierror) 486 IF (jf.eq. 16)518 IF (jf.eq.4) 487 519 $ CALL locwrite(cl_writ(jf),tauxx_v, imjm, 488 520 $ file_unit_field(jf),ierror) 489 IF (jf.eq. 17)521 IF (jf.eq.5) 490 522 $ CALL locwrite(cl_writ(jf),tauyy_v, imjm, 491 523 $ file_unit_field(jf),ierror) 492 IF (jf.eq. 18)524 IF (jf.eq.6) 493 525 $ CALL locwrite(cl_writ(jf),tauzz_v, imjm, 494 526 $ file_unit_field(jf),ierror) 527 c -- LOOP 528 IF (jf.eq.7) 529 CALL locwrite(cl_writ(jf),windsp, imjm, 530 $ file_unit_field(jf),ierror) 531 c -- LOOP 532 495 533 END DO 496 534 C … … 522 560 DO jn=1, jpflda2o1 + jpflda2o2 523 561 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)562 IF (jn.eq.8) CALL CLIM_Export(cl_writ(jn), kt, fsolice, info) 563 IF (jn.eq.9) CALL CLIM_Export(cl_writ(jn), kt, fsolwat, info) 564 IF (jn.eq.10) CALL CLIM_Export(cl_writ(jn), kt, fnsolice, info) 565 IF (jn.eq.11) CALL CLIM_Export(cl_writ(jn), kt, fnsolwat, info) 566 IF (jn.eq.12) CALL CLIM_Export(cl_writ(jn), kt, fnsicedt, info) 529 567 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)568 IF (jn.eq.13) CALL CLIM_Export(cl_writ(jn), kt, evice, info) 569 IF (jn.eq.14) CALL CLIM_Export(cl_writ(jn), kt, evwat, info) 570 IF (jn.eq.15) CALL CLIM_Export(cl_writ(jn), kt, lpre, info) 571 IF (jn.eq.16) CALL CLIM_Export(cl_writ(jn), kt, spre, info) 572 IF (jn.eq.17) CALL CLIM_Export(cl_writ(jn),kt,dirunoff, info) 573 IF (jn.eq.18) CALL CLIM_Export(cl_writ(jn),kt,rivrunoff,info) 574 IF (jn.eq.19) CALL CLIM_Export(cl_writ(jn),kt,calving,info) 537 575 c$$$ IF (jn.eq.13) CALL CLIM_Export(cl_writ(jn), kt, tauxu, info) 538 576 c$$$ IF (jn.eq.14) CALL CLIM_Export(cl_writ(jn), kt, tauxv, info) 539 577 c$$$ IF (jn.eq.15) CALL CLIM_Export(cl_writ(jn), kt, tauyv, info) 540 578 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 579 IF (jn.eq.1) CALL CLIM_Export(cl_writ(jn), kt, tauxx_u, info) 580 IF (jn.eq.2) CALL CLIM_Export(cl_writ(jn), kt, tauyy_u, info) 581 IF (jn.eq.3) CALL CLIM_Export(cl_writ(jn), kt, tauzz_u, info) 582 IF (jn.eq.4) CALL CLIM_Export(cl_writ(jn), kt, tauxx_v, info) 583 IF (jn.eq.5) CALL CLIM_Export(cl_writ(jn), kt, tauyy_v, info) 584 IF (jn.eq.6) CALL CLIM_Export(cl_writ(jn), kt, tauzz_v, info) 585 c -- LOOP 586 IF (jn.eq.7) CALL CLIM_Export(cl_writ(jn), kt, windsp, info) 587 c -- LOOP 548 588 IF (info .NE. CLIM_Ok) THEN 549 589 WRITE (nuout,*) 'STEP : Pb giving ',cl_writ(jn), ':',jn
Note: See TracChangeset
for help on using the changeset viewer.
