Skip to content
Snippets Groups Projects
Commit 84e045a4 authored by Gaelle Tanguy's avatar Gaelle Tanguy Committed by WAUTELET Philippe
Browse files

S.Bielli 24/11/2014 : bug when writing netcdf file with LTURB_DIAG=T and LTURB_FLX=F or conversely

parent 8a66ffee
No related branches found
No related tags found
No related merge requests found
...@@ -945,38 +945,6 @@ IF (CTURBDIM=='3DIM') THEN ...@@ -945,38 +945,6 @@ IF (CTURBDIM=='3DIM') THEN
END IF END IF
END IF END IF
!! !!
IF ( LTURB_DIAG ) THEN
! stores the mixing length
!
YRECFM ='LM'
YCOMMENT='X_Y_Z_LM (M)'
IGRID = 1
ILENCH=LEN(YCOMMENT)
CALL FMWRIT(HFMFILE,YRECFM,CLUOUT,'XY',ZSTORE_3D,IGRID,ILENCH,YCOMMENT,IRESP)
!!
IF (NRR /= 0) THEN
!
! stores the conservative potential temperature
!
YRECFM ='THLM'
YCOMMENT='X_Y_Z_THLM (KELVIN)'
IGRID = 1
ILENCH=LEN(YCOMMENT)
CALL FMWRIT(HFMFILE,YRECFM,CLUOUT,'XY',ZSTORE_3D,IGRID,ILENCH,YCOMMENT,IRESP)
!!
! stores the conservative mixing ratio
!
YRECFM ='RNPM'
YCOMMENT='X_Y_Z_RNPM (KG/KG)'
IGRID = 1
ILENCH=LEN(YCOMMENT)
CALL FMWRIT(HFMFILE,YRECFM,CLUOUT,'XY',ZSTORE_3D,IGRID,ILENCH, &
YCOMMENT,IRESP)
END IF
!!
!!
END IF
!!!!!! turb_hor_uw
IF ( LTURB_FLX ) THEN IF ( LTURB_FLX ) THEN
! stores <U W> ! stores <U W>
YRECFM ='UW_HFLX' YRECFM ='UW_HFLX'
...@@ -1014,6 +982,36 @@ IF ( LTURB_FLX ) THEN ...@@ -1014,6 +982,36 @@ IF ( LTURB_FLX ) THEN
END DO END DO
END IF END IF
!! !!
IF ( LTURB_DIAG ) THEN
! stores the dynamic production
YRECFM ='DP'
YCOMMENT='X_Y_Z_DP (M**2/S**3)'
IGRID = 1
ILENCH=LEN(YCOMMENT)
CALL FMWRIT(HFMFILE,YRECFM,CLUOUT,'XY',ZSTORE_3D,IGRID,ILENCH,YCOMMENT,IRESP)
!!
! stores the thermal production
YRECFM ='TP'
YCOMMENT='X_Y_Z_TP (M**2/S**3)'
IGRID = 1
ILENCH=LEN(YCOMMENT)
CALL FMWRIT(HFMFILE,YRECFM,CLUOUT,'XY',ZSTORE_3D,IGRID,ILENCH,YCOMMENT,IRESP)
!!
! stores the whole turbulent transport
YRECFM ='TR'
YCOMMENT='X_Y_Z_TR (M**2/S**3)'
IGRID = 1
ILENCH=LEN(YCOMMENT)
CALL FMWRIT(HFMFILE,YRECFM,CLUOUT,'XY',ZSTORE_3D,IGRID,ILENCH,YCOMMENT,IRESP)
!!
! stores the dissipation of TKE
YRECFM ='DISS'
YCOMMENT='X_Y_Z_DISS (M**2/S**3)'
IGRID = 1
ILENCH=LEN(YCOMMENT)
CALL FMWRIT(HFMFILE,YRECFM,CLUOUT,'XY',ZSTORE_3D,IGRID,ILENCH,YCOMMENT,IRESP)
END IF
!
IF ( LTURB_DIAG ) THEN IF ( LTURB_DIAG ) THEN
! stores the mixing length ! stores the mixing length
YRECFM ='LM' YRECFM ='LM'
...@@ -1061,35 +1059,35 @@ IF ( LTURB_DIAG ) THEN ...@@ -1061,35 +1059,35 @@ IF ( LTURB_DIAG ) THEN
!! !!
END IF END IF
!!!! tke_eps_sources !!!! tke_eps_sources
IF ( LTURB_DIAG ) THEN !IF ( LTURB_DIAG ) THEN
! stores the dynamic production ! ! stores the dynamic production
YRECFM ='DP' ! YRECFM ='DP'
YCOMMENT='X_Y_Z_DP (M**2/S**3)' ! YCOMMENT='X_Y_Z_DP (M**2/S**3)'
IGRID = 1 ! IGRID = 1
ILENCH=LEN(YCOMMENT) ! ILENCH=LEN(YCOMMENT)
CALL FMWRIT(HFMFILE,YRECFM,CLUOUT,'XY',ZSTORE_3D,IGRID,ILENCH,YCOMMENT,IRESP) ! CALL FMWRIT(HFMFILE,YRECFM,CLUOUT,'XY',ZSTORE_3D,IGRID,ILENCH,YCOMMENT,IRESP)
!! !!!
! stores the thermal production ! ! stores the thermal production
YRECFM ='TP' ! YRECFM ='TP'
YCOMMENT='X_Y_Z_TP (M**2/S**3)' ! YCOMMENT='X_Y_Z_TP (M**2/S**3)'
IGRID = 1 ! IGRID = 1
ILENCH=LEN(YCOMMENT) ! ILENCH=LEN(YCOMMENT)
CALL FMWRIT(HFMFILE,YRECFM,CLUOUT,'XY',ZSTORE_3D,IGRID,ILENCH,YCOMMENT,IRESP) ! CALL FMWRIT(HFMFILE,YRECFM,CLUOUT,'XY',ZSTORE_3D,IGRID,ILENCH,YCOMMENT,IRESP)
!! !!!
! stores the whole turbulent transport ! ! stores the whole turbulent transport
YRECFM ='TR' ! YRECFM ='TR'
YCOMMENT='X_Y_Z_TR (M**2/S**3)' ! YCOMMENT='X_Y_Z_TR (M**2/S**3)'
IGRID = 1 ! IGRID = 1
ILENCH=LEN(YCOMMENT) ! ILENCH=LEN(YCOMMENT)
CALL FMWRIT(HFMFILE,YRECFM,CLUOUT,'XY',ZSTORE_3D,IGRID,ILENCH,YCOMMENT,IRESP) ! CALL FMWRIT(HFMFILE,YRECFM,CLUOUT,'XY',ZSTORE_3D,IGRID,ILENCH,YCOMMENT,IRESP)
!! !!!
! stores the dissipation of TKE ! ! stores the dissipation of TKE
YRECFM ='DISS' ! YRECFM ='DISS'
YCOMMENT='X_Y_Z_DISS (M**2/S**3)' ! YCOMMENT='X_Y_Z_DISS (M**2/S**3)'
IGRID = 1 ! IGRID = 1
ILENCH=LEN(YCOMMENT) ! ILENCH=LEN(YCOMMENT)
CALL FMWRIT(HFMFILE,YRECFM,CLUOUT,'XY',ZSTORE_3D,IGRID,ILENCH,YCOMMENT,IRESP) ! CALL FMWRIT(HFMFILE,YRECFM,CLUOUT,'XY',ZSTORE_3D,IGRID,ILENCH,YCOMMENT,IRESP)
END IF !END IF
!! !!
!!!!!!!!! Shallow_mf_pack !!!!!!!!! Shallow_mf_pack
IF (CSCONV == 'EDKF') THEN IF (CSCONV == 'EDKF') THEN
......
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment