Changes between Version 7 and Version 8 of DevelopmentActivities/MergeHydro/Martial_notes_on_merge
- Timestamp:
- 2011-10-20T06:40:50+02:00 (13 years ago)
Legend:
- Unmodified
- Added
- Removed
- Modified
-
DevelopmentActivities/MergeHydro/Martial_notes_on_merge
v7 v8 36 36 }}} 37 37 Il me semble utile de rajouter un getin pour cela, non ? 38 38 * Concernant les erreurs de fermeture de bilan (cf variable waterbal_error), on indique plusieurs fois dans le code que l'on va l'arrêter : 39 {{{ 40 waterbal_error=.TRUE. 41 CALL ipslerr(2, 'hydrol_split_soil', 'We will STOP after hydrol_split_soil.',& 42 & 'check_CWRR','PRECISOL SPLIT FALSE') 43 }}} 44 sans le faire ! Est-ce à rajouter ? 39 45 == problèmes conversion (frac_bare/veget) versus (veget,veget_max) == 40 41 46 * dans hydrol_init : on passait veget (ie l'ancien veget_max) et on a pas changer le code par rapport 47 à l'ancienne version. On obtient peut-être des incohérences : 42 48 {{{ 43 49 IF ( MINVAL(resdist) .EQ. MAXVAL(resdist) .AND. MINVAL(resdist) .EQ. val_exp) THEN … … 51 57 ENDDO 52 58 }}} 53 54 59 Dois-je laisser le veget de la version standard ?? a priori, je passe à veget_max car vegtot fait référence à 1-frac_nobio et est souvent divisé par veget_max. 60 * dans hydrol_vegetupd on utilise vraiment frac_bare : 55 61 {{{ 56 62 DO jv = 1, nvm … … 60 66 END DO 61 67 }}} 62 68 On doit donc le recalculer à partir de la formule dans slowproc : 63 69 {{{ 64 70 frac_bare(:,:) = zero … … 70 76 ENDIF 71 77 }}} 72 78 et la définition actuelle du veget : 73 79 {{{ 74 80 ! Ajout Nouveau calcul (stomate-like) … … 82 88 [...] 83 89 }}} 84 90 On déduit alors en combinant ces deux définitions : 85 91 {{{ 86 92 frac_bare(:,1) = un … … 93 99 ENDDO 94 100 }}} 95 96 97 101 L'utilisation du undef_sechiba me paraît cohérente lorsque l'on a pas de présence du PFT sur le pixel.[[BR]] 102 On avait d'ailleurs une erreur avant avec le frac_bare car il était parfois définit (utilisé ?) lorsque veget_max(ji,jv) était nul. [[BR]] 103 Un problème surgit alors lorsque l'on définit : 98 104 {{{ 99 105 tot_frac_bare(:) = zero … … 105 111 ENDDO 106 112 }}} 107 113 car on ne teste pas veget_max > min_sechiba. Je le rajoute. 108 114 * hydrol_canop : je me pose la question du veget/veget_max ici car on définit qsintveg et precisol dans cette routine. 109 115 Si l'on passe veget_max, on ne tient plus compte du LAI pour calculer ces deux variables importantes. Je laisse donc veget.