Changeset 5463 for LMDZ6/trunk/libf/phylmd/Dust/read_dust.f90

Timestamp:

Dec 31, 2024, 5:53:47 PM (5 days ago)

Author:

fhourdin

Message:

Bug fix for nf90_open in Dust routines

File:

• LMDZ6/trunk/libf/phylmd/Dust/read_dust.f90 (modified) (2 diffs)

LMDZ6/trunk/libf/phylmd/Dust/read_dust.f90

@@ -21 +21 @@
   save ncid1, varid1, ncid2, varid2
 !$OMP THREADPRIVATE(ncid1, varid1, ncid2, varid2)
-  integer :: start(4),count(4), status
+  integer :: start_(4),count_(4)
   integer :: i, j, ig
   !
@@ -28 +28 @@
   if (debutphy) then
   !
-     ncid1=nf90_open('dust.nc',nf90_nowrite,rcode)
-     varid1=nf90_inq_varid(ncid1,'EMISSION',rcode)
+     rcode=nf90_open('dust.nc',nf90_nowrite,ncid1)
+     if ( rcode /= 0 ) then ; call abort_physic('LMDZ','open dust.nc dans read_vent',1) ; endif
+
+     rcode=nf90_inq_varid(ncid1,'EMISSION',varid1)
+     if ( rcode /= 0 ) then ; call abort_physic('LMDZ','inq varid EMISSION dans read_vent',1) ; endif
   !
   endif
   !
-  start(1)=1
-  start(2)=1
-  start(4)=0
+  start_(1)=1
+  start_(2)=1
+  start_(3)=step
+  start_(4)=0
 
-   ! count(1)=iip1
-  count(1)=nbp_lon+1
-   ! count(2)=jjp1
-  count(2)=nbp_lat
-  count(3)=1
-  count(4)=0
+   ! count_(1)=iip1
+  count_(1)=nbp_lon+1
+   ! count_(2)=jjp1
+  count_(2)=nbp_lat
+  count_(3)=1
+  count_(4)=0
   !
-  start(3)=step
   !
-  status = nf90_get_var(ncid1, varid1, dust_nc_glo, start, count)
+  rcode = nf90_get_var(ncid1, varid1, dust_nc_glo, start_, count_)
+  if ( rcode /= 0 ) then ; call abort_physic('LMDZ','get EMISSION dans read_vent',1) ; endif
 
   !