Changeset 4179 for LMDZ6/branches/LMDZ_ECRad/libf/phylmd/regr_lat_time_coefoz_m.F90

Timestamp:

Jun 18, 2022, 1:32:25 AM (2 years ago)

Author:

lguez

Message:

Use an external updated NetCDF95 library

Remove NetCDF95 from source files. We want to use an up-to-date
NetCDF95 library to read a NetCDF file containing groups, for aerosol
optical properties. It seems complicated to keep the NetCDF95 library
inside LMDZ because:

• NetCDF95 now also needs a C compiler. I do not know how to make this work with FCM.

• NetCDF95 cannot be compiled with the -r8 option: some specific procedures in a generic interface become identical.

• Secondarily, we would have to change the names of files to adhere to the LMDZ standard. We are not glad to do that every time we update.

For now, we can compile using the options -include and -link of
makelmdz_fcm.

As we use an updated NetCDF95 library, we have to update some of the
calls in LMDZ. Those are the calls to nf95_inquire_variable and
nf95_gw_var which used to take a pointer argument and now take an
allocatable argument.

File:

• LMDZ6/branches/LMDZ_ECRad/libf/phylmd/regr_lat_time_coefoz_m.F90 (modified) (4 diffs)

LMDZ6/branches/LMDZ_ECRad/libf/phylmd/regr_lat_time_coefoz_m.F90

@@ -55 +55 @@
     integer n_lat! number of latitudes in the input data
 
-    real, pointer:: latitude(:)
+    real, allocatable:: latitude(:)
     ! (of input data, converted to rad, sorted in strictly ascending order)
 
@@ -62 +62 @@
     ! ascending order)
 
-    real, pointer:: plev(:) ! pressure level of input data
+    real, allocatable:: plev(:) ! pressure level of input data
     logical desc_lat ! latitude in descending order in the input file
 
@@ -172 +172 @@
     forall (j = 2:n_lat) lat_in_edg(j) = (latitude(j - 1) + latitude(j)) / 2
     lat_in_edg(n_lat + 1) = pi / 2
-    deallocate(latitude) ! pointer
+    deallocate(latitude)
 
     call nf95_inq_varid(ncid_in, "plev", varid)
@@ -192 +192 @@
     call nf95_put_var(ncid_out, varid_plev, plev)
 
-    deallocate(plev) ! pointer
+    deallocate(plev)
 
     allocate(o3_par_in(n_lat, n_plev, 12))