Changeset 5072 for LMDZ6/trunk

Timestamp:

Jul 18, 2024, 3:12:27 PM (5 months ago)

Author:

abarral

Message:

Missing bits of r5071

File:

• LMDZ6/trunk/libf/dyn3dmem/guide_loc_mod.F90 (modified) (15 diffs)

LMDZ6/trunk/libf/dyn3dmem/guide_loc_mod.F90

@@ -9 +9 @@
   USE Write_Field_loc
   use lmdz_netcdf, only: nf90_nowrite, nf90_open, nf90_inq_varid, nf90_close, &
-          nf90_inq_dimid, nf90_inquire_dimension, nf_get_vara_rd, nf_inq_dimid, &
+          nf90_inq_dimid, nf90_inquire_dimension, nf_inq_dimid, &
           nf_inq_dimlen, nf_enddef, nf_def_dim, nf90_put_var, nf_noerr, nf_close, nf_inq_varid, &
           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
          IF (guide_plevs==1) THEN
-             status=nf_get_vara_rd(ncidpl,varidap,[1],[nlevnc],apnc)
-             status=nf_get_vara_rd(ncidpl,varidbp,[1],[nlevnc],bpnc)
+             status=nf90_put_var(ncidpl,varidap,apnc,[1],[nlevnc])
+             status=nf90_put_var(ncidpl,varidbp,bpnc,[1],[nlevnc])
          ELSEIF (guide_plevs==0) THEN
-             status=nf_get_vara_rd(ncidpl,varidpl,[1],[nlevnc],apnc)
+             status=nf90_put_var(ncidpl,varidpl,apnc,[1],[nlevnc])
 !FC Pour les corrections la pression est deja en Pascals on commente la ligne ci-dessous
              IF(convert_Pa) apnc=apnc*100.! conversion en Pascals
@@ -1815 +1815 @@
 ! Pression 
      if (guide_plevs==2) then
-         status=nf_get_vara_rd(ncidp,varidp,start,count,pnat2)
+         status=nf90_put_var(ncidp,varidp,pnat2,start,count)
          IF (invert_y) THEN
 !           PRINT*,"Invertion impossible actuellement"
@@ -1825 +1825 @@
 !  Vent zonal
      if (guide_u) then
-         status=nf_get_vara_rd(ncidu,varidu,start,count,unat2)
+         status=nf90_put_var(ncidu,varidu,unat2,start,count)
          IF (invert_y) THEN
 !           PRINT*,"Invertion impossible actuellement"
@@ -1837 +1837 @@
 !  Temperature
      if (guide_T) then
-         status=nf_get_vara_rd(ncidt,varidt,start,count,tnat2)
+         status=nf90_put_var(ncidt,varidt,tnat2,start,count)
          IF (invert_y) THEN
 !           PRINT*,"Invertion impossible actuellement"
@@ -1847 +1847 @@
 !  Humidite
      if (guide_Q) then
-         status=nf_get_vara_rd(ncidQ,varidQ,start,count,qnat2)
+         status=nf90_put_var(ncidQ,varidQ,qnat2,start,count)
          IF (invert_y) THEN
 !           PRINT*,"Invertion impossible actuellement"
@@ -1862 +1862 @@
          IF (invert_y) start(2)=jjm-jje_v+1
 
-         status=nf_get_vara_rd(ncidv,varidv,start,count,vnat2)
+         status=nf90_put_var(ncidv,varidv,vnat2,start,count)
          IF (invert_y) THEN
 !           PRINT*,"Invertion impossible actuellement"
@@ -1879 +1879 @@
          count(4)=0
          IF (invert_y) start(2)=jjp1-jje_u+1
-         status=nf_get_vara_rd(ncidps,varidps,start,count,psnat2)
+         status=nf90_put_var(ncidps,varidps,psnat2,start,count)
          IF (invert_y) THEN
 !           PRINT*,"Invertion impossible actuellement"
@@ -2039 +2039 @@
 ! Coefs ap, bp pour calcul de la pression aux differents niveaux
          if (guide_plevs==1) then
-             status=nf_get_vara_rd(ncidpl,varidap,[1],[nlevnc],apnc)
-             status=nf_get_vara_rd(ncidpl,varidbp,[1],[nlevnc],bpnc)
+             status=nf90_put_var(ncidpl,varidap,apnc,[1],[nlevnc])
+             status=nf90_put_var(ncidpl,varidbp,bpnc,[1],[nlevnc])
          elseif (guide_plevs==0) THEN
-             status=nf_get_vara_rd(ncidpl,varidpl,[1],[nlevnc],apnc)
+             status=nf90_put_var(ncidpl,varidpl,apnc,[1],[nlevnc])
              apnc=apnc*100.! conversion en Pascals
              bpnc(:)=0.
@@ -2067 +2067 @@
 !  Pression
      if (guide_plevs==2) then
-         status=nf_get_vara_rd(ncidp,varidp,start,count,zu)
+         status=nf90_put_var(ncidp,varidp,zu,start,count)
          DO i=1,iip1
              pnat2(i,:,:)=zu(:,:)
@@ -2080 +2080 @@
 !  Vent zonal
      if (guide_u) then
-         status=nf_get_vara_rd(ncidu,varidu,start,count,zu)
+         status=nf90_put_var(ncidu,varidu,zu,start,count)
          DO i=1,iip1
              unat2(i,:,:)=zu(:,:)
@@ -2095 +2095 @@
 !  Temperature
      if (guide_T) then
-         status=nf_get_vara_rd(ncidt,varidt,start,count,zu)
+         status=nf90_put_var(ncidt,varidt,zu,start,count)
          DO i=1,iip1
              tnat2(i,:,:)=zu(:,:)
@@ -2109 +2109 @@
 !  Humidite
      if (guide_Q) then
-         status=nf_get_vara_rd(ncidQ,varidQ,start,count,zu)
+         status=nf90_put_var(ncidQ,varidQ,zu,start,count)
          DO i=1,iip1
              qnat2(i,:,:)=zu(:,:)
@@ -2126 +2126 @@
          count(2)=jjnb_v
          IF (invert_y) start(2)=jjm-jje_v+1
-         status=nf_get_vara_rd(ncidv,varidv,start,count,zv)
+         status=nf90_put_var(ncidv,varidv,zv,start,count)
          DO i=1,iip1
              vnat2(i,:,:)=zv(:,:)
@@ -2148 +2148 @@
          count(4)=0
          IF (invert_y) start(2)=jjp1-jje_u+1
-         status=nf_get_vara_rd(ncidps,varidps,start,count,zu(:,1))
+         status=nf90_put_var(ncidps,varidps,zu(:,1),start,count)
          DO i=1,iip1
              psnat2(i,:)=zu(:,1)