source: LMDZ6/branches/LMDZ_ECRad/libf/phylmd/readchlorophyll.F90 @ 4203

Last change on this file since 4203 was 4179, checked in by lguez, 2 years ago

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 size: 3.9 KB
Line 
1!
2! $Id$
3!
4!--This routine is to be tested with MPI / OMP parallelism
5!--OB 26/03/2018
6
7SUBROUTINE readchlorophyll(debut)
8
9    USE netcdf95, ONLY: nf95_close, nf95_gw_var, nf95_inq_dimid, nf95_inq_varid, nf95_open
10    USE netcdf, ONLY: nf90_get_var, nf90_noerr, nf90_nowrite
11    USE phys_cal_mod, ONLY: mth_cur
12    USE mod_grid_phy_lmdz, ONLY: nbp_lon, nbp_lat, klon_glo, grid2dto1d_glo
13    USE mod_phys_lmdz_mpi_data, ONLY: is_mpi_root
14    USE mod_phys_lmdz_omp_data, ONLY: is_omp_root
15    USE mod_phys_lmdz_para, ONLY: scatter
16    USE phys_state_var_mod, ONLY: chl_con
17    USE print_control_mod, ONLY: prt_level,lunout
18
19    IMPLICIT NONE
20
21    INCLUDE "YOMCST.h"
22
23! Variable input
24    LOGICAL debut
25
26! Variables locales
27    INTEGER n_lat   ! number of latitudes in the input data
28    INTEGER n_lon   ! number of longitudes in the input data
29    INTEGER n_lev   ! number of levels in the input data
30    INTEGER n_month ! number of months in the input data
31    REAL, ALLOCATABLE :: latitude(:)
32    REAL, ALLOCATABLE :: longitude(:)
33    REAL, ALLOCATABLE :: time(:)
34    INTEGER i, k
35    INTEGER, SAVE :: mth_pre
36!$OMP THREADPRIVATE(mth_pre)
37
38! Champs reconstitues
39    REAL, ALLOCATABLE :: chlorocon(:, :, :)
40    REAL, ALLOCATABLE :: chlorocon_mois(:, :)
41    REAL, ALLOCATABLE :: chlorocon_mois_glo(:)
42
43! For NetCDF:
44    INTEGER ncid_in  ! IDs for input files
45    INTEGER varid, ncerr
46
47!--------------------------------------------------------
48    CHARACTER (len = 20)  :: modname = 'readchlorophyll'
49    CHARACTER (len = 80)  :: abort_message
50
51!--only read file if beginning of run or start of new month
52    IF (debut.OR.mth_cur.NE.mth_pre) THEN
53
54    IF (is_mpi_root.AND.is_omp_root) THEN
55
56    CALL nf95_open("chlorophyll.nc", nf90_nowrite, ncid_in)
57
58    CALL nf95_inq_varid(ncid_in, "lon", varid)
59    CALL nf95_gw_var(ncid_in, varid, longitude)
60    n_lon = size(longitude)
61    IF (n_lon.NE.nbp_lon) THEN
62       abort_message='Le nombre de lon n est pas egal a nbp_lon'
63       CALL abort_physic(modname,abort_message,1)
64    ENDIF
65
66    CALL nf95_inq_varid(ncid_in, "lat", varid)
67    CALL nf95_gw_var(ncid_in, varid, latitude)
68    n_lat = size(latitude)
69    IF (n_lat.NE.nbp_lat) THEN
70       abort_message='Le nombre de lat n est pas egal a jnbp_lat'
71       CALL abort_physic(modname,abort_message,1)
72    ENDIF
73
74    CALL nf95_inq_varid(ncid_in, "time", varid)
75    CALL nf95_gw_var(ncid_in, varid, time)
76    n_month = size(time)
77    IF (n_month.NE.12) THEN
78       abort_message='Le nombre de month n est pas egal a 12'
79       CALL abort_physic(modname,abort_message,1)
80    ENDIF
81
82    IF (.not.ALLOCATED(chlorocon))          ALLOCATE(chlorocon(n_lon, n_lat, n_month))
83    IF (.not.ALLOCATED(chlorocon_mois))     ALLOCATE(chlorocon_mois(n_lon, n_lat))
84    IF (.not.ALLOCATED(chlorocon_mois_glo)) ALLOCATE(chlorocon_mois_glo(klon_glo))
85
86!--reading stratospheric AOD at 550 nm
87    CALL nf95_inq_varid(ncid_in, "CHL", varid)
88    ncerr = nf90_get_var(ncid_in, varid, chlorocon)
89    WRITE(lunout,*)'code erreur readchlorophyll=', ncerr, varid
90
91    CALL nf95_close(ncid_in)
92
93!---select the correct month
94    IF (mth_cur.LT.1.OR.mth_cur.GT.12) THEN
95      WRITE(lunout,*)'probleme avec le mois dans readchlorophyll =', mth_cur
96    ENDIF
97    chlorocon_mois(:,:) = chlorocon(:,:,mth_cur)
98
99!---reduce to a klon_glo grid
100    CALL grid2dTo1d_glo(chlorocon_mois,chlorocon_mois_glo)
101
102    WRITE(lunout,*)"chrolophyll current month",mth_cur
103    DO i=1,klon_glo
104!      if(isnan(chlorocon_mois_glo(i)))then ! isnan() is not in the Fortran standard...
105!      Another way to check for NaN:
106       IF (chlorocon_mois_glo(i).NE.chlorocon_mois_glo(i)) chlorocon_mois_glo(i)=0.
107    ENDDO
108
109!    DEALLOCATE(chlorocon)
110!    DEALLOCATE(chlorocon_mois)
111!    DEALLOCATE(chlorocon_mois_glo)
112 
113    ENDIF !--is_mpi_root and is_omp_root
114
115!--scatter on all proc
116    CALL scatter(chlorocon_mois_glo,chl_con)
117
118!--keep memory of previous month
119    mth_pre=mth_cur
120
121    ENDIF !--debut ou nouveau mois
122
123END SUBROUTINE readchlorophyll
Note: See TracBrowser for help on using the repository browser.