Changeset 590 for LMDZ4/trunk/libf/phylmd/oasis.dummy
- Timestamp:
- Feb 7, 2005, 5:41:35 PM (19 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
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
Note: See TracChangeset
for help on using the changeset viewer.