- Timestamp:
- Jul 22, 2024, 9:29:09 PM (6 months ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
LMDZ6/branches/Amaury_dev/libf/phylmd/dyn1d/mod_1D_cases_read2.F90
r5088 r5099 1 ! 1 2 2 ! $Id: mod_1D_cases_read.F90 2373 2015-10-13 17:28:01Z jyg $ 3 ! 3 4 4 MODULE mod_1D_cases_read2 5 5 USE netcdf, ONLY: nf90_get_var,nf90_noerr,nf90_inq_varid,nf90_inquire_dimension,nf90_strerror,nf90_open,& … … 372 372 else 373 373 allocate(time_val(nt_cas)) 374 ierr = NF90_GET_VAR(nid,timeid,time_val)374 ierr = nf90_get_var(nid,timeid,time_val) 375 375 if(ierr/=nf90_noerr) then 376 376 print *,'Pb a la lecture de time cas.nc: ' … … 582 582 print *,'Dans read_cas2, on va lire ',var3didin(i),name_var(i) 583 583 if(i<=35) then 584 ierr = NF90_GET_VAR(nid,var3didin(i),resul, count = [1, 1, nlevel, ntime])584 ierr = nf90_get_var(nid,var3didin(i),resul, count = [1, 1, nlevel, ntime]) 585 585 print *,'Dans read_cas2, on a lu ',ierr,var3didin(i),name_var(i) 586 586 if(ierr/=nf90_noerr) then … … 590 590 else 591 591 print *,'Dans read_cas2, on a lu ',ierr,var3didin(i),name_var(i) 592 ierr = NF90_GET_VAR(nid,var3didin(i),resul1, count = [1, 1, ntime])592 ierr = nf90_get_var(nid,var3didin(i),resul1, count = [1, 1, ntime]) 593 593 if(ierr/=nf90_noerr) then 594 594 print *,'Pb a la lecture de cas.nc: ',name_var(i) … … 702 702 !----------------------------------------------------------------------- 703 703 if(i<=4) then ! Lecture des coord pression en (nlevelp1,lat,lon) 704 ierr = NF90_GET_VAR(nid,var3didin(i),apbp, count = [1, 1, nlevel + 1])704 ierr = nf90_get_var(nid,var3didin(i),apbp, count = [1, 1, nlevel + 1]) 705 705 print *,'read2_cas(apbp), on a lu ',i,name_var(i) 706 706 if(ierr/=nf90_noerr) then … … 710 710 !----------------------------------------------------------------------- 711 711 else if(i>4.and.i<=45) then ! Lecture des variables en (time,nlevel,lat,lon) 712 ierr = NF90_GET_VAR(nid,var3didin(i),resul, count = [1, 1, nlevel, ntime])712 ierr = nf90_get_var(nid,var3didin(i),resul, count = [1, 1, nlevel, ntime]) 713 713 print *,'read2_cas(resul), on a lu ',i,name_var(i) 714 714 if(ierr/=nf90_noerr) then … … 718 718 !----------------------------------------------------------------------- 719 719 else if (i>45.and.i<=51) then ! Lecture des variables en (time,lat,lon) 720 ierr = NF90_GET_VAR(nid,var3didin(i),resul2, count = [1, 1, ntime])720 ierr = nf90_get_var(nid,var3didin(i),resul2, count = [1, 1, ntime]) 721 721 print *,'read2_cas(resul2), on a lu ',i,name_var(i) 722 722 if(ierr/=nf90_noerr) then … … 726 726 !----------------------------------------------------------------------- 727 727 else ! Lecture des constantes (lat,lon) 728 ierr = NF90_GET_VAR(nid,var3didin(i),resul3)728 ierr = nf90_get_var(nid,var3didin(i),resul3) 729 729 print *,'read2_cas(resul3), on a lu ',i,name_var(i) 730 730 if(ierr/=nf90_noerr) then … … 877 877 !----------------------------------------------------------------------- 878 878 if(i<=4) then ! Lecture des coord pression en (nlevelp1,lat,lon) 879 ierr = NF90_GET_VAR(nid,var3didin(i),apbp)879 ierr = nf90_get_var(nid,var3didin(i),apbp) 880 880 print *,'read2_cas(apbp), on a lu ',i,name_var(i) 881 881 if(ierr/=nf90_noerr) then … … 885 885 !----------------------------------------------------------------------- 886 886 else if(i>4.and.i<=12) then ! Lecture des variables en (time,nlevel,lat,lon) 887 ierr = NF90_GET_VAR(nid,var3didin(i),resul1)887 ierr = nf90_get_var(nid,var3didin(i),resul1) 888 888 print *,'read2_cas(resul1), on a lu ',i,name_var(i) 889 889 if(ierr/=nf90_noerr) then … … 894 894 !----------------------------------------------------------------------- 895 895 else if(i>12.and.i<=54) then ! Lecture des variables en (time,nlevel,lat,lon) 896 ierr = NF90_GET_VAR(nid,var3didin(i),resul)896 ierr = nf90_get_var(nid,var3didin(i),resul) 897 897 print *,'read2_cas(resul), on a lu ',i,name_var(i) 898 898 if(ierr/=nf90_noerr) then … … 903 903 !----------------------------------------------------------------------- 904 904 else if (i>54.and.i<=65) then ! Lecture des variables en (time,lat,lon) 905 ierr = NF90_GET_VAR(nid,var3didin(i),resul2)905 ierr = nf90_get_var(nid,var3didin(i),resul2) 906 906 print *,'read2_cas(resul2), on a lu ',i,name_var(i) 907 907 if(ierr/=nf90_noerr) then … … 912 912 !----------------------------------------------------------------------- 913 913 else ! Lecture des constantes (lat,lon) 914 ierr = NF90_GET_VAR(nid,var3didin(i),resul3)914 ierr = nf90_get_var(nid,var3didin(i),resul3) 915 915 print *,'read2_cas(resul3), on a lu ',i,name_var(i) 916 916 if(ierr/=nf90_noerr) then … … 1042 1042 !--------------------------------------------------------------------------------------- 1043 1043 ! Time interpolation of a 2D field to the timestep corresponding to day 1044 ! 1044 1045 1045 ! day: current julian day (e.g. 717538.2) 1046 1046 ! day1: first day of the simulation … … 1235 1235 ,lat_cas,sens_cas,ustar_cas & 1236 1236 ,uw_cas,vw_cas,q1_cas,q2_cas,tke_cas & 1237 ! 1237 1238 1238 ,ts_prof_cas,plev_prof_cas,t_prof_cas,theta_prof_cas & 1239 1239 ,thv_prof_cas,thl_prof_cas,qv_prof_cas,ql_prof_cas,qi_prof_cas & … … 1251 1251 !--------------------------------------------------------------------------------------- 1252 1252 ! Time interpolation of a 2D field to the timestep corresponding to day 1253 ! 1253 1254 1254 ! day: current julian day (e.g. 717538.2) 1255 1255 ! day1: first day of the simulation
Note: See TracChangeset
for help on using the changeset viewer.