Changeset 5509 for LMDZ6/trunk/libf

Timestamp:

Jan 28, 2025, 5:21:18 PM (6 days ago)

Author:

fhourdin

Message:

Bug fix lecture netcdf nudging

File:

• LMDZ6/trunk/libf/dyn3dmem/guide_loc_mod.f90 (modified) (16 diffs)

LMDZ6/trunk/libf/dyn3dmem/guide_loc_mod.f90

@@ -1589 +1589 @@
 !=======================================================================
   SUBROUTINE guide_read(timestep)
-    USE netcdf, ONLY: nf90_put_var
+    USE netcdf, ONLY: nf90_get_var
     USE dimensions_mod, ONLY: iim, jjm, llm, ndm
 USE paramet_mod_h
@@ -1803 +1803 @@
          endif
 ! Coefs ap, bp pour calcul de la pression aux differents niveaux
+!function nf90_get_var(ncid, varid, values, start, count, stride, map)
+!  integer,                         intent( in) :: ncid, varid
+!  any valid type, scalar or array of any rank, &
+!                                   intent(out) :: values
+!  integer, dimension(:), optional, intent( in) :: start, count, stride, map
+!  integer                                      :: nf90_get_var
          IF (guide_plevs.EQ.1) THEN
-             status = nf90_put_var(ncidpl, varidap, apnc, [1], [nlevnc])
-             status = nf90_put_var(ncidpl, varidbp, bpnc, [1], [nlevnc])
+             status = nf90_get_var(ncidpl, varidap, apnc, [1], [nlevnc])
+             status = nf90_get_var(ncidpl, varidbp, bpnc, [1], [nlevnc])
          ELSEIF (guide_plevs.EQ.0) THEN
-             status = nf90_put_var(ncidpl, varidpl, apnc, [1], [nlevnc])
+             status = nf90_get_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
@@ -1833 +1839 @@
 ! Pression
      if (guide_plevs.EQ.2) then
-         status = nf90_put_var(ncidp, varidp, pnat2, start, count)
+         status = nf90_get_var(ncidp, varidp, pnat2, start, count)
          IF (invert_y) THEN
 !           PRINT*,"Invertion impossible actuellement"
@@ -1843 +1849 @@
 !  Vent zonal
      if (guide_u) then
-         status = nf90_put_var(ncidu, varidu, unat2, start, count)
+         status = nf90_get_var(ncidu, varidu, unat2, start, count)
          IF (invert_y) THEN
 !           PRINT*,"Invertion impossible actuellement"
@@ -1855 +1861 @@
 !  Temperature
      if (guide_T) then
-         status = nf90_put_var(ncidt, varidt, tnat2, start, count)
+         status = nf90_get_var(ncidt, varidt, tnat2, start, count)
          IF (invert_y) THEN
 !           PRINT*,"Invertion impossible actuellement"
@@ -1865 +1871 @@
 !  Humidite
      if (guide_Q) then
-         status = nf90_put_var(ncidQ, varidQ, qnat2, start, count)
+         status = nf90_get_var(ncidQ, varidQ, qnat2, start, count)
          IF (invert_y) THEN
 !           PRINT*,"Invertion impossible actuellement"
@@ -1879 +1885 @@
          count(2)=jjnb_v
          IF (invert_y) start(2)=jjm-jje_v+1
-         status = nf90_put_var(ncidv, varidv, vnat2, start, count)
+         status = nf90_get_var(ncidv, varidv, vnat2, start, count)
          IF (invert_y) THEN
 !           PRINT*,"Invertion impossible actuellement"
@@ -1896 +1902 @@
          count(4)=0
          IF (invert_y) start(2)=jjp1-jje_u+1
-         status = nf90_put_var(ncidps, varidps, psnat2, start, count)
+         status = nf90_get_var(ncidps, varidps, psnat2, start, count)
          IF (invert_y) THEN
 !           PRINT*,"Invertion impossible actuellement"
@@ -1908 +1914 @@
 !=======================================================================
   SUBROUTINE guide_read2D(timestep)
-    USE netcdf, ONLY: nf90_put_var
+    USE netcdf, ONLY: nf90_get_var
     USE dimensions_mod, ONLY: iim, jjm, llm, ndm
 USE paramet_mod_h
@@ -2057 +2063 @@
 ! Coefs ap, bp pour calcul de la pression aux differents niveaux
          if (guide_plevs.EQ.1) then
-             status = nf90_put_var(ncidpl, varidap, apnc, [1], [nlevnc])
-             status = nf90_put_var(ncidpl, varidbp, bpnc, [1], [nlevnc])
+             status = nf90_get_var(ncidpl, varidap, apnc, [1], [nlevnc])
+             status = nf90_get_var(ncidpl, varidbp, bpnc, [1], [nlevnc])
          elseif (guide_plevs.EQ.0) THEN
-             status = nf90_put_var(ncidpl, varidpl, apnc, [1], [nlevnc])
+             status = nf90_get_var(ncidpl, varidpl, apnc, [1], [nlevnc])
              apnc=apnc*100.! conversion en Pascals
              bpnc(:)=0.
@@ -2085 +2091 @@
 !  Pression
      if (guide_plevs.EQ.2) then
-         status = nf90_put_var(ncidp, varidp, zu, start, count)
+         status = nf90_get_var(ncidp, varidp, zu, start, count)
          DO i=1,iip1
              pnat2(i,:,:)=zu(:,:)
@@ -2098 +2104 @@
 !  Vent zonal
      if (guide_u) then
-         status = nf90_put_var(ncidu, varidu, zu, start, count)
+         status = nf90_get_var(ncidu, varidu, zu, start, count)
          DO i=1,iip1
              unat2(i,:,:)=zu(:,:)
@@ -2113 +2119 @@
 !  Temperature
      if (guide_T) then
-         status = nf90_put_var(ncidt, varidt, zu, start, count)
+         status = nf90_get_var(ncidt, varidt, zu, start, count)
          DO i=1,iip1
              tnat2(i,:,:)=zu(:,:)
@@ -2127 +2133 @@
 !  Humidite
      if (guide_Q) then
-         status = nf90_put_var(ncidQ, varidQ, zu, start, count)
+         status = nf90_get_var(ncidQ, varidQ, zu, start, count)
          DO i=1,iip1
              qnat2(i,:,:)=zu(:,:)
@@ -2144 +2150 @@
          count(2)=jjnb_v
          IF (invert_y) start(2)=jjm-jje_v+1
-         status = nf90_put_var(ncidv, varidv, zv, start, count)
+         status = nf90_get_var(ncidv, varidv, zv, start, count)
          DO i=1,iip1
              vnat2(i,:,:)=zv(:,:)
@@ -2166 +2172 @@
          count(4)=0
          IF (invert_y) start(2)=jjp1-jje_u+1
-         status = nf90_put_var(ncidps, varidps, zu(:, 1), start, count)
+         status = nf90_get_var(ncidps, varidps, zu(:, 1), start, count)
          DO i=1,iip1
              psnat2(i,:)=zu(:,1)