JOUNRAL DE MERGE VERSION ORCHIDEE ENTRE orchidee-stics_maj et et 1_9_4 ********************************************* NV src_stomate/stomate_turnover.f90 = DONE --------------------------------------------- => vérifier code concernant initialisation de carbres_max => codage en dur de num PFT 11 12 13, probleme a reprendre => modif VUICHARD 20/09/2006 : parametrer la fraction de biomasse aerienne a exporter => écrire tout le vecteur EXPORT ou au moins 11 12 13 ? ********************************************* NV src_stomate/stomate_soilcarbon.f90 = DONE --------------------------------------------- => On aurait besoin de modifier le temps de résidence pour les differents pools de Matiere Organique (Active, Slow, Passive) = PAS FAIT ********************************************* NV src_stomate/stomate_npp.f90 = DONE --------------------------------------------- => codage en dur de num PFT 11 12 13, probleme a reprendre ********************************************* NV src_stomate/stomate_lpj.f90 = DONE ********************************************* NV src_stomate/stomate_litter.f90 = DONE --------------------------------------------- => On aurait besoin de modifier le temps de résidence pour les differents pools de LITIERE (Struct,Metabolic) = PAS FAIT ********************************************* NV src_stomate/stomate_io.f90 = DONE ********************************************* NV src_stomate/stomate_data.f90 = DONE ********************************************* NV src_stomate/stomate_constants.f90 = DONE --------------------------------------------- Sous /home/users/vuichard/ Fichier test1.txt correspond a definition param sans appel a STICS dans /stics_orchidee_maj/... Fichier test2.txt correspond a definition param avec appel a STICS dans /stics_orchidee_maj/... Fichier test3.txt correspond a definition param sans appel a STICS dans /SOURCE_NICOLAS_VIOVY/commit/... VARIABLES pour lesquelles les versions avec STICS et sans different: => lai_max ?? for PFT 1 2 5 => pheno_type_tab => gdd_crit1_tab ?? for PFT 9 184.375 or 270 like in std versions ? => gdd_crit2_tab => gdd_crit3_tab => lowgpp_time_tab ?? for PFT 9 : why 30 instead of 0 ? => hum_min_time_tab => vcmax_opt ?? for PFT 13: why 80 ? 100 more appropriate ? => vjmax_opt => leafagecrit_tab ???? why 100 for PFT 5 when coupling with STICS ? and not 180 => senescence_type_tab => senescence_temp1_tab ??? why 5 for PFT 5 when coupling with STICS ? and not 12 => senescence_hum_tab => min_leaf_age_for_senescence_tab => pheno_model_tab => tphoto_min_c_tab => tphoto_opt_c_tab => tphoto_max_c_tab = Pour les parametres dont les valeurs different entre version pour les PFT naturels, j'ai retenu les valeurs de la version 1.9.4 = Toutes les valeurs de parametres qui avaient ete modifies par Gervois ou Smith pour les PFT crops ont ete introduites dans la 1.9.4 = Deux parametres (leaf_life_tab et leaf_agre_crit_tab) ont ete modifies pour les PFT crops au fil des versions d'orchidee = on a retenu les valeurs de la 1.9.4 (ces valeurs ne correspondent plus aux anciennes valeurs specifiees pour orchidee couple a STICS ********************************************* DL src_stomate/stomate.f90 = TO BE CHECKED --------------------------------------------- => pbme des nums de PFT en dur => stomate_var_xout : t2m_max_daily_fm et swdown_daily_fm sont rajoutées dans l'alloc de la v_commit MAIS cette routine n'existe pas dans la v_194 !!!! NV: je n'ai pas integre cette fonction => dans stomate_main, 5.2, le calcul de t2m_max_daily est rajouté, mais louche - AMAX1( t2m_min(:), t2m_max_daily(:) ) - A VERIFIER NV: verifie => stomate_main 6 : la structure entre la v_commit et la v_194 sont differentes : Il faut savoir si pas de pbme a rajouter l'appel a STICS et la lecture des lai sur ce nouveau pas de temps ne pose pas de probleme avec les op de conversion de surface. NV: pas de problemen a priori = ai rajoute la partie STICS associe Voir la structurepour les 2 versions : v_commit v_194 !! CALL STICS !! (d'ailleurs la limite fait apparaitre un nb en dur if (kjit*dtradia/dt_slow <= 1095)) 6.0 update lai update lai [seulement recalcul des spacenat et space agri [- calcul lai(:,j) = biomass(:,j,ileaf)*sla(j) si defor==1 et qu'on est en fin d'année] pour j=2..nvm et lai(:,0)= zero; ou lecture si forcage] (6.1) [vide] (6.2) [identiques] 6.3 transform GPP from gC/(m**2 of total ground)/day to STOMATE: allocation, phenology, etc. gC/(m**2 of nat/agri ground)/day [call stomate] [call natagritot x 2] 6.4 - transform spatial fractions from fraction of output : transform from dimension nvm to total space to fraction of natural/agricultural space dimension nvm [ . call natagritot(...,ito_natagri,..)x 2, - calculate veget, veget_max, . fait l'update lai comme dans v_194 from veget_cov and veget_cov_max mais * 1/(veget_max(:,j) - lai and height . update classes d'age des feuilles - photosynthesis parameters . #ifdef ctics_crops (mise à jour n_jour++, read lai from stics à njour, attribution au PFTs agris) #endif ?? et en cas de lai force ?? - call stomate [ pareil v_194 ] - transformation inverse des surfaces ********************************************* DL src_stomate/stomate_alloc.f90 = TO BE CHECKED --------------------------------------------- => pbme des nums de PFTs en dur => reste des 0. en initialisation => une modif de nico pour le 'rescaling' des allocations des reserves (en 3.3.3) : Les plantes agri gerees par STICS ne sont pas concernees NICO : A VERIF QUE J'AI BIEN FAIT LA MODIF (notamment le changement de l'indice limite : j < 11 pour v_commit, j < 12 pour v_194) NV: ok mais peut-etre a re-travailler ********************************************* DL src_sechiba/condveg.f90 = DONE --------------------------------------------- ********************************************* DL src_sechiba/diffuco.f90 = TO BE CHECKED --------------------------------------------- => pbme des nums de pfts en dur ********************************************* DL src_sechiba/hydrol.f90 = DONE --------------------------------------------- => Constante 200 renomme comme nblayers TO DO => Creer la contantes nblayers dans constants_soil.f90 ********************************************* DL src_sechiba/hydrolc.f90 = TO BE CHECKED --------------------------------------------- => Constante 200 renommer comme nblayers TO DO => Creer la contantes nblayers dans constants_soil.f90 => dans subroutine hydrols_soil NV+DL => code a regarder plus tart code 6.0 = implementer tel quue dans OS Smith GROS PROBLEME jv == 15 alors que jv ne va pas au-dela de nvp (=14) ??? Implementer tel quel. A REPRENDRE ********************************************* DL src_sechiba/sechiba.f90 = TO BE CHECKED --------------------------------------------- => sehiba_top n'est plus present dans la v_194, pourtant il y a des specificites STICS : (i) declaration swdown_r, sw_down_fm, resp_prec, (ii) rajout de l'attribut SAVE pour ein_jahr, ein_tag (iii) si deja alloue, desalloc swdown_r, swdown_fm; ensuite alloc de ces var (iv) rajout de la taille de swdown dans la variable totsize_1step (v) rajout des ier et des blocks de t2m_max et swdown (vi) swdown_r(:)=swdown_fm(:,iisf) (vii) changement de la liste d'argument de l'appel a slow_proc_main L 2850 (v_comit) swdown_r rajoute, ainsi que t2m_min_r remplace par Nico Vuichard en t2m_r (viii) allocation de resp_prec dans eq soil carbon mensuel, puis init à 0 Il y a aussi a la suite, des modifs de calcul de resp_hetero_tot, bonsoil_input, plus des modif Vuichard 050908 d'apres Viovy A PRIORI PAS A PRENDRE EN COMPTE ?! NICO, SI TU PEUX VERIF STP ... ********************************************* DL src_sechiba/slowproc.f90 = DONE --------------------------------------------- ********************************************* DL src_parameters/constantes_co2.f90 = DONE --------------------------------------------- ********************************************* DL src_parameters/constantes_inca.f90 = DONE --------------------------------------------- => inexistant dans v_194. param concernes, a retrouver ailleurs : splwORCH, em_factor_isoprene, em_factor monoterpene, em_factor_ORVOC, em_factor_OVOC, em_factor_MBO, em_factor_methanol, em_factor_acetone, em_factor_acetal, em_factor_formal, em_factor_acetic, em_factor_formic, em_factor_no_wet, em_factor_no_dry, k_pft, Larch NV+DL => diff ci-dessus pas pris en compte ********************************************* DL src_parameters/constantes_soil.f90 = DONE --------------------------------------------- TO DO => Creer la contantes nblayers dans constants_soil.f90 ********************************************* NV src_parameters/constantes_veg.f90 = DONE --------------------------------------------- COMPILATION rajouter -Dstics_crops dans le Makefile de src_sechiba et src_stomate rajouter interface_stics.f90 comme fichier a compiler dans le Makefile sous src_parameters rajouter "use constantes" dans interface_stics.f90 a la place de "use reqdprec" rajouter "-I/usr/local/install/openmpi-1.3.0/include/" dans Makefile de src_parameters rajouter un elt dans la variable throughfall_by_pft dans hydrol.f90 et hydrolc.f90 rajouter dams Makefile src_parameters: $(MODEL_LIB)(interface_stics.o) $(MODEL_LIB)(interface_stics.o): \