Changeset 719 for LMDZ4/trunk/libf/phylmd/physiq.F
- Timestamp:
- Oct 6, 2006, 12:34:58 PM (18 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
LMDZ4/trunk/libf/phylmd/physiq.F
r702 r719 12 12 #endif 13 13 . d_u, d_v, d_t, d_qx, d_ps 14 cIM Amip215 14 . , dudyn 16 15 . , PVteta) … … 38 37 c CLEFS CPP POUR LES IO 39 38 c ===================== 39 #define histins 40 40 #define histhf 41 41 #define histday 42 42 #define histmth 43 #define histins44 43 #define histmthNMC 45 44 #define histISCCP … … 69 68 c d_q_dyn-input-R-tendance dynamique pour "q" (kg/kg/s) 70 69 c omega---input-R-vitesse verticale en Pa/s 71 cIM comgeomphy.h BEG72 70 c cuphy----input-R-resolution des mailles en x (m) 73 71 c cvphy----input-R-resolution des mailles en y (m) 74 cIM comgeomphy.h END75 72 c d_u-----output-R-tendance physique de "u" (m/s/s) 76 73 c d_v-----output-R-tendance physique de "v" (m/s/s) … … 78 75 c d_qx----output-R-tendance physique de "qx" (kg/kg/s) 79 76 c d_ps----output-R-tendance physique de la pression au sol 80 cIM81 77 c PVteta--output-R-vorticite potentielle a des thetas constantes 82 78 c====================================================================== … … 127 123 cIM "slab" ocean 128 124 REAL tslab(klon) !Temperature du slab-ocean 129 cIM SAVE tslab130 125 REAL seaice(klon) !glace de mer (kg/m2) 131 cIM SAVE seaice132 126 REAL fluxo(klon) !flux turbulents ocean-glace de mer 133 127 REAL fluxg(klon) !flux turbulents ocean-atmosphere 134 cIM135 128 REAL amn, amx 136 129 c====================================================================== … … 206 199 REAL znivsig(klev) 207 200 REAL zsurf(nbsrf) 208 cIM209 201 INTEGER kinv 210 202 real pir 211 cMI 203 212 204 REAL u(klon,klev) 213 205 REAL v(klon,klev) … … 464 456 INTEGER iw, iwmax 465 457 REAL wmin, pas_w 466 c PARAMETER(wmin=-100.,pas_w=10.,iwmax=30)467 cIM 051005 PARAMETER(wmin=-200.,pas_w=10.,iwmax=40)468 458 PARAMETER(wmin=-100.,pas_w=10.,iwmax=20) 469 459 REAL o500(klon) … … 473 463 INTEGER nreg, nbregdyn 474 464 PARAMETER(nbregdyn=5) 475 cIM 051005 BEG476 c REAL histoW(iwmax,nbregdyn,napisccp)477 c REAL nhistoW(iwmax,nbregdyn,napisccp)478 c REAL histoWn(iwmax,nbregdyn)479 c REAL nhistoWn(iwmax,nbregdyn)480 cIM 090905 END481 465 482 466 INTEGER linv … … 486 470 c sorties ISCCP 487 471 488 c logical ok_isccp489 c real ecrit_isccp490 472 integer nid_isccp 491 c save ok_isccp, ecrit_isccp, nid_isccp492 473 save nid_isccp 493 494 #undef histISCCP495 #define histISCCP496 c data ok_isccp,ecrit_isccp/.true.,0.125/497 c data ok_isccp,ecrit_isccp/.true.,1./498 cIM 190504 data ok_isccp/.true./499 cIM 190504 #else500 cIM 190504 data ok_isccp/.false./501 cIM 190504 #endif502 474 503 475 REAL zx_tau(kmaxm1), zx_pc(lmaxm1), zx_o500(iwmax) 504 476 DATA zx_tau/0.0, 0.3, 1.3, 3.6, 9.4, 23., 60./ 505 cIM bad 151205 DATA zx_pc/50., 180., 310., 440., 560., 680., 800./506 477 DATA zx_pc/180., 310., 440., 560., 680., 800., 1000./ 507 478 508 479 c cldtopres pression au sommet des nuages 509 480 REAL cldtopres(lmaxm1), cldtopres3(lmax3) 510 cIM 151205 erreur DATA cldtopres/50., 180., 310., 440., 560., 680., 800./511 481 DATA cldtopres/180., 310., 440., 560., 680., 800., 1000./ 512 482 DATA cldtopres3/440., 680., 1000./ 513 cIM 051005 BEG514 483 REAL tmp_his1_3d(iwmax,kmaxm1,lmaxm1,nbregdyn,napisccp) 515 484 REAL tmp_his2_3d(iwmax,kmaxm1,lmaxm1,nbregdyn,napisccp) 516 485 REAL tmp_his3_3d(iwmax,kmaxm1,lmaxm1,nbregdyn,napisccp) 517 cIM 051005 END518 519 486 INTEGER komega, nhoriRD 520 487 … … 594 561 c QUESTION : noms de variables ? 595 562 596 #undef histhf597 #define histhf598 563 #ifdef histhf 599 cIM 130904 data ok_hf,ecrit_hf/.true.,0.25/600 564 data ok_hf/.true./ 601 565 #else … … 1615 1579 #endif 1616 1580 1617 #undef histmthNMC1618 #define histmthNMC1619 1581 #ifdef histmthNMC 1620 1582 #include "ini_histmthNMC.h"
Note: See TracChangeset
for help on using the changeset viewer.