Changeset 5072
- Timestamp:
- Jul 18, 2024, 3:12:27 PM (7 weeks ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
LMDZ6/trunk/libf/dyn3dmem/guide_loc_mod.F90
r5071 r5072 9 9 USE Write_Field_loc 10 10 use lmdz_netcdf, only: nf90_nowrite, nf90_open, nf90_inq_varid, nf90_close, & 11 nf90_inq_dimid, nf90_inquire_dimension, nf_ get_vara_rd, nf_inq_dimid, &11 nf90_inq_dimid, nf90_inquire_dimension, nf_inq_dimid, & 12 12 nf_inq_dimlen, nf_enddef, nf_def_dim, nf90_put_var, nf_noerr, nf_close, nf_inq_varid, & 13 13 nf_redef, nf_write, nf_unlimited, nf_float, nf_clobber, nf_64bit_offset, nf90_float, & … … 1786 1786 ! Coefs ap, bp pour calcul de la pression aux differents niveaux 1787 1787 IF (guide_plevs==1) THEN 1788 status=nf _get_vara_rd(ncidpl,varidap,[1],[nlevnc],apnc)1789 status=nf _get_vara_rd(ncidpl,varidbp,[1],[nlevnc],bpnc)1788 status=nf90_put_var(ncidpl,varidap,apnc,[1],[nlevnc]) 1789 status=nf90_put_var(ncidpl,varidbp,bpnc,[1],[nlevnc]) 1790 1790 ELSEIF (guide_plevs==0) THEN 1791 status=nf _get_vara_rd(ncidpl,varidpl,[1],[nlevnc],apnc)1791 status=nf90_put_var(ncidpl,varidpl,apnc,[1],[nlevnc]) 1792 1792 !FC Pour les corrections la pression est deja en Pascals on commente la ligne ci-dessous 1793 1793 IF(convert_Pa) apnc=apnc*100.! conversion en Pascals … … 1815 1815 ! Pression 1816 1816 if (guide_plevs==2) then 1817 status=nf _get_vara_rd(ncidp,varidp,start,count,pnat2)1817 status=nf90_put_var(ncidp,varidp,pnat2,start,count) 1818 1818 IF (invert_y) THEN 1819 1819 ! PRINT*,"Invertion impossible actuellement" … … 1825 1825 ! Vent zonal 1826 1826 if (guide_u) then 1827 status=nf _get_vara_rd(ncidu,varidu,start,count,unat2)1827 status=nf90_put_var(ncidu,varidu,unat2,start,count) 1828 1828 IF (invert_y) THEN 1829 1829 ! PRINT*,"Invertion impossible actuellement" … … 1837 1837 ! Temperature 1838 1838 if (guide_T) then 1839 status=nf _get_vara_rd(ncidt,varidt,start,count,tnat2)1839 status=nf90_put_var(ncidt,varidt,tnat2,start,count) 1840 1840 IF (invert_y) THEN 1841 1841 ! PRINT*,"Invertion impossible actuellement" … … 1847 1847 ! Humidite 1848 1848 if (guide_Q) then 1849 status=nf _get_vara_rd(ncidQ,varidQ,start,count,qnat2)1849 status=nf90_put_var(ncidQ,varidQ,qnat2,start,count) 1850 1850 IF (invert_y) THEN 1851 1851 ! PRINT*,"Invertion impossible actuellement" … … 1862 1862 IF (invert_y) start(2)=jjm-jje_v+1 1863 1863 1864 status=nf _get_vara_rd(ncidv,varidv,start,count,vnat2)1864 status=nf90_put_var(ncidv,varidv,vnat2,start,count) 1865 1865 IF (invert_y) THEN 1866 1866 ! PRINT*,"Invertion impossible actuellement" … … 1879 1879 count(4)=0 1880 1880 IF (invert_y) start(2)=jjp1-jje_u+1 1881 status=nf _get_vara_rd(ncidps,varidps,start,count,psnat2)1881 status=nf90_put_var(ncidps,varidps,psnat2,start,count) 1882 1882 IF (invert_y) THEN 1883 1883 ! PRINT*,"Invertion impossible actuellement" … … 2039 2039 ! Coefs ap, bp pour calcul de la pression aux differents niveaux 2040 2040 if (guide_plevs==1) then 2041 status=nf _get_vara_rd(ncidpl,varidap,[1],[nlevnc],apnc)2042 status=nf _get_vara_rd(ncidpl,varidbp,[1],[nlevnc],bpnc)2041 status=nf90_put_var(ncidpl,varidap,apnc,[1],[nlevnc]) 2042 status=nf90_put_var(ncidpl,varidbp,bpnc,[1],[nlevnc]) 2043 2043 elseif (guide_plevs==0) THEN 2044 status=nf _get_vara_rd(ncidpl,varidpl,[1],[nlevnc],apnc)2044 status=nf90_put_var(ncidpl,varidpl,apnc,[1],[nlevnc]) 2045 2045 apnc=apnc*100.! conversion en Pascals 2046 2046 bpnc(:)=0. … … 2067 2067 ! Pression 2068 2068 if (guide_plevs==2) then 2069 status=nf _get_vara_rd(ncidp,varidp,start,count,zu)2069 status=nf90_put_var(ncidp,varidp,zu,start,count) 2070 2070 DO i=1,iip1 2071 2071 pnat2(i,:,:)=zu(:,:) … … 2080 2080 ! Vent zonal 2081 2081 if (guide_u) then 2082 status=nf _get_vara_rd(ncidu,varidu,start,count,zu)2082 status=nf90_put_var(ncidu,varidu,zu,start,count) 2083 2083 DO i=1,iip1 2084 2084 unat2(i,:,:)=zu(:,:) … … 2095 2095 ! Temperature 2096 2096 if (guide_T) then 2097 status=nf _get_vara_rd(ncidt,varidt,start,count,zu)2097 status=nf90_put_var(ncidt,varidt,zu,start,count) 2098 2098 DO i=1,iip1 2099 2099 tnat2(i,:,:)=zu(:,:) … … 2109 2109 ! Humidite 2110 2110 if (guide_Q) then 2111 status=nf _get_vara_rd(ncidQ,varidQ,start,count,zu)2111 status=nf90_put_var(ncidQ,varidQ,zu,start,count) 2112 2112 DO i=1,iip1 2113 2113 qnat2(i,:,:)=zu(:,:) … … 2126 2126 count(2)=jjnb_v 2127 2127 IF (invert_y) start(2)=jjm-jje_v+1 2128 status=nf _get_vara_rd(ncidv,varidv,start,count,zv)2128 status=nf90_put_var(ncidv,varidv,zv,start,count) 2129 2129 DO i=1,iip1 2130 2130 vnat2(i,:,:)=zv(:,:) … … 2148 2148 count(4)=0 2149 2149 IF (invert_y) start(2)=jjp1-jje_u+1 2150 status=nf _get_vara_rd(ncidps,varidps,start,count,zu(:,1))2150 status=nf90_put_var(ncidps,varidps,zu(:,1),start,count) 2151 2151 DO i=1,iip1 2152 2152 psnat2(i,:)=zu(:,1)
Note: See TracChangeset
for help on using the changeset viewer.