Changeset 2408 for LMDZ5/branches/testing/libf/phylmd/readchlorophyll.F90

Timestamp:

Dec 14, 2015, 11:43:09 AM (9 years ago)

Author:

Laurent Fairhead

Message:

Merged trunk changes r2298:2396 into testing branch

Location:

LMDZ5/branches/testing

Files:

• . (modified) (1 prop)
• libf/phylmd/readchlorophyll.F90 (modified) (6 diffs)

LMDZ5/branches/testing/libf/phylmd/readchlorophyll.F90

@@ -10 +10 @@
 
     USE phys_cal_mod, ONLY : mth_cur
-    USE mod_grid_phy_lmdz
+    USE mod_grid_phy_lmdz, ONLY: nbp_lon, nbp_lat, klon_glo, &
+                                 grid2dto1d_glo
     USE mod_phys_lmdz_mpi_data, ONLY :  is_mpi_root
-    USE mod_phys_lmdz_para 
-    USE phys_state_var_mod
-    USE phys_local_var_mod
-    USE dimphy
+    USE mod_phys_lmdz_para, ONLY: scatter 
+    USE phys_state_var_mod, ONLY: chl_con
 
     implicit none
 
     include "YOMCST.h"
-    include "dimensions.h"
 
 ! Variable input
@@ -62 +60 @@
     n_lon = size(longitude)
 !    print *, 'LON chlorophyll=', n_lon, longitude
-    IF (n_lon.NE.iim) THEN
-       print *,'Le nombre de lon n est pas egal a iim'
+    IF (n_lon.NE.nbp_lon) THEN
+       print *,'Le nombre de lon n est pas egal a nbp_lon'
        STOP
     ENDIF
@@ -72 +70 @@
     n_lat = size(latitude)
 !    print *, 'LAT chlorophyll=', n_lat, latitude
-    IF (n_lat.NE.jjm+1) THEN 
-       print *,'Le nombre de lat n est pas egal a jjm+1'
+    IF (n_lat.NE.nbp_lat) THEN 
+       print *,'Le nombre de lat n est pas egal a jnbp_lat'
        STOP
     ENDIF
@@ -93 +91 @@
     CALL nf95_inq_varid(ncid_in, "CHL", varid)
     ncerr = nf90_get_var(ncid_in, varid, chlorocon)
-    print *,'code erreur readaerosolstrato=', ncerr, varid
+    print *,'code erreur readchlorophyll=', ncerr, varid
 
     CALL nf95_close(ncid_in)
@@ -99 +97 @@
 !---select the correct month
     IF (mth_cur.LT.1.OR.mth_cur.GT.12) THEN
-      print *,'probleme avec le mois dans readaerosolstrat =', mth_cur
+      print *,'probleme avec le mois dans readchlorophyll =', mth_cur
     ENDIF
     chlorocon_mois(:,:) = chlorocon(:,:,mth_cur)
@@ -109 +107 @@
     print*,"chrolophyll current month",mth_cur
     do i=1,klon_glo
-      if(isnan(chlorocon_mois_glo(i)))then
+!      if(isnan(chlorocon_mois_glo(i)))then ! isnan() is not in the Fortran standard...
+!      Another way to check for NaN:
+       if(chlorocon_mois_glo(i).ne.chlorocon_mois_glo(i)) then
          chlorocon_mois_glo(i)=0.
       endif