1 | subroutine read_surface(name,surfa) |
---|
2 | |
---|
3 | |
---|
4 | ! common |
---|
5 | ! ------ |
---|
6 | USE ioipsl |
---|
7 | USE comgeomphy |
---|
8 | USE dimphy |
---|
9 | USE mod_grid_phy_lmdz |
---|
10 | USE mod_phys_lmdz_para |
---|
11 | USE iophy |
---|
12 | ! USE netcdf |
---|
13 | IMPLICIT NONE |
---|
14 | |
---|
15 | include "netcdf.inc" |
---|
16 | #include "dimensions.h" |
---|
17 | #include "paramet.h" |
---|
18 | |
---|
19 | character*10 name |
---|
20 | character*10 varname |
---|
21 | ! |
---|
22 | real tmp_dyn(iip1,jjp1) |
---|
23 | REAL tmp_dyn_invers(iip1,jjp1) |
---|
24 | real tmp_fi(klon) |
---|
25 | real surfa(klon,5) |
---|
26 | ! |
---|
27 | integer ncid |
---|
28 | integer varid |
---|
29 | real rcode |
---|
30 | integer start(2),count(2),status |
---|
31 | integer i,j,l,ig |
---|
32 | character*1 str1 |
---|
33 | |
---|
34 | print*,'Lecture du fichier donnees_lisa.nc' |
---|
35 | ncid=NCOPN('donnees_lisa.nc',NCNOWRIT,rcode) |
---|
36 | |
---|
37 | DO i=1,5 |
---|
38 | write(str1,'(i1)') i |
---|
39 | varname=trim(name)//str1 |
---|
40 | print*,'lecture variable:',varname |
---|
41 | varid=NCVID(ncid,trim(varname),rcode) |
---|
42 | ! varid=NCVID(ncid,varname,rcode) |
---|
43 | |
---|
44 | ! dimensions pour les champs scalaires et le vent zonal |
---|
45 | ! ----------------------------------------------------- |
---|
46 | |
---|
47 | start(1)=1 |
---|
48 | start(2)=1 |
---|
49 | count(1)=iip1 |
---|
50 | count(2)=jjp1 |
---|
51 | |
---|
52 | ! mise a zero des tableaux |
---|
53 | ! ------------------------ |
---|
54 | tmp_dyn(:,:)=0.0 |
---|
55 | tmp_fi(:)=0.0 |
---|
56 | ! Lecture |
---|
57 | ! ----------------------- |
---|
58 | #ifdef NC_DOUBLE |
---|
59 | status=NF_GET_VARA_DOUBLE(ncid,varid,start,count,tmp_dyn) |
---|
60 | #else |
---|
61 | status=NF_GET_VARA_REAL(ncid,varid,start,count,tmp_dyn) |
---|
62 | #endif |
---|
63 | |
---|
64 | ! call dump2d(iip1,jjp1,tmp_dyn,'tmp_dyn ') |
---|
65 | DO j=1, jjp1 |
---|
66 | DO ig=1, iip1 |
---|
67 | tmp_dyn_invers(ig,j)=tmp_dyn(ig,jjp1-j+1) |
---|
68 | ENDDO |
---|
69 | ENDDO |
---|
70 | |
---|
71 | |
---|
72 | !JE20140522! call gr_dyn_fi_p(1, iip1, jjp1, klon, tmp_dyn_invers, tmp_fi) |
---|
73 | call gr_dyn_fi(1, iip1, jjp1, klon, tmp_dyn_invers, tmp_fi) |
---|
74 | ! call dump2d(iim,jjm-1,tmp_fi(2),'tmp_fi ') |
---|
75 | ! |
---|
76 | DO j=1,klon |
---|
77 | |
---|
78 | surfa(j,i)=tmp_fi(j) |
---|
79 | |
---|
80 | ENDDO ! Fin de recopie du tableau |
---|
81 | ! |
---|
82 | ENDDO ! Fin boucle 1 a 5 |
---|
83 | print*,'Passage Grille Dyn -> Phys' |
---|
84 | return |
---|
85 | end subroutine read_surface |
---|