From d3ef9cbeac31818b167345e5a73d0fda4533900d Mon Sep 17 00:00:00 2001
From: Juan Escobar <escj@aero.obs-mip.fr>
Date: Mon, 3 Jul 2017 15:45:13 +0200
Subject: [PATCH] Maud & Juan 1/07/2017: Add DIAG chimical surface fluxes

---
 src/MNH/modd_ch_flxn.f90 | 66 ++++++++++++++++++++++++++++++++++++++++
 1 file changed, 66 insertions(+)
 create mode 100644 src/MNH/modd_ch_flxn.f90

diff --git a/src/MNH/modd_ch_flxn.f90 b/src/MNH/modd_ch_flxn.f90
new file mode 100644
index 000000000..47c729a93
--- /dev/null
+++ b/src/MNH/modd_ch_flxn.f90
@@ -0,0 +1,66 @@
+!MNH_LIC Copyright 1994-2014 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.
+!-----------------------------------------------------------------
+!--------------- special set of characters for RCS information
+!-----------------------------------------------------------------
+! $Source: /home/cvsroot/MNH-VX-Y-Z/src/MNH/modd_ch_flxn.f90,v $ $Revision: 1.1 $
+! MASDEV5_2 modd 2016/06/27 14:05:40
+!-----------------------------------------------------------------
+!     #####################
+      MODULE MODD_CH_FLX_n
+!     ######################
+!
+!!
+!!    PURPOSE
+!!    -------
+!     Save the net surface flux at the surface 
+!     for output with diag
+!
+!!
+!!**  IMPLICIT ARGUMENTS
+!!    ------------------
+!!      None
+!!
+!
+!!    AUTHOR
+!!    ------
+!!  P. Tulet   *Meteo France*
+!!
+!!    MODIFICATIONS
+!!    -------------
+!!  12/07/16 (M. Leriche) keep only the flux
+!!  
+!------------------------------------------------------------------------------
+!
+!*       0.   DECLARATIONS
+!             ------------
+!
+USE MODD_PARAMETERS, ONLY: JPMODELMAX
+IMPLICIT NONE
+
+TYPE CH_FLX_t
+!
+  REAL, DIMENSION(:,:,:), POINTER :: XCHFLX=>NULL() ! chemical fluxes ppp.m/s at t
+!
+END TYPE CH_FLX_t
+
+TYPE(CH_FLX_t), DIMENSION(JPMODELMAX), TARGET, SAVE :: CH_FLX_MODEL
+
+REAL, DIMENSION(:,:,:), POINTER :: XCHFLX=>NULL()
+
+CONTAINS
+
+SUBROUTINE CH_FLX_GOTO_MODEL(KFROM, KTO)
+INTEGER, INTENT(IN) :: KFROM, KTO
+!
+! Save current state for allocated arrays
+CH_FLX_MODEL(KFROM)%XCHFLX=>XCHFLX
+!
+! Current model is set to model KTO
+XCHFLX=>CH_FLX_MODEL(KTO)%XCHFLX
+
+END SUBROUTINE CH_FLX_GOTO_MODEL
+
+END MODULE MODD_CH_FLX_n
-- 
GitLab