6 |
use dimens_m |
use dimens_m |
7 |
use paramet_m |
use paramet_m |
8 |
use comconst |
use comconst |
9 |
use comvert |
use disvert_m |
10 |
use comgeom |
use comgeom |
11 |
use q_sat_m, only: q_sat |
use q_sat_m, only: q_sat |
12 |
use guide_m |
use guide_m |
15 |
|
|
16 |
real u(iip1,jjp1,llm),v(iip1,jjm,llm) |
real u(iip1,jjp1,llm),v(iip1,jjm,llm) |
17 |
real t(iip1,jjp1,llm),pk(iip1,jjp1,llm),rh(iip1,jjp1,llm) |
real t(iip1,jjp1,llm),pk(iip1,jjp1,llm),rh(iip1,jjp1,llm) |
|
real ps(iip1,jjp1) |
|
|
|
|
18 |
real ucov(iip1,jjp1,llm),vcov(iip1,jjm,llm) |
real ucov(iip1,jjp1,llm),vcov(iip1,jjm,llm) |
19 |
real teta(iip1,jjp1,llm),q(iip1,jjp1,llm) |
real teta(iip1,jjp1,llm),q(iip1,jjp1,llm) |
20 |
|
|
63 |
enddo |
enddo |
64 |
enddo |
enddo |
65 |
|
|
|
! call dump2d(iip1,jjp1,teta,'TETA EN BAS ') |
|
|
! call dump2d(iip1,jjp1,teta(1,1,llm),'TETA EN HAUT ') |
|
|
|
|
66 |
! Humidite relative -> specifique |
! Humidite relative -> specifique |
67 |
! ------------------------------- |
! ------------------------------- |
68 |
if (1.eq.0) then |
if (1.eq.0) then |
79 |
q(:,:,:)=qsat(:,:,:)*rh(:,:,:) |
q(:,:,:)=qsat(:,:,:)*rh(:,:,:) |
80 |
print*,'calcul de q OK' |
print*,'calcul de q OK' |
81 |
|
|
|
call dump2d(iip1,jjp1,pres,'PRESSION PREMIERE COUCHE ') |
|
|
call dump2d(iip1,jjp1,q,'HUMIDITE SPECIFIQUE COUCHE 1 ') |
|
82 |
endif |
endif |
83 |
|
|
84 |
|
|