Changeset 3605 for LMDZ6/branches/Ocean_skin/libf/phylmd/readchlorophyll.F90

Timestamp:

Nov 21, 2019, 4:43:45 PM (4 years ago)

Author:

lguez

Message:

Merge revisions 3427:3600 of trunk into branch Ocean_skin

Location:

LMDZ6/branches/Ocean_skin

Files:

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

LMDZ6/branches/Ocean_skin/libf/phylmd/readchlorophyll.F90

@@ -15 +15 @@
     USE mod_phys_lmdz_para, ONLY: scatter 
     USE phys_state_var_mod, ONLY: chl_con
+    USE print_control_mod, ONLY: prt_level,lunout
 
     IMPLICIT NONE
@@ -45 +46 @@
 
 !--------------------------------------------------------
+    CHARACTER (len = 20)  :: modname = 'readchlorophyll'
+    CHARACTER (len = 80)  :: abort_message
 
 !--only read file if beginning of run or start of new month
@@ -56 +59 @@
     CALL nf95_gw_var(ncid_in, varid, longitude)
     n_lon = size(longitude)
-!    print *, 'LON chlorophyll=', n_lon, longitude
     IF (n_lon.NE.nbp_lon) THEN
-       print *,'Le nombre de lon n est pas egal a nbp_lon'
-       STOP
+       abort_message='Le nombre de lon n est pas egal a nbp_lon'
+       CALL abort_physic(modname,abort_message,1)
     ENDIF
 
@@ -65 +67 @@
     CALL nf95_gw_var(ncid_in, varid, latitude)
     n_lat = size(latitude)
-!    print *, 'LAT chlorophyll=', n_lat, latitude
     IF (n_lat.NE.nbp_lat) THEN 
-       print *,'Le nombre de lat n est pas egal a jnbp_lat'
-       STOP
+       abort_message='Le nombre de lat n est pas egal a jnbp_lat'
+       CALL abort_physic(modname,abort_message,1)
     ENDIF
 
@@ -74 +75 @@
     CALL nf95_gw_var(ncid_in, varid, time)
     n_month = size(time)
-!    print *, 'TIME aerosol strato=', n_month, time
     IF (n_month.NE.12) THEN 
-       print *,'Le nombre de month n est pas egal a 12'
-       STOP
+       abort_message='Le nombre de month n est pas egal a 12'
+       CALL abort_physic(modname,abort_message,1)
     ENDIF
 
@@ -87 +87 @@
     CALL nf95_inq_varid(ncid_in, "CHL", varid)
     ncerr = nf90_get_var(ncid_in, varid, chlorocon)
-    print *,'code erreur readchlorophyll=', ncerr, varid
+    WRITE(lunout,*)'code erreur readchlorophyll=', ncerr, varid
 
     CALL nf95_close(ncid_in)
@@ -93 +93 @@
 !---select the correct month
     IF (mth_cur.LT.1.OR.mth_cur.GT.12) THEN
-      print *,'probleme avec le mois dans readchlorophyll =', mth_cur
+      WRITE(lunout,*)'probleme avec le mois dans readchlorophyll =', mth_cur
     ENDIF
     chlorocon_mois(:,:) = chlorocon(:,:,mth_cur)
@@ -100 +100 @@
     CALL grid2dTo1d_glo(chlorocon_mois,chlorocon_mois_glo)
 
-    print *,"chrolophyll current month",mth_cur
+    WRITE(lunout,*)"chrolophyll current month",mth_cur
     DO i=1,klon_glo
 !      if(isnan(chlorocon_mois_glo(i)))then ! isnan() is not in the Fortran standard...