phylmd/Interface_surf/read_sst.f

module read_sst_m

  implicit none

contains

  SUBROUTINE read_sst(julien, knindex, tsurf)

    ! From interfoce_lim

    use conf_gcm_m, only: lmt_pas
    USE dimphy, ONLY: klon
    USE netcdf, ONLY: nf90_nowrite
    use netcdf95, only: NF95_CLOSE, nf95_get_var, NF95_INQ_VARID, nf95_open
    use nr_util, only: assert
    use time_phylmdz, only: itap

    integer, intent(IN):: julien ! jour a lire dans l'annee

    integer, intent(in):: knindex(:) ! (knon)
    ! index des points de la surface a traiter

    real, intent(out):: tsurf(:) ! (knon)
    ! SST lues dans le fichier de conditions aux limites

    ! Local:

    ! Champ lu dans le fichier de conditions aux limites :
    real, save:: sst_lu(klon)

    integer ncid, varid ! pour NetCDF

    ! --------------------------------------------------

    call assert(size(knindex) == size(tsurf), "read_sst knon")

    ! Tester d'abord si c'est le moment de lire le fichier
    if (mod(itap - 1, lmt_pas) == 0) then
       call NF95_OPEN ('limit.nc', NF90_NOWRITE, ncid)

       call NF95_INQ_VARID(ncid, 'SST', varid)
       call NF95_GET_VAR(ncid, varid, sst_lu, start = (/1, julien/))

       call NF95_CLOSE(ncid)
    endif

    tsurf = sst_lu(knindex)

  END SUBROUTINE read_sst

end module read_sst_m
1	guez	202	module read_sst_m
2
3			implicit none
4
5			contains
6
7	guez	222	SUBROUTINE read_sst(julien, knindex, tsurf)
8	guez	202
9			! From interfoce_lim
10
11	guez	221	use conf_gcm_m, only: lmt_pas
12	guez	202	USE dimphy, ONLY: klon
13			USE netcdf, ONLY: nf90_nowrite
14			use netcdf95, only: NF95_CLOSE, nf95_get_var, NF95_INQ_VARID, nf95_open
15			use nr_util, only: assert
16			use time_phylmdz, only: itap
17
18	guez	221	integer, intent(IN):: julien ! jour a lire dans l'annee
19	guez	202
20			integer, intent(in):: knindex(:) ! (knon)
21			! index des points de la surface a traiter
22
23	guez	222	real, intent(out):: tsurf(:) ! (knon)
24	guez	202	! SST lues dans le fichier de conditions aux limites
25
26			! Local:
27
28			! Champ lu dans le fichier de conditions aux limites :
29	guez	221	real, save:: sst_lu(klon)
30	guez	202
31			integer ncid, varid ! pour NetCDF
32
33			! --------------------------------------------------
34
35	guez	222	call assert(size(knindex) == size(tsurf), "read_sst knon")
36	guez	202
37			! Tester d'abord si c'est le moment de lire le fichier
38	guez	221	if (mod(itap - 1, lmt_pas) == 0) then
39	guez	202	call NF95_OPEN ('limit.nc', NF90_NOWRITE, ncid)
40
41			call NF95_INQ_VARID(ncid, 'SST', varid)
42	guez	221	call NF95_GET_VAR(ncid, varid, sst_lu, start = (/1, julien/))
43	guez	202
44			call NF95_CLOSE(ncid)
45			endif
46
47	guez	222	tsurf = sst_lu(knindex)
48	guez	202
49			END SUBROUTINE read_sst
50
51			end module read_sst_m