Changeset 24


Ignore:
Timestamp:
12/08/15 15:22:18 (8 years ago)
Author:
dumas
Message:

Correction Bug dans le calcul de la temperature : reecriture sans les modules-procedures

Location:
trunk/SOURCES
Files:
6 added
6 deleted
4 edited

Legend:

Unmodified
Added
Removed

  • trunk/SOURCES/Makefile.grisli.inc

    r22 r24  
    6464        massb_perturb_mois.o massb_perturb_mois_lapsecouche.o  
    6565 
    66 mod_communs =   util_recovery.o \ 
    67         geom_type.o temperature_type.o ice_flow_type.o mask_flgz_type.o \ 
    68         interface_input.o deformation_type.o autre_pr_temp_type.o  \ 
     66mod_communs = util_recovery.o \ 
     67        interface_input.o \ 
    6968        iso_declar_mod-0.3.o isostasie_mod-0.3.o noisostasie_mod-0.3.o \ 
    70         deformation_mod_2lois.o deformation_mod_2lois_isotherme.o     \ 
    71         sliding_Bindshadler_mod.o no_sliding.o                      \ 
     69        deformation_mod_2lois.o deformation_mod_2lois_isotherme.o     \ 
     70        sliding_Bindshadler_mod.o no_sliding.o                      \ 
    7271        dragging_hwatermax_0.2_mod.o dragging_calc_beta_mod.o  \ 
    7372        dragging_hwat-contigu_mod.o dragging_hwat_contmaj_mod.o \ 
    … …  
    8281        tridiag_mod-0.3.o out_horiz_mod.o bmelt-seuil-profondeur_mod.o  
    8382 
     83        #    geom_type.o temperature_type.o ice_flow_type.o mask_flgz_type.o \ 
     84        #deformation_type.o autre_pr_temp_type.o  \ 
     85         
    8486mod_communs_RGL =   util_recovery.o \ 
    85         geom_type.o temperature_type.o ice_flow_type.o mask_flgz_type.o \ 
    86         interface_input.o deformation_type.o autre_pr_temp_type.o  \ 
     87        interface_input.o \ 
    8788        iso_declar_mod-0.3.o isostasie_mod-0.3.o noisostasie_mod-0.3.o \ 
    8889        deformation_mod_2lois.o deformation_mod_2lois_isotherme.o      \ 
    … …  
    100101        tridiag_mod-0.3.o out_horiz_mod.o bmelt-seuil-profondeur_mod.o 
    101102 
    102  
    103  
    104 # Modules communs Post geographiques  
    105  
    106 mod_post_geo = step_temp_type.o 
     103#        geom_type.o temperature_type.o ice_flow_type.o mask_flgz_type.o deformation_type.o autre_pr_temp_type.o \ 
    107104 
    108105 
    … …  
    133130# enleve  : flow_general-0.3.o flowlaw-0.3.o  Cat 23 mars 2012 
    134131 
    135 routines_communes = tracebug.o ablation_bord.o allocate_types.o \ 
    136         icetemp_declar-0.5_mod.o bmelt-grounded.o calceps2-0.2.o  \ 
     132routines_communes = tracebug.o ablation_bord.o \ 
     133        icetemp_declar_mod.o bmelt-grounded.o calceps2-0.2.o  \ 
    137134        diffusiv-polyn-0.6.o \ 
    138135        flottab-rescue-1.o flottab2-0.7.o   \ 
    139         interface_advec_icetemp.o \ 
    140         interface_prop_th_icetemp.o \ 
    141         interface_Qprod_icetemp.o \ 
    142         interface_temp_col.o \ 
    143         interface_icetempmod.o \ 
     136        advec_icetemp.o \ 
     137        prop_th_icetemp.o \ 
     138        Qprod_icetemp.o \ 
     139        temp_col.o \ 
     140        icetemp_mod.o \ 
    144141        initial-0.3.o initial-phy-2.o initial2-0.4.o \ 
    145142        lastoutput-0.2.o lect-eis.o  readinput.o lect_datfile.o write_datfile.o \ 
    … …  
    151148        tab-litho-0.3.o taubed-0.3.o \ 
    152149        velocities-polyn-0.3.o \ 
    153         vtoutput-0.2.o  courbures.o  
     150        vtoutput-0.2.o  courbures.o 
     151 
     152# allocate_types.o \ 
    154153#       vtoutput-0.2.o steps_time_loop.o courbures.o  
    155154 
    … …  
    471470        $(FT) 3D-physique-gen_mod.f90 
    472471 
    473 # types Hassine 
    474472param_phy_mod.o: param_phy_mod.f90 
    475473        $(FT) param_phy_mod.f90 
    476  
    477 geom_type.o: Types/geom_type.f90 
    478         $(FT) Types/geom_type.f90 
    479  
    480 temperature_type.o: Types/temperature_type.f90 
    481         $(FT) Types/temperature_type.f90 
    482  
    483 ice_flow_type.o: Types/ice_flow_type.f90 
    484         $(FT) Types/ice_flow_type.f90 
    485  
    486 mask_flgz_type.o: Types/mask_flgz_type.f90 
    487         $(FT) Types/mask_flgz_type.f90 
    488  
    489 deformation_type.o: Types/deformation_type.f90 
    490         $(FT) Types/deformation_type.f90 
    491  
    492 autre_pr_temp_type.o: Types/autre_pr_temp_type.f90 
    493         $(FT) Types/autre_pr_temp_type.f90 
    494  
    495 step_temp_type.o: Types/step_temp_type.f90 
    496         $(FT) Types/step_temp_type.f90 
    497  
    498 allocate_types.o:  Types/allocate_types.f90 
    499         $(FT) Types/allocate_types.f90 
    500  
    501 # fin types Hassine 
    502474 
    503475lect-anteis_mod.o : Ant40_files/lect-anteis_mod.f90 
    … …  
    13201292 
    13211293 
    1322 # Routines Temperature Hassine 
    1323  
    1324 icetemp_declar-0.5_mod.o : Temperature-routines/icetemp_declar-0.5_mod.f90 
    1325         $(FT) Temperature-routines/icetemp_declar-0.5_mod.f90 
    1326  
    1327 interface_advec_icetemp.o :  Temperature-routines/interface_advec_icetemp.f90 
    1328         $(FT) Temperature-routines/interface_advec_icetemp.f90 
    1329  
    1330 interface_icetempmod.o :  Temperature-routines/interface_icetempmod.f90 
    1331         $(FT) Temperature-routines/interface_icetempmod.f90 
    1332  
    1333 interface_prop_th_icetemp.o:  Temperature-routines/interface_prop_th_icetemp.f90 
    1334         $(FT) Temperature-routines/interface_prop_th_icetemp.f90 
    1335  
    1336 interface_Qprod_icetemp.o:  Temperature-routines/interface_Qprod_icetemp.f90 
    1337         $(FT) Temperature-routines/interface_Qprod_icetemp.f90 
    1338  
    1339 interface_temp_col.o :  Temperature-routines/interface_temp_col.f90 
    1340         $(FT) Temperature-routines/interface_temp_col.f90 
     1294# Routines Temperature 
     1295 
     1296icetemp_declar_mod.o : Temperature-routines/icetemp_declar_mod.f90 
     1297        $(FT) Temperature-routines/icetemp_declar_mod.f90 
     1298 
     1299advec_icetemp.o :  Temperature-routines/advec_icetemp.f90 
     1300        $(FT) Temperature-routines/advec_icetemp.f90 
     1301 
     1302icetemp_mod.o :  Temperature-routines/icetemp_mod.f90 
     1303        $(FT) Temperature-routines/icetemp_mod.f90 
     1304 
     1305prop_th_icetemp.o:  Temperature-routines/prop_th_icetemp.f90 
     1306        $(FT) Temperature-routines/prop_th_icetemp.f90 
     1307 
     1308Qprod_icetemp.o:  Temperature-routines/Qprod_icetemp.f90 
     1309        $(FT) Temperature-routines/Qprod_icetemp.f90 
     1310 
     1311temp_col.o :  Temperature-routines/temp_col.f90 
     1312        $(FT) Temperature-routines/temp_col.f90 
    13411313 
    13421314# Routines Agrif 
    … …  
    15561528        $(mod_tracers) \ 
    15571529        $(mod_ell) $(Liste_ANT15-LBq) \ 
    1558         $(mod_post_geo) \ 
    15591530        $(diagnoshelf) \ 
    15601531        $(Liste_Netcdf) \ 
    … …  
    15681539        $(mod_ell) \ 
    15691540        $(Liste_ANT15-LBq) \ 
    1570         $(mod_post_geo) \ 
    15711541        $(diagnoshelf) $(Liste_Netcdf) \ 
    15721542        $(routines_communes) steps_time_loop.o $(routine_elliptiques) $(NCDF_LIB)  $(MKL_LIB) 
    … …  
    15781548        $(mod_tracers) \ 
    15791549        $(mod_ell) $(Liste_ANT15-LBq) \ 
    1580         $(mod_post_geo) \ 
    15811550        $(diagnoshelf) \ 
    15821551        $(Liste_Netcdf) \ 
    … …  
    15901559        $(mod_ell) \ 
    15911560        $(Liste_ANT15-LBq) \ 
    1592         $(mod_post_geo) \ 
    15931561        $(diagnoshelf) $(Liste_Netcdf) \ 
    15941562        $(routines_communes) steps_time_loop.o $(routine_elliptiques) $(NCDF_LIB)  $(MKL_LIB) 
    … …  
    16131581        $(mod_tracers) \ 
    16141582        $(mod_ell) $(Liste_GrIce2sea) \ 
    1615         $(mod_post_geo) \ 
    16161583        $(diagnoshelf) \ 
    16171584        $(Liste_Netcdf) \ 
    … …  
    16251592        $(mod_tracers) \ 
    16261593        $(mod_ell) $(Liste_GrIce2sea) \ 
    1627         $(mod_post_geo) \ 
    16281594        $(diagnoshelf) \ 
    16291595        $(Liste_Netcdf) \ 
    … …  
    16371603        $(mod_tracers) \ 
    16381604        $(mod_ell) $(Liste_GrIce2sea) \ 
    1639         $(mod_post_geo) \ 
    16401605        $(diagnoshelf) \ 
    16411606        $(Liste_Netcdf) \ 
    … …  
    16501615        $(mod_tracers) \ 
    16511616        $(mod_ell) $(Liste_GrIce2sea) \ 
    1652         $(mod_post_geo) \ 
    16531617        $(diagnoshelf) \ 
    16541618        $(Liste_Netcdf) \ 
    … …  
    16621626        $(mod_no_tracers) \ 
    16631627        $(mod_ell) $(Liste_hemin40) \ 
    1664         $(mod_post_geo) \ 
    16651628        $(diagnoshelf) \ 
    16661629        $(Liste_Netcdf) \ 
    … …  
    16761639        $(mod_no_tracers) \ 
    16771640        $(mod_ell) $(Liste_hemin40) \ 
    1678         $(mod_post_geo) \ 
    16791641        $(diagnoshelf) \ 
    16801642        $(Liste_Netcdf) \ 
    … …  
    16881650        $(mod_no_tracers) \ 
    16891651        $(mod_ell) $(Liste_hemin15) \ 
    1690         $(mod_post_geo) \ 
    16911652        $(diagnoshelf) \ 
    16921653        $(Liste_Netcdf) \ 
    … …  
    17011662        $(mod_no_tracers) \ 
    17021663        $(mod_ell) $(Liste_hemin15) \ 
    1703         $(mod_post_geo) \ 
    17041664        $(diagnoshelf) \ 
    17051665        $(Liste_Netcdf) \ 
    … …  
    17131673        $(mod_no_tracers) \ 
    17141674        $(mod_ell) $(Liste_greeneem15) \ 
    1715         $(mod_post_geo) \ 
    17161675        $(diagnoshelf) \ 
    17171676        $(Liste_Netcdf) \ 
    … …  
    17261685        $(mod_no_tracers) \ 
    17271686        $(mod_ell) $(Liste_greeneem15) \ 
    1728         $(mod_post_geo) \ 
    17291687        $(diagnoshelf) \ 
    17301688        $(Liste_Netcdf) \ 

  • trunk/SOURCES/initial-0.3.f90

    r11 r24  
    1212!! @note    - use module3D_phy 
    1313!! @note    - use module_choix 
    14 !! @note    - use interface_icetempmod 
    15 !! @note    - use interface_prop_th_icetemp 
     14!! @note    - use icetempmod 
    1615!! @note    - use diagno_mod  
    1716!! @note    - use resolmeca_SIA_L1 
    … …  
    2524  use module3D_phy 
    2625  use module_choix 
    27   use interface_icetempmod 
    28   use interface_prop_th_icetemp 
     26  use icetempmod 
    2927  use diagno_mod  
    3028  use resolmeca_SIA_L1 
    … …  
    8886  if (geoplace(1:5).ne.'mism3') then 
    8987     call init_icetemp(Num_rep_42)        ! initialise le calcul de temperature 
    90      call thermal_prop_icetemp(T,TPMP,H,debug_3D)  ! appel aux proprietes thermiques 
     88     call thermal_prop_icetemp            ! appel aux proprietes thermiques 
    9189     call init_bmelt                      ! Appel aux tableaux d'initialisation des ice shelves 
    9290 
    … …  
    137135  call init_calving 
    138136 
    139   call allocate_types  ! cette routine alloue tous les types derive  
    140137! 
    141138!------------------------------------------------------------------------------------- 

  • trunk/SOURCES/main3D-0.4-40km.f90

    r19 r24  
    5555!! @note use module_choix  
    5656!! @note use flottab_mod 
    57 !! @note use interface_icetempmod 
     57!! @note use icetempmod 
    5858!! @note use sorties_ncdf_grisli 
    5959!! @note use diagno_mod  
    … …  
    8585!! @note use module_choix  
    8686!! @note use flottab_mod 
    87 !! @note use interface_icetempmod 
     87!! @note use icetempmod 
    8888!! @note use sorties_ncdf_grisli 
    8989!! @note use diagno_mod  
    … …  
    107107  !  de declaration des packages 
    108108  use flottab_mod 
    109   use interface_icetempmod 
     109  use icetempmod 
    110110  use sorties_ncdf_grisli 
    111111  use diagno_mod  
    … …  
    147147  !  de declaration des packages 
    148148  use flottab_mod 
    149   use interface_icetempmod 
     149  use icetempmod 
    150150  use sorties_ncdf_grisli 
    151151  use util_recovery 

  • trunk/SOURCES/steps_time_loop.f90

    r19 r24  
    1616  use module3d_phy 
    1717 
    18   use geom_type 
    19   use temperature_type 
    20   use ice_flow_type 
    21   use mask_flgz_type 
    22   use deformation_type 
    23   use autre_pr_temp_type 
    24  
    2518  use module_choix ! module de choix du type de run 
    2619  !  module_choix donne acces a tous les modules 
    … …  
    2922  use sorties_ncdf_grisli 
    3023  use flottab_mod 
    31   use interface_icetempmod 
     24  use icetempmod 
    3225  use diagno_mod  
    3326  !  use track_debug  
    … …  
    199192  !  module_choix donne acces a tous les modules 
    200193  !  de declaration des packages 
    201   use interface_icetempmod 
     194  use icetempmod 
    202195  use sorties_ncdf_grisli 
    203196  use flottab_mod 
    … …  
    205198  use resolmeca_SIA_L1 
    206199!  use track_debug  
    207  
    208   use geom_type 
    209   use temperature_type 
    210   use ice_flow_type 
    211   use mask_flgz_type 
    212   use deformation_type 
    213   use autre_pr_temp_type 
    214   use step_temp_type 
    215200 
    216201 
    … …  
    303288 
    304289     calc_temp:     if ((nt.gt.2).and.(geoplace(1:5).ne.'mism3')) then 
    305         ! Hassine ca va etre modifier bien sur apres: l appel sera beaucoup moin long que ca  
    306         call Update_geom(Geom_g)                
    307         call Update_Temperature(Temperature_g) 
    308         call Update_Ice_flow(Ice_flow_g) 
    309         call Update_Mask_flgz(Mask_flgz_g) 
    310         call Update_Deformation(Deformation_g) 
    311         call Update_autre_pr_temp(Autre_pr_temp_g) 
    312         call Update_step_temp(step_temp_g) 
    313  
    314290        if (itracebug.eq.1)  call tracebug('avant appel icetemp') 
    315  
    316         call icetemp(geom_g,temperature_g, & 
    317              Ice_flow_g,Mask_flgz_g,& 
    318              deformation_g,Autre_pr_temp_g,& 
    319              step_temp_g,num_tracebug,debug_3D) 
    320  
    321         ! update arrays T, Tpmp 
    322         T= Temperature_g%Temperature  
    323         Tpmp= Temperature_g%Temp_melting 
    324         Tbdot= Autre_pr_temp_g%Tbdot 
    325         Bmelt= Autre_pr_temp_g%BMELT 
    326         Ibase=Autre_pr_temp_g%IBASE 
    327         Phid=Autre_pr_temp_g%Phid 
    328  
    329         !      write(6,*) time,'call icetemp' 
     291        call icetemp 
    330292        ! subroutines pour le calcul de la fusion basale 
    331293        call bmeltshelf 
Note: See TracChangeset for help on using the changeset viewer.