Changeset 5075 for LMDZ6/trunk/libf/phylmd/tracco2i_mod.F90

Timestamp:

Jul 19, 2024, 10:05:57 AM (4 months ago)

Author:

abarral

Message:

[continued & end] replace netcdf by lmdz_netcdf.F90 wrapper
"use netcdf" is now only used in lmdz_netcdf.F90 (except ecrad and obsolete/)
<include "netcdf.inc"> is now likewise only used in lmdz_netcdf.F90.

systematically specify explicitely <USE lmdz_netcdf, ONLY:> (probably left some missing, to correct later on)

Further replacement of nf_put_* by nf90_put_* (same for _get_)

[minor] replace deprecated boolean operators along the way

File:

• LMDZ6/trunk/libf/phylmd/tracco2i_mod.F90 (modified) (8 diffs)

LMDZ6/trunk/libf/phylmd/tracco2i_mod.F90

@@ -104 +104 @@
 
 ! Initialization of tr_seri(id_CO2) If it is not initialized
-      IF (MAXVAL(tr_seri(:,:,id_CO2)).LT.1.e-15) THEN 
+      IF (MAXVAL(tr_seri(:,:,id_CO2))<1.e-15) THEN
         tr_seri(:,:,id_CO2)=co2_ppm0*1.e-6/RMD*RMCO2 !--initialised from co2_ppm0 in rdem
       ENDIF 
@@ -299 +299 @@
 !--for every timestep comment out the IF ENDIF statements 
 !--otherwise this is updated every day
-    IF (debutphy.OR.day_cur.NE.day_pre) THEN
+    IF (debutphy.OR.day_cur/=day_pre) THEN
 
       CALL gather(tr_seri(:,:,id_CO2),co2_glo)
@@ -351 +351 @@
 
     USE netcdf95, ONLY: nf95_close, nf95_gw_var, nf95_inq_varid, nf95_open
-    USE netcdf, ONLY: nf90_get_var, nf90_noerr, nf90_nowrite
+    USE lmdz_netcdf, ONLY: nf90_get_var, nf90_noerr, nf90_nowrite
 
     USE carbon_cycle_mod, ONLY : fco2_ff, fco2_bb, fco2_land, fco2_ocean
@@ -401 +401 @@
         CALL nf95_gw_var(ncid_in, varid, vector)
         n_glo = size(vector)
-        IF (n_glo.NE.klon_glo) THEN
+        IF (n_glo/=klon_glo) THEN
            abort_message='sflx_lmdz_co2_ff: le nombre de points n est pas egal a klon_glo'
            CALL abort_physic(modname,abort_message,1)
@@ -409 +409 @@
         CALL nf95_gw_var(ncid_in, varid, time)
         n_month = size(time)
-        IF (n_month.NE.12) THEN
+        IF (n_month/=12) THEN
            abort_message='sflx_lmdz_co2_ff: le nombre de month n est pas egal a 12'
            CALL abort_physic(modname,abort_message,1)
@@ -434 +434 @@
       CALL nf95_gw_var(ncid_in, varid, vector)
       n_glo = size(vector)
-      IF (n_glo.NE.klon_glo) THEN
+      IF (n_glo/=klon_glo) THEN
          abort_message='sflx_lmdz_co2_bb: le nombre de points n est pas egal a klon_glo'
          CALL abort_physic(modname,abort_message,1)
@@ -442 +442 @@
       CALL nf95_gw_var(ncid_in, varid, time)
       n_month = size(time)
-      IF (n_month.NE.12) THEN
+      IF (n_month/=12) THEN
          abort_message='sflx_lmdz_co2_bb: le nombre de month n est pas egal a 12'
          CALL abort_physic(modname,abort_message,1)
@@ -474 +474 @@
 
 !---select the correct month
-  IF (mth_cur.LT.1.OR.mth_cur.GT.12) THEN
+  IF (mth_cur<1.OR.mth_cur>12) THEN
     PRINT *,'probleme avec le mois dans co2_ini =', mth_cur
   ENDIF