[4] | 1 | !> \file massb-Ice2sea_fixe.f90 |
---|
| 2 | !! Calcule le mass balance en mode perturbation |
---|
| 3 | !! avec la temperature parametree |
---|
| 4 | !< |
---|
| 5 | |
---|
| 6 | !> SUBROUTINE: massb_Ice2sea_fixe |
---|
| 7 | !! \author Cat |
---|
| 8 | !! \date 31 oct 2011 |
---|
| 9 | !! @note Routine qui calcule impose un smb fixe |
---|
| 10 | !! et la temperature en surface avec un lapse rate |
---|
| 11 | !! version groenland voir plus tard si convient Antarctique |
---|
| 12 | !! |
---|
| 13 | !! @note Used modules: |
---|
| 14 | !! @note - use module3D_phy |
---|
| 15 | !! @note - use climat_perturb_mod |
---|
| 16 | !< |
---|
| 17 | |
---|
| 18 | |
---|
| 19 | ! version pour l'antarctique --------------------------------------------------------------- |
---|
| 20 | ! simple copie nettoyee de massb_anteis_perturb |
---|
| 21 | |
---|
| 22 | subroutine massb_Ice2sea_fixe ! calcule le mass balance |
---|
| 23 | |
---|
| 24 | use module3D_phy |
---|
| 25 | use climat_Grice2sea_mod |
---|
| 26 | |
---|
| 27 | implicit none |
---|
| 28 | |
---|
| 29 | ! surface temperature et accumulation |
---|
| 30 | |
---|
| 31 | Tann (:,:) = Ta0 (:,:) + T_lapse_rate * (S(:,:)-S0(:,:)) |
---|
| 32 | Ts(:,:) = Tann(:,:) |
---|
| 33 | |
---|
| 34 | ! bm est inchange |
---|
| 35 | |
---|
| 36 | |
---|
| 37 | ! ablation (fonction de T et acc) est maintenant appelee dans le main |
---|
| 38 | ! et dans steps_time_loop. |
---|
| 39 | ! On fait donc une routine fake ici |
---|
| 40 | |
---|
| 41 | ! massb_Ice2sea_fixe est appeles par forclim du module climat_Grice2sea_mod |
---|
| 42 | |
---|
| 43 | |
---|
| 44 | debug_3D(:,:,29)=Tann(:,:)-Ta0(:,:) |
---|
| 45 | debug_3D(:,:,30)=Acc(:,:)-precip(:,:) |
---|
| 46 | !debug_3D(:,:,31)=Acc(:,:)-bmelt(:,:) |
---|
| 47 | |
---|
| 48 | return |
---|
| 49 | end subroutine massb_Ice2sea_fixe |
---|