Main changes in Surfex v7_3
- ISBA coupling with the TOPMODEL library (B. Vincendon)
- Chemical part: NO fluxes (C. Delon)
- if LCH_NO_FLUX (NAM_CH_ISBAn), NO flux is caculated in ISBA part
- if LCH_NO_FLUX (NAM_CH_ISBAn), NO flux is written in NOFLUX output
- with LCH_NO_FLUX, PH and FERT fields are given at PGD step (NAM_ISBA)
- Chemical part: SNAP emissions (S. Queguiner)
- to use other types of cadastres for the chemical emissions (smaller input files)
- Contribution from GMGEC (B. Decharme)
- New initialization for ALBEDO (D. Carrer & C. Meurey)
- Like for the LAI, soil and vegetation ALBEDO are derived from satellite products (SAF) and defined for each cover and each vegetation type in ECOCLIMAP.
- Contribution from CEN (M. Lafaysse)
- possibility during PREP to initialize TG, WG and WGI with NETCDF files in ISBA part.
- modification of the management of CROCUS optimal thicknesses for very thick coverings of snow.
- partial merge of routines SNOWCRORAD and SNOWCROALB to supress the code that was duplicated
- possibility to simulate glaciers with crocus (albedo forced in the 3 bands and different roughness for the apparent ice)
- all constants concerning albedo routines are put in modd_snow_par to be declared only once.
- new constants for ice albedo and roughness can be defined in namelist NAM_SURF_CSTS.
Default values : ! Roughness length ratio between ice and snow XZ0ICEZ0SNOW=10. ! 3 bands spectral albedo for glacier ice (CROCUS) ! Default values from Lejeune et al 2009 (Zongo, Bolivia) XALBICE1=0.38 XALBICE2=0.23 XALBICE3=0.08 ! Density threshold for ice detection kg.m-3 XRHOTHRESHOLD_ICE=850.
- variable PSNOWAGE is now an age and no longer a date (what allows to do easier spinup when CROCUS is activated)
- print option for CROCUS debugging, conditioned by the definition of environment variables
- Modifications for flake fluxes and 1D ocean model (P. Le Moigne & P. Peyrillé)
- possibility to activate ECUME in Flake (CFLK_FLUX=’ECUME’ in NAM_FLAKEn) + associated options LPRECIP and LPWEBB
- add of relaxation for the ocean model: temperature, salinity, currents
- possibility to correct fluxes with a coefficient given in namelist (XCORFLX)
- possibility to set surface fluxes to 0
- flag diapycnal : calculation of mixing coefficients at internal waves and wind shearing
- Contribution for the assimilation part (SODA) (Trygve Aspelien)
- changes for the snow assimilation
- New TEB (V. Masson, A. Lemonsu, G. Pigeon, C. De Munck)
- Parallelization of the OFFLINE driver (F. Suzat & S. Faroux)
- parallelization possible with MPI and / or OPEN-MP
- possibility to unactivate the MPI implementation (export VER_MPI=NOMPI during configure part of installation)
- For Méso-NH (M. Tomasini):
- Add an index in the second dimension of the ISBA array variables for BILIN interpolation routine to not bug in case 2D (this is not the more beautiful method; the BILIN routine should better be adapted)
- Side developments:
- all fields names have now 12 characters and no longer 16.
=> In NETCDF forcing files, FORC_TIME_STEP becomes FRC_TIME_STP. To change it easily, do:
ncrename -v FORC_TIME_STEP,FRC_TIME_STP FORCING.nc
- all fields names have now 12 characters and no longer 16.