| 1 | JOUNRAL DE MERGE VERSION ORCHIDEE ENTRE orchidee-stics_maj et et 1_9_4 |
|---|
| 2 | |
|---|
| 3 | ********************************************* |
|---|
| 4 | NV src_stomate/stomate_turnover.f90 = DONE |
|---|
| 5 | --------------------------------------------- |
|---|
| 6 | => vérifier code concernant initialisation de carbres_max |
|---|
| 7 | => codage en dur de num PFT 11 12 13, probleme a reprendre |
|---|
| 8 | => modif VUICHARD 20/09/2006 : parametrer la fraction de biomasse aerienne a exporter |
|---|
| 9 | => écrire tout le vecteur EXPORT ou au moins 11 12 13 ? |
|---|
| 10 | ********************************************* |
|---|
| 11 | NV src_stomate/stomate_soilcarbon.f90 = DONE |
|---|
| 12 | --------------------------------------------- |
|---|
| 13 | => On aurait besoin de modifier le temps de résidence pour les differents pools de Matiere Organique (Active, Slow, Passive) = PAS FAIT |
|---|
| 14 | ********************************************* |
|---|
| 15 | NV src_stomate/stomate_npp.f90 = DONE |
|---|
| 16 | --------------------------------------------- |
|---|
| 17 | => codage en dur de num PFT 11 12 13, probleme a reprendre |
|---|
| 18 | ********************************************* |
|---|
| 19 | NV src_stomate/stomate_lpj.f90 = DONE |
|---|
| 20 | ********************************************* |
|---|
| 21 | NV src_stomate/stomate_litter.f90 = DONE |
|---|
| 22 | --------------------------------------------- |
|---|
| 23 | => On aurait besoin de modifier le temps de résidence pour les differents pools de LITIERE (Struct,Metabolic) = PAS FAIT |
|---|
| 24 | ********************************************* |
|---|
| 25 | NV src_stomate/stomate_io.f90 = DONE |
|---|
| 26 | ********************************************* |
|---|
| 27 | NV src_stomate/stomate_data.f90 = DONE |
|---|
| 28 | ********************************************* |
|---|
| 29 | NV src_stomate/stomate_constants.f90 = DONE |
|---|
| 30 | --------------------------------------------- |
|---|
| 31 | Sous /home/users/vuichard/ |
|---|
| 32 | Fichier test1.txt correspond a definition param sans appel a STICS dans /stics_orchidee_maj/... |
|---|
| 33 | Fichier test2.txt correspond a definition param avec appel a STICS dans /stics_orchidee_maj/... |
|---|
| 34 | Fichier test3.txt correspond a definition param sans appel a STICS dans /SOURCE_NICOLAS_VIOVY/commit/... |
|---|
| 35 | VARIABLES pour lesquelles les versions avec STICS et sans different: |
|---|
| 36 | => lai_max ?? for PFT 1 2 5 |
|---|
| 37 | => pheno_type_tab |
|---|
| 38 | => gdd_crit1_tab ?? for PFT 9 184.375 or 270 like in std versions ? |
|---|
| 39 | => gdd_crit2_tab |
|---|
| 40 | => gdd_crit3_tab |
|---|
| 41 | => lowgpp_time_tab ?? for PFT 9 : why 30 instead of 0 ? |
|---|
| 42 | => hum_min_time_tab |
|---|
| 43 | => vcmax_opt ?? for PFT 13: why 80 ? 100 more appropriate ? |
|---|
| 44 | => vjmax_opt |
|---|
| 45 | => leafagecrit_tab ???? why 100 for PFT 5 when coupling with STICS ? and not 180 |
|---|
| 46 | => senescence_type_tab |
|---|
| 47 | => senescence_temp1_tab ??? why 5 for PFT 5 when coupling with STICS ? and not 12 |
|---|
| 48 | => senescence_hum_tab |
|---|
| 49 | => min_leaf_age_for_senescence_tab |
|---|
| 50 | => pheno_model_tab |
|---|
| 51 | => tphoto_min_c_tab |
|---|
| 52 | => tphoto_opt_c_tab |
|---|
| 53 | => tphoto_max_c_tab |
|---|
| 54 | = Pour les parametres dont les valeurs different entre version pour les PFT naturels, j'ai retenu les valeurs de la version 1.9.4 |
|---|
| 55 | = 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 |
|---|
| 56 | = 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 |
|---|
| 57 | |
|---|
| 58 | ********************************************* |
|---|
| 59 | DL src_stomate/stomate.f90 = TO BE CHECKED |
|---|
| 60 | --------------------------------------------- |
|---|
| 61 | => pbme des nums de PFT en dur |
|---|
| 62 | => stomate_var_xout : t2m_max_daily_fm et swdown_daily_fm |
|---|
| 63 | sont rajoutées dans l'alloc de la v_commit |
|---|
| 64 | MAIS cette routine n'existe pas dans la v_194 !!!! |
|---|
| 65 | NV: je n'ai pas integre cette fonction |
|---|
| 66 | |
|---|
| 67 | => dans stomate_main, 5.2, le calcul de t2m_max_daily est rajouté, |
|---|
| 68 | mais louche - AMAX1( t2m_min(:), t2m_max_daily(:) ) - |
|---|
| 69 | A VERIFIER |
|---|
| 70 | NV: verifie |
|---|
| 71 | |
|---|
| 72 | => stomate_main 6 : la structure entre la v_commit et la v_194 sont differentes : |
|---|
| 73 | Il faut savoir si pas de pbme a rajouter l'appel a STICS et la lecture des lai |
|---|
| 74 | sur ce nouveau pas de temps ne pose pas de probleme avec les op de conversion de surface. |
|---|
| 75 | NV: pas de problemen a priori = ai rajoute la partie STICS associe |
|---|
| 76 | |
|---|
| 77 | Voir la structurepour les 2 versions : |
|---|
| 78 | |
|---|
| 79 | v_commit v_194 |
|---|
| 80 | |
|---|
| 81 | !! CALL STICS !! |
|---|
| 82 | (d'ailleurs la limite fait apparaitre un nb en dur |
|---|
| 83 | if (kjit*dtradia/dt_slow <= 1095)) |
|---|
| 84 | |
|---|
| 85 | 6.0 update lai update lai |
|---|
| 86 | [seulement recalcul des spacenat et space agri [- calcul lai(:,j) = biomass(:,j,ileaf)*sla(j) |
|---|
| 87 | si defor==1 et qu'on est en fin d'année] pour j=2..nvm et lai(:,0)= zero; ou lecture si |
|---|
| 88 | forcage] |
|---|
| 89 | |
|---|
| 90 | (6.1) [vide] (6.2) [identiques] |
|---|
| 91 | |
|---|
| 92 | |
|---|
| 93 | 6.3 transform GPP from gC/(m**2 of total ground)/day to STOMATE: allocation, phenology, etc. |
|---|
| 94 | gC/(m**2 of nat/agri ground)/day [call stomate] |
|---|
| 95 | [call natagritot x 2] |
|---|
| 96 | |
|---|
| 97 | 6.4 - transform spatial fractions from fraction of output : transform from dimension nvm to |
|---|
| 98 | total space to fraction of natural/agricultural space dimension nvm |
|---|
| 99 | [ . call natagritot(...,ito_natagri,..)x 2, - calculate veget, veget_max, |
|---|
| 100 | . fait l'update lai comme dans v_194 from veget_cov and veget_cov_max |
|---|
| 101 | mais * 1/(veget_max(:,j) - lai and height |
|---|
| 102 | . update classes d'age des feuilles - photosynthesis parameters |
|---|
| 103 | . #ifdef ctics_crops |
|---|
| 104 | (mise à jour n_jour++, read lai from stics à njour, |
|---|
| 105 | attribution au PFTs agris) |
|---|
| 106 | #endif |
|---|
| 107 | ?? et en cas de lai force ?? |
|---|
| 108 | - call stomate [ pareil v_194 ] |
|---|
| 109 | - transformation inverse des surfaces |
|---|
| 110 | |
|---|
| 111 | ********************************************* |
|---|
| 112 | DL src_stomate/stomate_alloc.f90 = TO BE CHECKED |
|---|
| 113 | --------------------------------------------- |
|---|
| 114 | => pbme des nums de PFTs en dur |
|---|
| 115 | => reste des 0. en initialisation |
|---|
| 116 | => une modif de nico pour le 'rescaling' des allocations des reserves (en 3.3.3) : |
|---|
| 117 | Les plantes agri gerees par STICS ne sont pas concernees |
|---|
| 118 | NICO : A VERIF QUE J'AI BIEN FAIT LA MODIF |
|---|
| 119 | (notamment le changement de l'indice limite : j < 11 pour v_commit, j < 12 pour v_194) |
|---|
| 120 | NV: ok mais peut-etre a re-travailler |
|---|
| 121 | |
|---|
| 122 | ********************************************* |
|---|
| 123 | DL src_sechiba/condveg.f90 = DONE |
|---|
| 124 | --------------------------------------------- |
|---|
| 125 | |
|---|
| 126 | ********************************************* |
|---|
| 127 | DL src_sechiba/diffuco.f90 = TO BE CHECKED |
|---|
| 128 | --------------------------------------------- |
|---|
| 129 | => pbme des nums de pfts en dur |
|---|
| 130 | |
|---|
| 131 | ********************************************* |
|---|
| 132 | DL src_sechiba/hydrol.f90 = DONE |
|---|
| 133 | --------------------------------------------- |
|---|
| 134 | => Constante 200 renomme comme nblayers |
|---|
| 135 | TO DO => Creer la contantes nblayers dans constants_soil.f90 |
|---|
| 136 | |
|---|
| 137 | |
|---|
| 138 | ********************************************* |
|---|
| 139 | DL src_sechiba/hydrolc.f90 = TO BE CHECKED |
|---|
| 140 | --------------------------------------------- |
|---|
| 141 | => Constante 200 renommer comme nblayers |
|---|
| 142 | TO DO => Creer la contantes nblayers dans constants_soil.f90 |
|---|
| 143 | => dans subroutine hydrols_soil |
|---|
| 144 | 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 |
|---|
| 145 | |
|---|
| 146 | ********************************************* |
|---|
| 147 | DL src_sechiba/sechiba.f90 = TO BE CHECKED |
|---|
| 148 | --------------------------------------------- |
|---|
| 149 | |
|---|
| 150 | |
|---|
| 151 | => sehiba_top n'est plus present dans la v_194, pourtant il y a des |
|---|
| 152 | specificites STICS : |
|---|
| 153 | (i) declaration swdown_r, sw_down_fm, resp_prec, |
|---|
| 154 | (ii) rajout de l'attribut SAVE pour ein_jahr, ein_tag |
|---|
| 155 | (iii) si deja alloue, desalloc swdown_r, swdown_fm; ensuite alloc de ces var |
|---|
| 156 | (iv) rajout de la taille de swdown dans la variable totsize_1step |
|---|
| 157 | (v) rajout des ier et des blocks de t2m_max et swdown |
|---|
| 158 | (vi) swdown_r(:)=swdown_fm(:,iisf) |
|---|
| 159 | (vii) changement de la liste d'argument de l'appel a slow_proc_main L 2850 (v_comit) |
|---|
| 160 | swdown_r rajoute, ainsi que t2m_min_r remplace par Nico Vuichard en t2m_r |
|---|
| 161 | (viii) allocation de resp_prec dans eq soil carbon mensuel, puis init à 0 |
|---|
| 162 | |
|---|
| 163 | Il y a aussi a la suite, des modifs de calcul de resp_hetero_tot, bonsoil_input, |
|---|
| 164 | plus des modif Vuichard 050908 d'apres Viovy |
|---|
| 165 | A PRIORI PAS A PRENDRE EN COMPTE ?! |
|---|
| 166 | |
|---|
| 167 | NICO, SI TU PEUX VERIF STP ... |
|---|
| 168 | |
|---|
| 169 | ********************************************* |
|---|
| 170 | DL src_sechiba/slowproc.f90 = DONE |
|---|
| 171 | --------------------------------------------- |
|---|
| 172 | |
|---|
| 173 | ********************************************* |
|---|
| 174 | DL src_parameters/constantes_co2.f90 = DONE |
|---|
| 175 | --------------------------------------------- |
|---|
| 176 | |
|---|
| 177 | ********************************************* |
|---|
| 178 | DL src_parameters/constantes_inca.f90 = DONE |
|---|
| 179 | --------------------------------------------- |
|---|
| 180 | => inexistant dans v_194. |
|---|
| 181 | param concernes, a retrouver ailleurs : |
|---|
| 182 | splwORCH, em_factor_isoprene, em_factor monoterpene, |
|---|
| 183 | em_factor_ORVOC, em_factor_OVOC, em_factor_MBO, |
|---|
| 184 | em_factor_methanol, em_factor_acetone, em_factor_acetal, |
|---|
| 185 | em_factor_formal, em_factor_acetic, em_factor_formic, |
|---|
| 186 | em_factor_no_wet, em_factor_no_dry, k_pft, Larch |
|---|
| 187 | NV+DL => diff ci-dessus pas pris en compte |
|---|
| 188 | |
|---|
| 189 | ********************************************* |
|---|
| 190 | DL src_parameters/constantes_soil.f90 = DONE |
|---|
| 191 | --------------------------------------------- |
|---|
| 192 | TO DO => Creer la contantes nblayers dans constants_soil.f90 |
|---|
| 193 | |
|---|
| 194 | ********************************************* |
|---|
| 195 | NV src_parameters/constantes_veg.f90 = DONE |
|---|
| 196 | --------------------------------------------- |
|---|
| 197 | |
|---|
| 198 | |
|---|
| 199 | COMPILATION |
|---|
| 200 | |
|---|
| 201 | rajouter -Dstics_crops dans le Makefile de src_sechiba et src_stomate |
|---|
| 202 | rajouter interface_stics.f90 comme fichier a compiler dans le Makefile sous src_parameters |
|---|
| 203 | rajouter "use constantes" dans interface_stics.f90 a la place de "use reqdprec" |
|---|
| 204 | rajouter "-I/usr/local/install/openmpi-1.3.0/include/" dans Makefile de src_parameters |
|---|
| 205 | rajouter un elt dans la variable throughfall_by_pft dans hydrol.f90 et hydrolc.f90 |
|---|
| 206 | rajouter dams Makefile src_parameters: |
|---|
| 207 | $(MODEL_LIB)(interface_stics.o) |
|---|
| 208 | $(MODEL_LIB)(interface_stics.o): \ |
|---|
