4 |
!=========================================================================== |
!=========================================================================== |
5 |
|
|
6 |
use dimens_m |
use dimens_m |
|
use dump2d_m, only: dump2d |
|
7 |
use paramet_m |
use paramet_m |
8 |
use comconst |
use comconst |
9 |
use disvert_m |
use disvert_m |
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 |
|
|