- Timestamp:
- Jul 22, 2024, 9:29:09 PM (2 months ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
LMDZ6/branches/Amaury_dev/libf/phylmd/Dust/read_surface.F90
r5088 r5099 5 5 ! ------ 6 6 USE ioipsl 7 ! USE comgeomphy8 7 USE dimphy 9 8 USE mod_grid_phy_lmdz 10 9 USE mod_phys_lmdz_para 11 10 USE iophy 12 USE netcdf, ONLY:nf90_inq_varid,nf90_noerr,nf90_get_var 11 USE netcdf, ONLY:nf90_inq_varid,nf90_noerr,nf90_get_var,nf90_nowrite,nf90_inq_varid,nf90_open 13 12 IMPLICIT NONE 14 13 … … 18 17 character*10 name 19 18 character*10 varname 20 ! 19 21 20 real tmp_dyn(iip1,jjp1) 22 21 real tmp_dyn_glo(nbp_lon+1,nbp_lat) 23 ! real tmp_dyn_glo(nbp_lon,nbp_lat)24 22 REAL tmp_dyn_invers(iip1,jjp1) 25 23 real tmp_dyn_invers_glo(nbp_lon+1,nbp_lat) 26 ! real tmp_dyn_invers_glo(nbp_lon,nbp_lat)27 24 real tmp_fi(klon) 28 25 real tmp_fi_glo(klon_glo) 29 26 real surfa(klon,5) 30 27 real surfa_glo(klon_glo,5) 31 ! 28 32 29 integer ncid 33 30 integer varid 34 realrcode31 integer rcode 35 32 integer start(2),count(2),status 36 33 integer i,j,l,ig … … 49 46 50 47 print*,'Lecture du fichier donnees_lisa.nc' 51 ncid= NCOPN('donnees_lisa.nc',NCNOWRIT,rcode)48 ncid=nf90_open('donnees_lisa.nc',nf90_nowrite,rcode) 52 49 53 50 !JE20140526<<: check if are inversed or not the latitude grid in donnes_lisa … … 72 69 ! endj(1)=jjp1 73 70 endj(1)=nbp_lat 74 varid= NCVID(ncid,latstr,rcode)71 varid=nf90_inq_varid(ncid,latstr,rcode) 75 72 76 73 status=nf90_get_var(ncid,varid,lats_glo,startj,endj) … … 89 86 varname=trim(name)//str1 90 87 print*,'lecture variable:',varname 91 varid=NCVID(ncid,trim(varname),rcode) 92 ! varid=NCVID(ncid,varname,rcode) 88 varid=nf90_inq_varid(ncid,trim(varname),rcode) 93 89 94 90 ! dimensions pour les champs scalaires et le vent zonal … … 135 131 !JE20140526>> 136 132 ! call dump2d(iim,jjm-1,tmp_fi(2),'tmp_fi ') 137 ! 133 138 134 DO j=1,klon_glo 139 135 … … 141 137 142 138 ENDDO ! Fin de recopie du tableau 143 ! 139 144 140 ENDDO ! Fin boucle 1 a 5 145 141 print*,'Passage Grille Dyn -> Phys'
Note: See TracChangeset
for help on using the changeset viewer.