Changeset 403
- Timestamp:
- Oct 15, 2002, 4:28:41 PM (22 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
LMDZ.3.3/branches/rel-LF/libf/dyn3d/gcm.F
r353 r403 126 126 LOGICAL first 127 127 REAL alpha(ip1jmp1,llm),beta(ip1jmp1,llm) 128 c+jld variables test conservation energie 129 REAL ecin(ip1jmp1,llm),ecin0(ip1jmp1,llm) 130 C Tendance de la temp. potentiel d (theta)/ d t due a la 131 C tansformation d'energie cinetique en energie thermique 132 C cree par la dissipation 133 REAL dhecdt(ip1jmp1,llm) 134 REAL vcont(ip1jm,llm),ucont(ip1jmp1,llm) 135 REAL d_h_vcol, d_qt, d_qw, d_ql, d_ec 136 CHARACTER*15 ztit 137 INTEGER ip_ebil ! PRINT level for energy conserv. diag. 138 SAVE ip_ebil 139 DATA ip_ebil/1/ 140 c-jld 128 141 129 142 LOGICAL offline ! Controle du stockage ds "fluxmass" … … 274 287 CALL inidissip( lstardis, nitergdiv, nitergrot, niterh , 275 288 * tetagdiv, tetagrot , tetatemp ) 289 c+jld 290 C initialisation constantes thermo utilisees dans diagedyn 291 IF (ip_ebil.ge.1 ) THEN 292 CALL suphec 293 ENDIF 294 c-jld 276 295 c 277 296 … … 327 346 328 347 1 CONTINUE 348 349 329 350 330 351 if (ok_nudge) then … … 475 496 c rajout debug 476 497 c lafin = .true. 477 498 c+jld 499 IF (ip_ebil.ge.1 ) THEN 500 ztit='bil dyn' 501 CALL diagedyn(ztit,2,1,1,dtphys 502 e , ucov , vcov , ps, p ,pk , teta , q(:,:,1), q(:,:,2),aire 503 s , d_h_vcol , d_qt, d_qw, d_ql, d_ec) 504 ENDIF 505 c-jld 478 506 CALL calfis( nqmx, lafin ,rdayvrai,rday_ecri,time , 479 507 $ ucov,vcov,teta,q,masse,ps,p,pk,phis,phi , … … 486 514 $ dufi, dvfi, dhfi , dqfi ,dpfi ) 487 515 c 516 c+jld 517 IF (ip_ebil.ge.1 ) THEN 518 ztit='bil phys' 519 CALL diagedyn(ztit,2,1,1,dtphys 520 e , ucov , vcov , ps, p ,pk , teta , q(:,:,1), q(:,:,2),aire 521 s , d_h_vcol , d_qt, d_qw, d_ql, d_ec) 522 ENDIF 523 c-jld 488 524 ENDIF 489 525 … … 498 534 499 535 IF(apdiss) THEN 536 c+jld 537 IF (ip_ebil.ge.2 ) THEN 538 ztit='avant dissip' 539 CALL diagedyn(ztit,2,2,2,dtvr 540 e , ucov , vcov , ps, p ,pk , teta , q(:,:,1), q(:,:,2),aire 541 s , d_h_vcol , d_qt, d_qw, d_ql, d_ec) 542 ENDIF 543 CALL covcont ( llm , ucov , vcov , ucont, vcont ) 544 CALL enercin ( vcov , ucov , vcont , ucont , ecin0 ) 545 c-jld 500 546 501 547 CALL dissip(vcov,ucov,teta,p,dvdis,dudis,dhdis) … … 503 549 CALL addit( ijp1llm,ucov ,dudis,ucov ) 504 550 CALL addit( ijmllm ,vcov ,dvdis,vcov ) 551 c+jld 552 CALL covcont ( llm , ucov , vcov , ucont, vcont ) 553 CALL enercin ( vcov , ucov , vcont , ucont , ecin ) 554 C 555 C On rajoute la tendance due a la transform. Ec -> E therm. cree 556 C lors de la dissipation 557 DO l = 1, llm 558 DO ij = 1, ip1jmp1 559 dhecdt(ij,l)= (ecin0(ij,l)-ecin(ij,l))/ pk(ij,l) 560 dhdis(ij,l) = dhdis(ij,l) + dhecdt(ij,l) 561 ENDDO 562 ENDDO 563 c-jld 505 564 CALL addit( ijp1llm,teta ,dhdis,teta ) 506 565 c+jld 566 c 567 IF (ip_ebil.ge.2 ) THEN 568 ztit='apres dissip' 569 CALL diagedyn(ztit,2,2,2,dtdiss 570 e , ucov , vcov , ps, p ,pk , teta , q(:,:,1), q(:,:,2),aire 571 s , d_h_vcol , d_qt, d_qw, d_ql, d_ec) 572 ENDIF 573 c-jld 507 574 508 575 c ....... P. Le Van ( ajout le 17/04/96 ) ........... … … 570 637 571 638 IF( itau. EQ. itaufinp1 ) then 639 572 640 abort_message = 'Simulation finished' 573 641 call abort_gcm(modname,abort_message,0) … … 585 653 CALL writedynav(histaveid, nqmx, itau,vcov , 586 654 , ucov,teta,pk,phi,q,masse,ps,phis) 655 cccIM cf. FH 656 call bilan_dyn (2,dtvr*iperiod,dtvr*day_step*periodav, 657 , ps,masse,pk,pbaru,pbarv,teta,phi,ucov,vcov,q) 658 587 659 ENDIF 588 660 … … 676 748 CALL writedynav(histaveid, nqmx, itau,vcov , 677 749 , ucov,teta,pk,phi,q,masse,ps,phis) 750 cccIM cf. FH 751 call bilan_dyn (2,dtvr*iperiod,dtvr*day_step*periodav, 752 , ps,masse,pk,pbaru,pbarv,teta,phi,ucov,vcov,q) 678 753 679 754 ENDIF
Note: See TracChangeset
for help on using the changeset viewer.