Newer
Older
!MNH_LIC Copyright 1994-2018 CNRS, Meteo-France and Universite Paul Sabatier
!MNH_LIC This is part of the Meso-NH software governed by the CeCILL-C licence
!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_WRITE_DESFM_n
! #########################
!
INTERFACE
!
WAUTELET Philippe
committed
SUBROUTINE WRITE_DESFM_n(KMI,TPDATAFILE)
!
USE MODD_IO_ll, ONLY: TFILEDATA
!
INTEGER, INTENT(IN) :: KMI ! Model index
TYPE(TFILEDATA), INTENT(IN) :: TPDATAFILE ! Datafile
!
END SUBROUTINE WRITE_DESFM_n
!
END INTERFACE
!
END MODULE MODI_WRITE_DESFM_n
!
!
! ###################################################
WAUTELET Philippe
committed
SUBROUTINE WRITE_DESFM_n(KMI,TPDATAFILE)
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
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
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
129
130
131
132
133
! ###################################################
!
!!**** *WRITE_DESFM_n * - routine to write a descriptor file ( DESFM )
!!
!! PURPOSE
!! -------
! The purpose of this routine is to write the descriptive part of a Mesonh
! file (FM-file). The resulting file is called DESFM.
!
!!
!!** METHOD
!! ------
!!
!! This routine writes in the file HDESFM, previously opened, the group of
!! all the namelists used to specify a Mesonh simulation.
!! If verbose option is high enough : NVERB>=5, the variables in descriptor
!! file are printed on the right output-listing corresponding tomodel _n.
!!
!! EXTERNAL
!! --------
!!
!!
!! IMPLICIT ARGUMENTS
!! ------------------
!!
!! Module MODN_LUNIT_n : contains declarations of namelist NAM_LUNITn
!! and module MODD_LUNIT_n
!!
!!
!! Module MODN_CONF_n : contains declaration of namelist NAM_CONFn and
!! uses module MODD_CONF1 (configuration variables
!! for model _n )
!!
!! Module MODN_DYN_n : contains declaration of namelist NAM_DYNn and
!! uses module MODD_DYN_n (dynamic control variables
!! for model _n )
!!
!! Module MODN_ADV_n : contains declaration of namelist NAM_ADVn and
!! uses module MODD_ADV_n (control variables for the
!! advection scheme for model _n )
!!
!! Module MODN_PARAM_n : contains declaration of namelist NAM_PARAMn and
!! uses module MODD_PARAM_n (names of the physical
!! parameterizations for model _n )
!!
!! Module MODN_PARAM_RAD_n : contains declaration of the control parameters
!! for calling the radiation scheme
!!
!! Module MODN_PARAM_KAFR_n : contains declaration of control parameters
!! for calling the deep convection scheme
!!
!! Module MODN_LBC_n : contains declaration of namelis NAM_LBCn and
!! uses module MODD_LBC_n (lateral boundary conditions)
!!
!!
!! Module MODN_TURB_n : contains declaration of turbulence scheme options
!! present in the namelist
!!
!! Module MODN_CONF : contains declaration of namelist NAM_CONF and
!! uses module MODD_CONF (configuration variables)
!!
!! Module MODN_DYN : contains the declaration of namelist NAM_DYN and
!! uses module MODD_DYN (dynamic control variables)
!!
!! Module MODN_BUDGET : contains declaration of all the namelists
!! related to the budget computations
!!
!! Module MODN_LES : contains declaration of the control parameters
!! for Large Eddy Simulations' storages
!! Module MODN_BLANK: contains declaration of MesoNH developper variables
!! for test and debugging purposes.
!!
!!
!! REFERENCE
!! ---------
!! None
!!
!!
!! AUTHOR
!! ------
!! V. Ducrocq * Meteo France *
!!
!! MODIFICATIONS
!! -------------
!! Original 07/06/94
!! Updated V.Ducrocq 06/09/94
!! Updated J.Stein 20/10/94 to include NAM_OUTn
!! Updated J.Stein 24/10/94 change routine name
!! Updated J.Stein 26/10/94 add the OWRIGET argument
!! Updated J.Stein 06/12/94 add the LS fields
!! Updated J.Stein 09/01/95 add the turbulence scheme
!! Updated J.Stein 09/01/95 add the 1D switch
!! Updated J.Stein 20/03/95 remove R from the historical var.
!! Updated Ph.Hereil 20/06/95 add the budgets
!! Updated J.-P. Pinty 15/09/95 add the radiations
!! Updated J.Vila 06/02/96 implementation of scalar
!! advection schemes
!! Updated J.Stein 20/02/96 cleaning + add the LES namelist
!! Modifications 25/04/96 (Suhre) add NAM_BLANK
!! Modifications 25/04/96 (Suhre) add NAM_FRC
!! Modifications 25/04/96 (Suhre) add NAM_CH_MNHCn and NAM_CH_SOLVER
!! Modifications 11/04/96 (Pinty) add the ice concentration
!! Modifications 11/01/97 (Pinty) add the deep convection
!! Temporary Modification (Masson 06/09/96) manual write of the first and
!! third namelists because of compiler version.
!! Modifications J.-P. Lafore 22/07/96 gridnesting implementation
!! Modifications J.-P. Lafore 29/07/96 add NAM_FMOUT (renamed in NAM_OUTPUT/NAM_BACKUP)
!! Modifications V. Masson 10/07/97 add NAM_PARAM_GROUNDn
!! Modifications V. Masson 28/07/97 supress LSTEADY_DMASS
!! Modifications P. Jabouille 03/10/01 LHORELAX_ modifications
!! Modifications P. Jabouille 12/03/02 conditional writing of namelists
!! Modifications J.-P. Pinty 29/11/02 add C3R5, ICE2, ICE4, CELEC
!! Modification V. Masson 01/2004 removes surface (externalization)
!! Modification P. Tulet 01/2005 add dust, orilam
!! Modification 05/2006 Remove EPS and OWRIGET
!! 02/2018 Q.Libois ECRAD
!! Philippe Wautelet: 05/2016-04/2018: new data structures and calls for I/O
!-------------------------------------------------------------------------------
!
!* 0. DECLARATIONS
! ------------
USE MODD_CONF
WAUTELET Philippe
committed
USE MODD_DYN_n, ONLY: LHORELAX_SVLIMA
USE MODD_IO_ll, ONLY: TFILEDATA
USE MODD_LUNIT_n, ONLY: TLUOUT
USE MODD_PARAMETERS
!
USE MODE_MSG
USE MODN_BACKUP
USE MODN_CONF
USE MODN_DYN
USE MODN_NESTING
WAUTELET Philippe
committed
USE MODN_OUTPUT
USE MODN_BUDGET
USE MODN_LES
USE MODN_DYN_n
USE MODN_ADV_n
USE MODN_PARAM_n
USE MODN_PARAM_RAD_n
USE MODN_PARAM_ECRAD_n
USE MODN_PARAM_KAFR_n
USE MODN_PARAM_MFSHALL_n
USE MODN_PARAM_ICE
USE MODN_CONF_n
USE MODN_LUNIT_n
USE MODN_LBC_n
USE MODN_NUDGING_n
USE MODN_TURB_n
USE MODN_BLANK
USE MODN_FRC
USE MODN_CH_MNHC_n
USE MODN_CH_SOLVER_n
USE MODN_PARAM_C2R2
USE MODN_PARAM_C1R3
USE MODN_ELEC
USE MODN_SERIES
USE MODN_SERIES_n
USE MODN_TURB_CLOUD
USE MODN_TURB
USE MODN_CH_ORILAM
USE MODN_DUST
USE MODN_SALT
USE MODN_PASPOL
USE MODN_CONDSAMP
USE MODN_2D_FRC
USE MODN_LATZ_EDFLX
#ifdef MNH_FOREFIRE
USE MODN_FOREFIRE
USE MODD_FOREFIRE_n, ONLY : FFCOUPLING
#endif
!
IMPLICIT NONE
!
!* 0.1 declarations of arguments
!
WAUTELET Philippe
committed
INTEGER, INTENT(IN) :: KMI ! Model index
TYPE(TFILEDATA), INTENT(IN) :: TPDATAFILE ! Datafile
!
!* 0.2 declarations of local variables
!
INTEGER :: ILUSEG ! logical unit number of EXSEG file
INTEGER :: ILUOUT ! Logical unit number for output-listing TLUOUT file
!
LOGICAL :: GHORELAX_UVWTH, &
GHORELAX_RV, GHORELAX_RC, GHORELAX_RR, &
GHORELAX_RI, GHORELAX_RS, GHORELAX_RG, &
GHORELAX_TKE, GHORELAX_SVC2R2, GHORELAX_SVPP, &
GHORELAX_SVCS, GHORELAX_SVCHIC, &
#ifdef MNH_FOREFIRE
GHORELAX_SVFF, &
#endif
GHORELAX_SVCHEM, GHORELAX_SVC1R3, &
GHORELAX_SVELEC, GHORELAX_SVLIMA
LOGICAL :: GHORELAX_SVDST, GHORELAX_SVSLT, GHORELAX_SVAER
LOGICAL, DIMENSION(JPSVMAX) :: GHORELAX_SV
!
!-------------------------------------------------------------------------------
!
!* 1. UPDATE DESFM FILE
! -----------------
!
WAUTELET Philippe
committed
CALL PRINT_MSG(NVERB_DEBUG,'IO','WRITE_DESFM_n','called for '//TRIM(TPDATAFILE%CNAME))
!
IF (.NOT.ASSOCIATED(TPDATAFILE%TDESFILE)) &
CALL PRINT_MSG(NVERB_FATAL,'IO','WRITE_DESFM_n','TDESFILE not associated for '//TRIM(TPDATAFILE%CNAME))
WAUTELET Philippe
committed
ILUSEG = TPDATAFILE%TDESFILE%NLU
!
CALL INIT_NAM_LUNITn
WRITE(UNIT=ILUSEG,NML=NAM_LUNITn)
IF (CPROGRAM/='MESONH') THEN
LUSECI=.FALSE.
ENDIF
CALL INIT_NAM_CONFn
WRITE(UNIT=ILUSEG,NML=NAM_CONFn)
!
!
CALL INIT_NAM_DYNn
IF (CPROGRAM/='MESONH') THEN ! impose default value for next simulation
GHORELAX_UVWTH = LHORELAX_UVWTH
GHORELAX_RV = LHORELAX_RV
GHORELAX_RC = LHORELAX_RC
GHORELAX_RR = LHORELAX_RR
GHORELAX_RI = LHORELAX_RI
GHORELAX_RS = LHORELAX_RS
GHORELAX_RG = LHORELAX_RG
GHORELAX_TKE = LHORELAX_TKE
GHORELAX_SV(:) = LHORELAX_SV(:)
GHORELAX_SVC2R2= LHORELAX_SVC2R2
GHORELAX_SVC1R3= LHORELAX_SVC1R3
GHORELAX_SVELEC= LHORELAX_SVELEC
GHORELAX_SVCHEM= LHORELAX_SVCHEM
GHORELAX_SVCHIC= LHORELAX_SVCHIC
GHORELAX_SVDST = LHORELAX_SVDST
GHORELAX_SVSLT = LHORELAX_SVSLT
GHORELAX_SVPP = LHORELAX_SVPP
#ifdef MNH_FOREFIRE
GHORELAX_SVFF = LHORELAX_SVFF
#endif
GHORELAX_SVCS = LHORELAX_SVCS
GHORELAX_SVAER = LHORELAX_SVAER
!
LHORELAX_UVWTH = .FALSE.
LHORELAX_RV = .FALSE.
LHORELAX_RC = .FALSE.
LHORELAX_RR = .FALSE.
LHORELAX_RI = .FALSE.
LHORELAX_RS = .FALSE.
LHORELAX_RG = .FALSE.
LHORELAX_TKE = .FALSE.
LHORELAX_SV(:) = .FALSE.
LHORELAX_SVC2R2= .FALSE.
LHORELAX_SVC1R3= .FALSE.
LHORELAX_SVELEC= .FALSE.
LHORELAX_SVCHEM= .FALSE.
LHORELAX_SVCHIC= .FALSE.
LHORELAX_SVLG = .FALSE.
LHORELAX_SVPP = .FALSE.
#ifdef MNH_FOREFIRE
LHORELAX_SVFF = .FALSE.
#endif
LHORELAX_SVCS = .FALSE.
LHORELAX_SVDST= .FALSE.
LHORELAX_SVSLT= .FALSE.
LHORELAX_SVAER= .FALSE.
ELSE !return to namelist meaning of LHORELAX_SV
GHORELAX_SV(:) = LHORELAX_SV(:)
LHORELAX_SV(NSV_USER+1:)=.FALSE.
END IF
WRITE(UNIT=ILUSEG,NML=NAM_DYNn)
!
CALL INIT_NAM_ADVn
WRITE(UNIT=ILUSEG,NML=NAM_ADVn)
IF (CPROGRAM/='MESONH') THEN
CTURB = 'NONE'
CRAD = 'NONE'
CCLOUD = 'NONE'
CDCONV = 'NONE'
CSCONV = 'NONE'
CELEC = 'NONE'
CACTCCN = 'NONE'
END IF
CALL INIT_NAM_PARAMn
WRITE(UNIT=ILUSEG,NML=NAM_PARAMn)
!
CALL INIT_NAM_PARAM_RADn
IF(CRAD /= 'NONE') WRITE(UNIT=ILUSEG,NML=NAM_PARAM_RADn)
#ifdef MNH_ECRAD
CALL INIT_NAM_PARAM_ECRADn
IF(CRAD /= 'NONE') WRITE(UNIT=ILUSEG,NML=NAM_PARAM_ECRADn)
#endif
!
CALL INIT_NAM_PARAM_KAFRn
IF(CDCONV /= 'NONE' .OR. CSCONV == 'KAFR') &
WRITE(UNIT=ILUSEG,NML=NAM_PARAM_KAFRn)
!
CALL INIT_NAM_PARAM_MFSHALLn
IF (CSCONV == 'EDKF' ) WRITE(UNIT=ILUSEG,NML=NAM_PARAM_MFSHALLn)
!
CALL INIT_NAM_LBCn
WRITE(UNIT=ILUSEG,NML=NAM_LBCn)
!
CALL INIT_NAM_NUDGINGn
WRITE(UNIT=ILUSEG,NML=NAM_NUDGINGn)
!
CALL INIT_NAM_TURBn
IF(CTURB /= 'NONE') WRITE(UNIT=ILUSEG,NML=NAM_TURBn)
!
!IF (CPROGRAM/='MESONH') THEN
! LUSECHEM = .FALSE.
! LORILAM = .FALSE.
! LDEPOS_AER = .FALSE.
! LDUST = .FALSE.
! LDEPOS_DST = .FALSE.
! LSALT = .FALSE.
! LDEPOS_SLT = .FALSE.
! LPASPOL = .FALSE.
! LCONDSAMP = .FALSE.
!END IF
CALL INIT_NAM_CH_MNHCn
IF(LUSECHEM .OR. LCH_CONV_LINOX .OR. LCH_CONV_SCAV) &
WRITE(UNIT=ILUSEG,NML=NAM_CH_MNHCn)
!
CALL INIT_NAM_CH_SOLVERn
IF(LUSECHEM) WRITE(UNIT=ILUSEG,NML=NAM_CH_SOLVERn)
!
IF(LDUST) WRITE(UNIT=ILUSEG,NML=NAM_DUST)
IF(LSALT) WRITE(UNIT=ILUSEG,NML=NAM_SALT)
IF(LPASPOL) WRITE(UNIT=ILUSEG,NML=NAM_PASPOL)
#ifdef MNH_FOREFIRE
IF(FFCOUPLING) WRITE(UNIT=ILUSEG,NML=NAM_FOREFIRE)
#endif
IF(LCONDSAMP) WRITE(UNIT=ILUSEG,NML=NAM_CONDSAMP)
IF(LORILAM.AND.LUSECHEM) WRITE(UNIT=ILUSEG,NML=NAM_CH_ORILAM)
!
CALL INIT_NAM_SERIESn
IF(LSERIES) WRITE(UNIT=ILUSEG,NML=NAM_SERIESn)
IF(L2D_ADV_FRC .OR. L2D_REL_FRC) WRITE(UNIT=ILUSEG,NML=NAM_2D_FRC)
!
IF (LUV_FLX .OR. LTH_FLX) WRITE(UNIT=ILUSEG,NML=NAM_LATZ_EDFLX)
!
IF (CPROGRAM/='MESONH') THEN
LLG = .FALSE.
END IF
WRITE(UNIT=ILUSEG,NML=NAM_CONF)
WRITE(UNIT=ILUSEG,NML=NAM_DYN)
WRITE(UNIT=ILUSEG,NML=NAM_NESTING)
!WRITE(UNIT=ILUSEG,NML=NAM_BACKUP)
WAUTELET Philippe
committed
!WRITE(UNIT=ILUSEG,NML=NAM_OUTPUT)
IF(CBUTYPE /= 'NONE') THEN
IF(CBUTYPE=='SKIP') CBUTYPE='CART'
WRITE(UNIT=ILUSEG,NML=NAM_BUDGET)
END IF
IF(LBU_RU) WRITE(UNIT=ILUSEG,NML=NAM_BU_RU)
IF(LBU_RV) WRITE(UNIT=ILUSEG,NML=NAM_BU_RV)
IF(LBU_RW) WRITE(UNIT=ILUSEG,NML=NAM_BU_RW)
IF(LBU_RTH) WRITE(UNIT=ILUSEG,NML=NAM_BU_RTH)
IF(LBU_RTKE) WRITE(UNIT=ILUSEG,NML=NAM_BU_RTKE)
IF(LBU_RRV) WRITE(UNIT=ILUSEG,NML=NAM_BU_RRV)
IF(LBU_RRC) WRITE(UNIT=ILUSEG,NML=NAM_BU_RRC)
IF(LBU_RRR) WRITE(UNIT=ILUSEG,NML=NAM_BU_RRR)
IF(LBU_RRI) WRITE(UNIT=ILUSEG,NML=NAM_BU_RRI)
IF(LBU_RRS) WRITE(UNIT=ILUSEG,NML=NAM_BU_RRS)
IF(LBU_RRG) WRITE(UNIT=ILUSEG,NML=NAM_BU_RRG)
IF(LBU_RRH) WRITE(UNIT=ILUSEG,NML=NAM_BU_RRH)
IF(LBU_RSV) WRITE(UNIT=ILUSEG,NML=NAM_BU_RSV)
IF(LLES_MEAN .OR. LLES_RESOLVED .OR. LLES_SUBGRID .OR. LLES_UPDRAFT &
.OR. LLES_DOWNDRAFT .OR. LLES_SPECTRA) WRITE(UNIT=ILUSEG,NML=NAM_LES)
WRITE(UNIT=ILUSEG,NML=NAM_BLANK)
IF(LFORCING .OR. LTRANS) WRITE(UNIT=ILUSEG,NML=NAM_FRC)
IF(CCLOUD(1:3) == 'ICE') WRITE(UNIT=ILUSEG,NML=NAM_PARAM_ICE)
IF(CCLOUD == 'C2R2' .OR. CCLOUD == 'C3R5' .OR. CCLOUD == 'KHKO') &
WRITE(UNIT=ILUSEG,NML=NAM_PARAM_C2R2)
IF(CCLOUD == 'C3R5' ) WRITE(UNIT=ILUSEG,NML=NAM_PARAM_C1R3)
IF(CELEC /= 'NONE') WRITE(UNIT=ILUSEG,NML=NAM_ELEC)
IF(LSERIES) WRITE(UNIT=ILUSEG,NML=NAM_SERIES)
IF(NMODEL_CLOUD/=NUNDEF) WRITE(UNIT=ILUSEG,NML=NAM_TURB_CLOUD)
IF(CTURB /= 'NONE') WRITE(UNIT=ILUSEG,NML=NAM_TURB)
!
!
!
!-------------------------------------------------------------------------------
!
!* 2. WRITE UPDATED DESFM ON OUTPUT LISTING
! -------------------------------------
!
IF (NVERB >= 5) THEN
!
ILUOUT = TLUOUT%NLU
!
WRITE(UNIT=ILUOUT,FMT="(/,'DESCRIPTOR OF SEGMENT FOR MODEL ',I2)") KMI
WRITE(UNIT=ILUOUT,FMT="( '------------------------------- ' )")
!
WRITE(UNIT=ILUOUT,FMT="('********** LOGICAL UNITSn **********')")
WRITE(UNIT=ILUOUT,NML=NAM_LUNITn)
!
WRITE(UNIT=ILUOUT,FMT="('********** CONFIGURATIONn **********')")
WRITE(UNIT=ILUOUT,NML=NAM_CONFn)
!
!
WRITE(UNIT=ILUOUT,FMT="('********** DYNAMICn ****************')")
WRITE(UNIT=ILUOUT,NML=NAM_DYNn)
!
WRITE(UNIT=ILUOUT,FMT="('********** ADVECTIONn **************')")
WRITE(UNIT=ILUOUT,NML=NAM_ADVn)
!
WRITE(UNIT=ILUOUT,FMT="('********** PARAMETERIZATIONSn ******')")
WRITE(UNIT=ILUOUT,NML=NAM_PARAMn)
!
WRITE(UNIT=ILUOUT,FMT="('********** RADIATIONn **************')")
WRITE(UNIT=ILUOUT,NML=NAM_PARAM_RADn)
#ifdef MNH_ECRAD
WRITE(UNIT=ILUOUT,FMT="('********** ECRADn **************')")
WRITE(UNIT=ILUOUT,NML=NAM_PARAM_ECRADn)
#endif
450
451
452
453
454
455
456
457
458
459
460
461
462
463
464
465
466
467
468
469
470
471
472
473
474
475
476
477
478
479
480
481
482
483
484
485
486
487
!
WRITE(UNIT=ILUOUT,FMT="('********** CONVECTIONn *************')")
WRITE(UNIT=ILUOUT,NML=NAM_PARAM_KAFRn)
!
WRITE(UNIT=ILUOUT,FMT="('************ PARAM_MFSHALLn *******')")
WRITE(UNIT=ILUOUT,NML=NAM_PARAM_MFSHALLn)
!
WRITE(UNIT=ILUOUT,FMT="('********** LBCn ********************')")
WRITE(UNIT=ILUOUT,NML=NAM_LBCn)
!
WRITE(UNIT=ILUOUT,FMT="('********** NUDGINGn*****************')")
WRITE(UNIT=ILUOUT,NML=NAM_NUDGINGn)
!
WRITE(UNIT=ILUOUT,FMT="('********** TURBn *******************')")
WRITE(UNIT=ILUOUT,NML=NAM_TURBn)
!
WRITE(UNIT=ILUOUT,FMT="('********** CHEMICAL MONITORn *******')")
WRITE(UNIT=ILUOUT,NML=NAM_CH_MNHCn)
!
WRITE(UNIT=ILUOUT,FMT="('************ CHEMICAL SOLVERn ******************')")
WRITE(UNIT=ILUOUT,NML=NAM_CH_SOLVERn)
!
WRITE(UNIT=ILUOUT,FMT="('************ TEMPORAL SERIESn ******************')")
WRITE(UNIT=ILUOUT,NML=NAM_SERIESn)
!
IF (KMI==1) THEN
WRITE(UNIT=ILUOUT,FMT="(/,'PART OF SEGMENT FILE COMMON TO ALL THE MODELS')")
WRITE(UNIT=ILUOUT,FMT="( '---------------------------------------------')")
!
WRITE(UNIT=ILUOUT,FMT="('************ CONFIGURATION ********************')")
WRITE(UNIT=ILUOUT,NML=NAM_CONF)
!
WRITE(UNIT=ILUOUT,FMT="('************ DYNAMIC **************************')")
WRITE(UNIT=ILUOUT,NML=NAM_DYN)
!
WRITE(UNIT=ILUOUT,FMT="(/,'********** NESTING **************************')")
WRITE(UNIT=ILUOUT,NML=NAM_NESTING)
!
! WRITE(UNIT=ILUOUT,FMT="(/,'********** BACKUP ***************************')")
! WRITE(UNIT=ILUOUT,NML=NAM_BACKUP)
!
WAUTELET Philippe
committed
! WRITE(UNIT=ILUOUT,FMT="(/,'********** OUTPUT ***************************')")
! WRITE(UNIT=ILUOUT,NML=NAM_OUTPUT)
!
494
495
496
497
498
499
500
501
502
503
504
505
506
507
508
509
510
511
512
513
514
515
516
517
518
519
520
521
522
523
524
525
526
527
528
529
530
531
532
533
534
535
536
537
538
539
540
541
542
543
544
545
546
547
548
549
550
551
552
553
WRITE(UNIT=ILUOUT,FMT="('************ BUDGET ***************************')")
WRITE(UNIT=ILUOUT,NML=NAM_BUDGET)
!
WRITE(UNIT=ILUOUT,FMT="('************ U BUDGET *************************')")
WRITE(UNIT=ILUOUT,NML=NAM_BU_RU)
!
WRITE(UNIT=ILUOUT,FMT="('************ V BUDGET *************************')")
WRITE(UNIT=ILUOUT,NML=NAM_BU_RV)
!
WRITE(UNIT=ILUOUT,FMT="('************ W BUDGET *************************')")
WRITE(UNIT=ILUOUT,NML=NAM_BU_RW)
!
WRITE(UNIT=ILUOUT,FMT="('************ TH BUDGET ************************')")
WRITE(UNIT=ILUOUT,NML=NAM_BU_RTH)
!
WRITE(UNIT=ILUOUT,FMT="('************ TKE BUDGET ***********************')")
WRITE(UNIT=ILUOUT,NML=NAM_BU_RTKE)
!
WRITE(UNIT=ILUOUT,FMT="('************ RV BUDGET ************************')")
WRITE(UNIT=ILUOUT,NML=NAM_BU_RRV)
!
WRITE(UNIT=ILUOUT,FMT="('************ RC BUDGET ************************')")
WRITE(UNIT=ILUOUT,NML=NAM_BU_RRC)
!
WRITE(UNIT=ILUOUT,FMT="('************ RR BUDGET ************************')")
WRITE(UNIT=ILUOUT,NML=NAM_BU_RRR)
!
WRITE(UNIT=ILUOUT,FMT="('************ RI BUDGET ************************')")
WRITE(UNIT=ILUOUT,NML=NAM_BU_RRI)
!
WRITE(UNIT=ILUOUT,FMT="('************ RS BUDGET ************************')")
WRITE(UNIT=ILUOUT,NML=NAM_BU_RRS)
!
WRITE(UNIT=ILUOUT,FMT="('************ RG BUDGET ************************')")
WRITE(UNIT=ILUOUT,NML=NAM_BU_RRG)
!
WRITE(UNIT=ILUOUT,FMT="('************ RH BUDGET ************************')")
WRITE(UNIT=ILUOUT,NML=NAM_BU_RRH)
!
WRITE(UNIT=ILUOUT,FMT="('************ SVx BUDGET ***********************')")
WRITE(UNIT=ILUOUT,NML=NAM_BU_RSV)
!
WRITE(UNIT=ILUOUT,FMT="('************ LES ******************************')")
WRITE(UNIT=ILUOUT,NML=NAM_LES)
!
WRITE(UNIT=ILUOUT,FMT="('************ BLANK ****************************')")
WRITE(UNIT=ILUOUT,NML=NAM_BLANK)
!
WRITE(UNIT=ILUOUT,FMT="('************ FORCING **************************')")
WRITE(UNIT=ILUOUT,NML=NAM_FRC)
!
WRITE(UNIT=ILUOUT,FMT="('************ ICE SCHEME ***********************')")
WRITE(UNIT=ILUOUT,NML=NAM_PARAM_ICE)
!
WRITE(UNIT=ILUOUT,FMT="('********** DUST SCHEME ************************')")
WRITE(UNIT=ILUOUT,NML=NAM_DUST)
!
WRITE(UNIT=ILUOUT,FMT="('********** PASPOL *****************************')")
WRITE(UNIT=ILUOUT,NML=NAM_PASPOL)
!
#ifdef MNH_FOREFIRE
WRITE(UNIT=ILUOUT,FMT="('********** FOREFIRE *****************************')")
WRITE(UNIT=ILUOUT,NML=NAM_FOREFIRE)
!
#endif
WRITE(UNIT=ILUOUT,FMT="('********** CONDSAMP****************************')")
WRITE(UNIT=ILUOUT,NML=NAM_CONDSAMP)
!
WRITE(UNIT=ILUOUT,FMT="('********** SALT SCHEME ************************')")
WRITE(UNIT=ILUOUT,NML=NAM_SALT)
!
WRITE(UNIT=ILUOUT,FMT="('************ ORILAM SCHEME ********************')")
WRITE(UNIT=ILUOUT,NML=NAM_CH_ORILAM)
!
IF( CCLOUD == 'C2R2' .OR. CCLOUD == 'C3R5') THEN
WRITE(UNIT=ILUOUT,FMT="('*********** C2R2 SCHEME *********************')")
WRITE(UNIT=ILUOUT,NML=NAM_PARAM_C2R2)
IF( CCLOUD == 'C3R5' ) THEN
WRITE(UNIT=ILUOUT,FMT="('*********** C1R3 SCHEME *********************')")
WRITE(UNIT=ILUOUT,NML=NAM_PARAM_C1R3)
END IF
END IF
!
IF( CCLOUD == 'LIMA' ) THEN
WRITE(UNIT=ILUOUT,FMT="('*********** LIMA SCHEME *********************')")
WRITE(UNIT=ILUOUT,NML=NAM_PARAM_LIMA)
END IF
!
IF( CCLOUD == 'KHKO' ) THEN
583
584
585
586
587
588
589
590
591
592
593
594
595
596
597
598
599
600
601
602
603
604
605
606
607
608
609
610
611
612
613
WRITE(UNIT=ILUOUT,FMT="('*********** KHKO SCHEME *********************')")
WRITE(UNIT=ILUOUT,NML=NAM_PARAM_C2R2)
END IF
!
IF( CELEC /= 'NONE' ) THEN
WRITE(UNIT=ILUOUT,FMT="('*********** ELEC SCHEME *********************')")
WRITE(UNIT=ILUOUT,NML=NAM_ELEC)
END IF
!
WRITE(UNIT=ILUOUT,FMT="('************ TEMPORAL SERIES ****************')")
WRITE(UNIT=ILUOUT,NML=NAM_SERIES)
!
WRITE(UNIT=ILUOUT,FMT="('************ MIXING LENGTH FOR CLOUD ***********')")
WRITE(UNIT=ILUOUT,NML=NAM_TURB_CLOUD)
!
END IF
!
END IF
!
IF (CPROGRAM /='MESONH') THEN !return to previous LHORELAX_
LHORELAX_UVWTH = GHORELAX_UVWTH
LHORELAX_RV = GHORELAX_RV
LHORELAX_RC = GHORELAX_RC
LHORELAX_RR = GHORELAX_RR
LHORELAX_RI = GHORELAX_RI
LHORELAX_RS = GHORELAX_RS
LHORELAX_RG = GHORELAX_RG
LHORELAX_TKE = GHORELAX_TKE
LHORELAX_SV(:) = GHORELAX_SV(:)
LHORELAX_SVC2R2= GHORELAX_SVC2R2
LHORELAX_SVC1R3= GHORELAX_SVC1R3
LHORELAX_SVELEC= GHORELAX_SVELEC
LHORELAX_SVCHEM= GHORELAX_SVCHEM
LHORELAX_SVCHIC= GHORELAX_SVCHIC
LHORELAX_SVLG = .FALSE.
LHORELAX_SVDST = GHORELAX_SVDST
LHORELAX_SVSLT = GHORELAX_SVSLT
LHORELAX_SVPP = GHORELAX_SVPP
#ifdef MNH_FOREFIRE
LHORELAX_SVFF = GHORELAX_SVFF
#endif
LHORELAX_SVCS = GHORELAX_SVCS
LHORELAX_SVAER = GHORELAX_SVAER
ELSE
LHORELAX_SV(:) = GHORELAX_SV(:)
ENDIF
CALL UPDATE_NAM_DYNn
!------------------------------------------------------------------------------
!
END SUBROUTINE WRITE_DESFM_n