14 |
! ----------------------------------------------------------------- |
! ----------------------------------------------------------------- |
15 |
use dimens_m |
use dimens_m |
16 |
use paramet_m |
use paramet_m |
17 |
use comvert |
use disvert_m |
18 |
use comgeom |
use comgeom |
19 |
use conf_guide_m |
use conf_guide_m |
20 |
use netcdf |
use netcdf |
24 |
! common |
! common |
25 |
! ------ |
! ------ |
26 |
|
|
|
include "netcdf.inc" |
|
|
|
|
|
|
|
27 |
! arguments |
! arguments |
28 |
! --------- |
! --------- |
29 |
integer nlevnc |
integer nlevnc |
120 |
! Niveaux de pression |
! Niveaux de pression |
121 |
print*,'WARNING!!! Il n y a pas de test de coherence' |
print*,'WARNING!!! Il n y a pas de test de coherence' |
122 |
print*,'sur le nombre de niveaux verticaux dans le fichier nc' |
print*,'sur le nombre de niveaux verticaux dans le fichier nc' |
123 |
status=NF_GET_VARA_REAL(ncidpl,varidpl,1,nlevnc,pl) |
status=NF90_GET_VAR(ncidpl,varidpl,pl) |
124 |
! passage en pascal |
! passage en pascal |
125 |
pl(:)=100.*pl(:) |
pl(:)=100.*pl(:) |
126 |
if (first) then |
if (first) then |
158 |
|
|
159 |
if (guide_u) then |
if (guide_u) then |
160 |
print*,'avant la lecture de UNCEP nd de niv:',nlevnc |
print*,'avant la lecture de UNCEP nd de niv:',nlevnc |
161 |
status=NF_GET_VARA_REAL(ncidu,varidu,start,count,unc) |
status=NF90_GET_VAR(ncidu,varidu,unc,start,count) |
162 |
! call dump2d(iip1,jjp1,unc,'VENT NCEP ') |
! call dump2d(iip1,jjp1,unc,'VENT NCEP ') |
163 |
! call dump2d(iip1,40,unc(1,1,nlevnc),'VENT NCEP ') |
! call dump2d(iip1,40,unc(1,1,nlevnc),'VENT NCEP ') |
164 |
print*,'WARNING!!! Correction bidon pour palier a un ' |
print*,'WARNING!!! Correction bidon pour palier a un ' |
173 |
print*,'ncidt=',ncidt,'varidt=',varidt,'start=',start |
print*,'ncidt=',ncidt,'varidt=',varidt,'start=',start |
174 |
print*,'count=',count |
print*,'count=',count |
175 |
if (guide_T) then |
if (guide_T) then |
176 |
status=NF_GET_VARA_REAL(ncidt,varidt,start,count,tnc) |
status=NF90_GET_VAR(ncidt,varidt,tnc,start,count) |
177 |
call dump2d(iip1,jjp1,tnc,'TNC COUCHE 1 AAA ') |
call dump2d(iip1,jjp1,tnc,'TNC COUCHE 1 AAA ') |
178 |
call correctbid(iim,jjp1*nlevnc,tnc) |
call correctbid(iim,jjp1*nlevnc,tnc) |
179 |
call dump2d(iip1,jjp1,tnc,'TNC COUCHE 1 BBB ') |
call dump2d(iip1,jjp1,tnc,'TNC COUCHE 1 BBB ') |
183 |
! -------- |
! -------- |
184 |
|
|
185 |
if (guide_Q) then |
if (guide_Q) then |
186 |
status=NF_GET_VARA_REAL(ncidQ,varidQ,start,count,Qnc) |
status=NF90_GET_VAR(ncidQ,varidQ,Qnc,start,count) |
187 |
call correctbid(iim,jjp1*nlevnc,Qnc) |
call correctbid(iim,jjp1*nlevnc,Qnc) |
188 |
call dump2d(iip1,jjp1,Qnc,'QNC COUCHE 1 ') |
call dump2d(iip1,jjp1,Qnc,'QNC COUCHE 1 ') |
189 |
endif |
endif |
193 |
! ------------- |
! ------------- |
194 |
|
|
195 |
if (guide_v) then |
if (guide_v) then |
196 |
status=NF_GET_VARA_REAL(ncidv,varidv,start,count,vnc) |
status=NF90_GET_VAR(ncidv,varidv,vnc,start,count) |
197 |
call correctbid(iim,jjm*nlevnc,vnc) |
call correctbid(iim,jjm*nlevnc,vnc) |
198 |
call dump2d(iip1,jjm,vnc,'VNC COUCHE 1 ') |
call dump2d(iip1,jjm,vnc,'VNC COUCHE 1 ') |
199 |
endif |
endif |
208 |
! ------------------- |
! ------------------- |
209 |
|
|
210 |
if (guide_P) then |
if (guide_P) then |
211 |
status=NF_GET_VARA_REAL(ncidps,varidps,start,count,psnc) |
status=NF90_GET_VAR(ncidps,varidps,psnc,start,count) |
212 |
call dump2d(iip1,jjp1,psnc,'PSNC COUCHE 1 ') |
call dump2d(iip1,jjp1,psnc,'PSNC COUCHE 1 ') |
213 |
call correctbid(iim,jjp1,psnc) |
call correctbid(iim,jjp1,psnc) |
214 |
endif |
endif |