Newer
Older
!MNH_LIC Copyright 1996-2023 CNRS, Meteo-France and Universite Paul Sabatier
!MNH_LIC This is part of the Meso-NH software governed by the CeCILL-C licence
WAUTELET Philippe
committed
!MNH_LIC version 1. See LICENSE, CeCILL-C_V1-en.txt and CeCILL-C_V1-fr.txt
!MNH_LIC for details. version 1.
!-----------------------------------------------------------------
!######################
MODULE MODI_READ_EXSPA
!######################
!
INTERFACE
!
WAUTELET Philippe
committed
SUBROUTINE READ_EXSPA(HINIFILE,HINIFILEPGD, &
KXOR,KYOR,KXSIZE,KYSIZE,KDXRATIO,KDYRATIO,&
OBAL_ONLY,HDOMAIN,HSPAFILE,HSPANBR, &
HDADINIFILE,HDADSPAFILE,HSONFILE)
!
!* 0.1 Dummy arguments
!
CHARACTER (LEN=28), INTENT(OUT) :: HINIFILE ! Name of the model 1 FM file
CHARACTER (LEN=28), INTENT(OUT) :: HINIFILEPGD
INTEGER, INTENT(OUT) :: KXOR ! horizontal position (i,j) of the ORigin
INTEGER, INTENT(OUT) :: KYOR ! of the model 2 domain, relative to model 1
INTEGER, INTENT(OUT) :: KDXRATIO ! x and y-direction resolution RATIO
INTEGER, INTENT(OUT) :: KDYRATIO ! between model 2 and model 1
INTEGER, INTENT(OUT) :: KXSIZE ! number of model 1 grid points in x and y-directions
INTEGER, INTENT(OUT) :: KYSIZE ! in the model 2 physical domain
!
CHARACTER (LEN=28), INTENT(OUT) :: HDOMAIN ! input fm-file for grid definition
CHARACTER (LEN=28), INTENT(OUT) :: HSPAFILE ! possible name of the output FM-file
CHARACTER (LEN= 2), INTENT(OUT) :: HSPANBR ! NumBeR associated to the SPAwned file
CHARACTER (LEN=28), INTENT(OUT) :: HSONFILE ! Name of the SON 1 file
CHARACTER (LEN=28), INTENT(OUT) :: HDADINIFILE ! name of the initial dad file
WAUTELET Philippe
committed
! of the model 1 only for
! spawning 1 with GBAL_ONLY
CHARACTER (LEN=28), INTENT(OUT) :: HDADSPAFILE ! name of the dad file of the
! model 2 only for spawning 1
! with GBAL_ONLY
LOGICAL, INTENT(OUT) :: OBAL_ONLY ! compute anelastique balance
! without change in the file
! definition
!
END SUBROUTINE READ_EXSPA
!
END INTERFACE
!
END MODULE MODI_READ_EXSPA
!
!
WAUTELET Philippe
committed
! #################################################################
SUBROUTINE READ_EXSPA(HINIFILE,HINIFILEPGD, &
KXOR,KYOR,KXSIZE,KYSIZE,KDXRATIO,KDYRATIO,&
OBAL_ONLY,HDOMAIN,HSPAFILE,HSPANBR, &
HDADINIFILE,HDADSPAFILE,HSONFILE)
WAUTELET Philippe
committed
! #################################################################
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
!
!!**** *READ_EXSPA * - subroutine to read spawning namelist
!!
!! PURPOSE
!! -------
!!
!!
!!** METHOD
!! ------
!!
!!
!!
!! EXTERNAL
!! --------
!!
!!
!!
!! IMPLICIT ARGUMENTS
!! ------------------
!!
!!
!! REFERENCE
!! ---------
!!
!!
!!
!! AUTHOR
!! ------
!!
!! V. Masson * METEO-FRANCE *
!!
!! MODIFICATIONS
!! -------------
!!
!! Original 19/03/96
!! Modification 14/12/00 (P.Jabouille) do not close the EXSPA file
!! Modification 15/10/01 (I.Mallet) allow namelists in different orders
!! Modification 08/04/04 (G.Jaubert) spawning 1 for anelastic balance only
!! Modification 07/07/05 (D.Barbary) spawn with 2 input files (father+son1)
!! Modification 30/03/12 (S.Bielli) add NAM_NCOUT for netcdf output (removed 08/07/2016)
!! Modification 08/07/2016 (P.Wautelet) removed MNH_NCWRIT define
!! Philippe Wautelet: 05/2016-04/2018: new data structures and calls for I/O
WAUTELET Philippe
committed
! P. Wautelet 14/02/2019: remove CLUOUT/CLUOUT0 and associated variables
!
!-------------------------------------------------------------------------------
!
!* 0. DECLARATIONS
! ------------
!
USE MODD_CONF
WAUTELET Philippe
committed
USE MODD_IO, ONLY : TFILEDATA,TFILE_OUTPUTLISTING
WAUTELET Philippe
committed
USE MODD_LUNIT_n, ONLY : LUNIT_MODEL
WAUTELET Philippe
committed
USE MODE_IO, only: IO_Config_set
USE MODE_IO_FILE, only: IO_File_open, IO_File_close
USE MODE_IO_MANAGE_STRUCT, only: IO_File_add2list
USE MODE_POS
USE MODE_MODELN_HANDLER
!
WAUTELET Philippe
committed
USE MODN_CONFIO, ONLY : NAM_CONFIO
!
IMPLICIT NONE
!
!* 0.1 Dummy arguments
!
CHARACTER (LEN=28), INTENT(OUT) :: HINIFILE ! Name of the model 1 FM file
CHARACTER (LEN=28), INTENT(OUT) :: HINIFILEPGD
INTEGER, INTENT(OUT) :: KXOR ! horizontal position (i,j) of the ORigin
INTEGER, INTENT(OUT) :: KYOR ! of the model 2 domain, relative to model 1
INTEGER, INTENT(OUT) :: KDXRATIO ! x and y-direction resolution RATIO
INTEGER, INTENT(OUT) :: KDYRATIO ! between model 2 and model 1
INTEGER, INTENT(OUT) :: KXSIZE ! number of model 1 grid points in x and y-directions
INTEGER, INTENT(OUT) :: KYSIZE ! in the model 2 physical domain
!
CHARACTER (LEN=28), INTENT(OUT) :: HDOMAIN ! input fm-file for grid definition
CHARACTER (LEN=28), INTENT(OUT) :: HSPAFILE ! possible name of the output FM-file
CHARACTER (LEN= 2), INTENT(OUT) :: HSPANBR ! NumBeR associated to the SPAwned file
CHARACTER (LEN=28), INTENT(OUT) :: HSONFILE ! Name of the SON 1 file
CHARACTER (LEN=28), INTENT(OUT) :: HDADINIFILE ! name of the initial dad file
WAUTELET Philippe
committed
! of the model 1 only for
! spawning 1 with GBAL_ONLY
CHARACTER (LEN=28), INTENT(OUT) :: HDADSPAFILE ! name of the dad file of the
! model 2 only for spawning 1
! with GBAL_ONLY
LOGICAL, INTENT(OUT) :: OBAL_ONLY ! compute anelastique balance
! without change in the file
! definition
!
!
!* 0.2 Local variables
!
INTEGER :: IRESP ! Return codes in FM routines
INTEGER :: ILUOUT ! Logical unit number for the output listing
INTEGER :: ILUSPA ! Logical unit number for the EXSPA file
LOGICAL :: GFOUND ! Return code when searching namelist
!
!* prefixes in the namelists are such because of code history
!
INTEGER :: IXOR ! horizontal position (i,j) of the ORigin
INTEGER :: IYOR ! of the model 2 domain, relative to model 1
INTEGER :: IDXRATIO ! x and y-direction resolution RATIO
INTEGER :: IDYRATIO ! between model 2 and model 1
INTEGER :: IXSIZE,IYSIZE ! number of model 1 grid points in x and y-directions
! in the model 2 physical domain
LOGICAL :: GBAL_ONLY ! compute only anelastique balance
WAUTELET Philippe
committed
CHARACTER(LEN=28) :: YDOMAIN ! input fm-file for grid definition
CHARACTER(LEN=28) :: YSPAFILE ! possible name of the output FM-file
CHARACTER(LEN= 2) :: YSPANBR ! NumBeR associated to the SPAwned file
CHARACTER(LEN=28) :: YDADINIFILE ! Name of dad model for model 1
CHARACTER(LEN=28) :: YDADSPAFILE ! Name of dad model for model 2
CHARACTER(LEN=28) :: CINIFILE ! re-declaration because of namelist
CHARACTER(LEN=28) :: CINIFILEPGD ! re-declaration because of namelist
CHARACTER (LEN=28) :: YSONFILE = ' ' ! Name of SON input file
WAUTELET Philippe
committed
TYPE(TFILEDATA),POINTER :: TZNMLFILE ! Namelist file
171
172
173
174
175
176
177
178
179
180
181
182
183
184
185
186
187
188
189
190
191
192
193
194
195
196
!
!* 0.3 Namelist declarations
!
NAMELIST/NAM_GRID2_SPA/ IXOR,IYOR, &! horizontal position (i,j) of the origin
! of the model 2 domain, relative to model 1
IDXRATIO, &! x and y-direction resolution RATIO between
IDYRATIO, &! model 2 (to be spawned) and model 1
IXSIZE, &! number of grid meshes of model 1 covered
IYSIZE, &! by physical model 2 in each direction
GBAL_ONLY ! compute only anelastique balance
NAMELIST/NAM_LUNIT2_SPA/ CINIFILE, &! In file name (model 1)
CINIFILEPGD,&! pgd file associated with In file name (model 1)
YDOMAIN, &! fm file in which is read the grid
! (if used, previous namelist is obsolete)
YSPAFILE, &! Name proposed for spawned file
YSPANBR, &! NumBeR associated to the SPAwned file
! (used if YSPAFILE is not specified)
NVERB, &! verbosity flag
YDADINIFILE, & ! Name of dad model for model 1
YDADSPAFILE, & ! Name of dad model for model 2
YSONFILE ! Name of SON input file
!
!-------------------------------------------------------------------------------
!
!* 1. initialize logical unit number of the EXSPA file :
!
WAUTELET Philippe
committed
TZNMLFILE => NULL()
!
WAUTELET Philippe
committed
YDOMAIN = ' '
YSPAFILE = ' '
YSPANBR = '00'
YDADINIFILE = ' '
YDADSPAFILE = ' '
WAUTELET Philippe
committed
CALL IO_File_add2list(LUNIT_MODEL(2)%TLUOUT,'OUTPUT_LISTING2','OUTPUTLISTING','WRITE')
CALL IO_File_open(LUNIT_MODEL(2)%TLUOUT)
!
!Set output file for PRINT_MSG
TFILE_OUTPUTLISTING => LUNIT_MODEL(2)%TLUOUT
!
ILUOUT=LUNIT_MODEL(2)%TLUOUT%NLU
WAUTELET Philippe
committed
!
WAUTELET Philippe
committed
CALL IO_File_add2list(TZNMLFILE,'SPAWN1.nam','NML','READ')
CALL IO_File_open(TZNMLFILE)
WAUTELET Philippe
committed
ILUSPA = TZNMLFILE%NLU
!
!
!* 2. read the EXSPA file :
!
IDXRATIO=NUNDEF
IDYRATIO=NUNDEF
IXSIZE=NUNDEF
IYSIZE=NUNDEF
IXOR=NUNDEF
IYOR=NUNDEF
GBAL_ONLY=.FALSE.
!
CALL POSNAM( TZNMLFILE, 'NAM_GRID2_SPA', GFOUND )
IF (GFOUND) READ(ILUSPA,NAM_GRID2_SPA)
WAUTELET Philippe
committed
!
CINIFILE = LUNIT_MODEL(1)%CINIFILE !To respect default values
CINIFILEPGD = LUNIT_MODEL(1)%CINIFILEPGD !To respect default values
!
CALL POSNAM( TZNMLFILE, 'NAM_LUNIT2_SPA', GFOUND )
IF (GFOUND) READ(ILUSPA,NAM_LUNIT2_SPA)
WAUTELET Philippe
committed
LUNIT_MODEL(2)%CINIFILE = CINIFILE
LUNIT_MODEL(2)%CINIFILEPGD = CINIFILEPGD
!
CALL POSNAM( TZNMLFILE, 'NAM_CONFIO', GFOUND )
IF (GFOUND) READ(ILUSPA,NAM_CONFIO)
WAUTELET Philippe
committed
!
WAUTELET Philippe
committed
CALL IO_Config_set()
CALL IO_File_close(TZNMLFILE)
WAUTELET Philippe
committed
!* 3. model 1 and SON1 FM file name (passed as arguments)
WAUTELET Philippe
committed
HINIFILE = CINIFILE
HINIFILEPGD = CINIFILEPGD
HSONFILE = YSONFILE
!
!* 4. CINIFILE value is also used for model 2 (cf SPAWN_MODEL2)
!
!
!
!* 5. other file names
!
HDOMAIN = YDOMAIN
HSPAFILE = YSPAFILE
HSPANBR = YSPANBR
HDADINIFILE=YDADINIFILE
HDADSPAFILE=YDADSPAFILE
!
!* 6. zoom characteristics
!
KDXRATIO = IDXRATIO
KDYRATIO = IDYRATIO
KXSIZE = IXSIZE
KYSIZE = IYSIZE
KXOR = IXOR
KYOR = IYOR
!
OBAL_ONLY=GBAL_ONLY
!
WAUTELET Philippe
committed
!-------------------------------------------------------------------------------