Changeset 4356

Timestamp:

Nov 28, 2022, 6:27:07 PM (2 years ago)

Author:

asima

Message:

Bug fix: out_dim(2) instead of out_dim(3).
And a nice story:
The bug was there from "immemorial" times, older than the "Initial revision" on svn, 19 mai 2004 : LMDZ4/trunk/libf/dyn3d/grilles_gcm_netcdf.F@524).
It only started to act in rev4253, when replacing nf_def_var by nf90_def_var:
with the original "out_dim(3)", the 2D variables got a supplementary dimension "lonu" with NaN values, making grilles_gcm.nc impossible to open with cdo and ferret. But some regridding scripts using NCO and lat*&lon* variables could probably still work.
So the bug remained undetected until trying to use grille_s to regrid ERA 10m-winds for SPLA ("dust" module).

12. Fairhead (lfairhead) helped with the "bug hunt" and L. Guez (lguez) independently got the same solution.

File:

• LMDZ6/trunk/libf/dyn3d_common/grilles_gcm_netcdf_sub.F90 (modified) (1 diff)

LMDZ6/trunk/libf/dyn3d_common/grilles_gcm_netcdf_sub.F90

@@ -30 +30 @@
   INTEGER out_varid
   INTEGER out_lonudim,out_lonvdim
-  INTEGER out_latudim,out_latvdim,out_dim(3)
+  INTEGER out_latudim,out_latvdim,out_dim(2)
   INTEGER out_levdim