Changeset 34


Ignore:
Timestamp:
01/19/16 15:53:01 (8 years ago)
Author:
dumas
Message:

Ajout du module Ant40_files/bmelt-ant-regions-oce_mod.f90 pour forcer la fonte basale avec des fichiers de temperature oceanique (M2 Ning) | lecture fichier netcdf pour Antarctique | Ajout dans les sortie de ghf : flux geothermique | climat_GrIce2sea_years_mod.f90 : lecture de smb et Tann au lieu de z | correction bug sorties Netdf dans initial-0.3.f90

Location:
trunk/SOURCES
Files:
2 added
12 edited

Legend:

Unmodified
Added
Removed

  • trunk/SOURCES/Ant40_files/lect-anteis_mod.f90

    r27 r34  
    115115!    endwhere 
    116116 
     117 
     118    where (BSOC(:,:).LT.-9999.) 
     119       BSOC(:,:)=-9999. 
     120    endwhere 
     121 
    117122    sealevel0=0.  ! voir a passer dans le fichier parametre 
    118123    S0(:,:)=max(S0(:,:),sealevel0)   ! pour etre au niveau des mers : ATTENTION si SEALEV <0 
    … …  
    231236    open(88,file=ghf_fich) 
    232237 
    233     write(num_rep_42,*) 'flux geothermique Shapiro : ', TRIM(DIRNAMEINP)//'ijphi-40km-ant.dat' 
    234  
    235     do k=1,nx*ny 
    236        read(88,*) i,j,ghf(i,j) 
    237     end do 
    238     close(88) 
     238!    write(num_rep_42,*) 'flux geothermique Shapiro : ', TRIM(DIRNAMEINP)//'ijphi-40km-ant.dat' 
     239!    do k=1,nx*ny 
     240!       read(88,*) i,j,ghf(i,j) 
     241!    end do 
     242!    close(88) 
     243 
     244    call lect_input(1,'ghf',1,ghf,ghf_fich,file_ncdf)    ! socle 
     245 
     246 
    239247 
    240248    ! pour passer les flux des mW/m2 au J/m2/an       

  • trunk/SOURCES/Ant40_files/module_choix-antar40-0.4.f90

    r27 r34  
    4545 
    4646!--------------Lecture climat forcage--------------- 
    47 !use climat_forcage_mois_mod ! forcage mensuel GCM 1 Snapshot Fev 2015 
    48 use climat_perturb_mod   ! climat perturbe a reverifier Dec 2015 
     47use climat_forcage_mois_mod ! forcage mensuel GCM 1 Snapshot Fev 2015 
     48!use climat_perturb_mod   ! climat perturbe a reverifier Dec 2015 
     49!use climat_Grice2sea_years_mod  ! climat force par fichier SMB directement (grice2sea) 
    4950 
    5051!use climat_forcage_mod 
    … …  
    105106 
    106107!--------------Fusion basale------------------------ 
    107 use bmelt_ant_regions ! pour l'Antarctique avec régions 
     108!use bmelt_ant_regions ! pour l'Antarctique avec régions 
     109use bmelt_ant_regions_oce ! pour Ant region et temp ocean 
    108110!use bmelt_nor_regions ! pour le nord avec régions 
    109111!use bmelt_nor_depth   ! pour le nord avec profondeur d'eau 

  • trunk/SOURCES/Fichiers-parametres/anteis1_LISTE-VAR-NETCDF.dat

    r27 r34  
    1441441       1       1                        
    145145----------------------------------------------------------- 
     146ghf 
     14723 
     1481       1      -1 
     149----------------------------------------------------------- 
    146150phid             
    14715124       

  • trunk/SOURCES/Fichiers-parametres/greeneem15_LISTE-VAR-NETCDF.dat

    r19 r34  
    1441441       1       1                        
    145145----------------------------------------------------------- 
     146ghf 
     14723 
     1481       1      -1 
     149----------------------------------------------------------- 
    146150phid             
    14715124       

  • trunk/SOURCES/Fichiers-parametres/hemin15_LISTE-VAR-NETCDF.dat

    r14 r34  
    1441441       1       1                        
    145145----------------------------------------------------------- 
     146ghf 
     14723 
     1481       1      -1 
     149----------------------------------------------------------- 
    146150phid             
    14715124       

  • trunk/SOURCES/Fichiers-parametres/hemin40_LISTE-VAR-NETCDF.dat

    r12 r34  
    1441441       1       1                        
    145145----------------------------------------------------------- 
     146ghf 
     14723 
     1481       1      -1 
     149----------------------------------------------------------- 
    146150phid             
    14715124       

  • trunk/SOURCES/GrIce2sea_files/climat_GrIce2sea_years_mod.f90

    r33 r34  
    1616 
    1717 
    18 use module3d_phy,only: nx,ny,S,S0,H,Bsoc,acc,abl,BM,Tann,Tjuly,Ts,time,dt,num_param,num_rep_42,dirnameinp 
     18use module3d_phy,only: nx,ny,S,S0,H,Bsoc,acc,abl,BM,Tann,Tjuly,Ts,time,dt,num_param,num_rep_42,dirnameinp,coefbmshelf 
    1919!use lect_climref_Ice2sea 
    2020use netcdf 
    … …  
    104104  smb_file  = trim(dirnameinp)//trim(smb_file) 
    105105 
    106   call Read_Ncdf_var('z',smb_file,tab) 
     106  call Read_Ncdf_var('smb',smb_file,tab) 
    107107 
    108108!  call lect_input(3,'smb',1,bm,smb_file,trim(dirnameinp)//trim(runname)//'.nc') 
    … …  
    132132  temp_annual_file = trim(dirnameinp)//trim(temp_annual_file) 
    133133 
    134   call Read_Ncdf_var('z',temp_annual_file,tab) 
     134  call Read_Ncdf_var('Tann',temp_annual_file,tab) 
    135135  Tann(:,:)  = tab(:,:) 
    136136!  call lect_input(3,'Tann',1,Tann,temp_annual_file,trim(dirnameinp)//trim(runname)//'.nc') 
    … …  
    284284     call init_grad_smb 
    285285  endif 
     286! coefbmshelf coefficient pour la fusion basale sous les ice shelves 
     287! ici pas d'evolution temporelle donc coefbmshelf est fixe 
     288  coefbmshelf=1. 
    286289 
    287290  return 

  • trunk/SOURCES/Makefile

    r33 r34  
    1010 
    1111# Choice for compiler (ifort = 1, gfortran = 0) 
    12 ifort ?= 0 
     12ifort ?= 1 
    1313 
    1414# Choice for librairies (mkl_c = 1 = MKL,  = 0 = BLAS) 

  • trunk/SOURCES/Makefile.grisli.inc

    r29 r34  
    214214        lect-anteis_mod.o \ 
    215215        bmelt-ant-regions_mod.o \ 
     216        bmelt-ant-regions-oce_mod.o \ 
    216217        fake-routines-ant_mod.o module_choix-antar40-0.4.o \ 
    217218        track_ant40_mod.o 
    … …  
    515516 
    516517bmelt-ant-regions_mod.o : Ant40_files/bmelt-ant-regions_mod.f90 
    517         $(FT)  Ant40_files/bmelt-ant-regions_mod.f90 
     518        $(FT) Ant40_files/bmelt-ant-regions_mod.f90 
     519 
     520bmelt-ant-regions-oce_mod.o : Ant40_files/bmelt-ant-regions-oce_mod.f90 
     521        $(FT) Ant40_files/bmelt-ant-regions-oce_mod.f90 
    518522 
    519523bmelt-nor-regions_mod.o : Heminord_files/bmelt-nor-regions_mod.f90 

  • trunk/SOURCES/Netcdf-routines/Description_Variables.dat

    r4 r34  
    225225"celsius"        
    226226"Basal temperature " 
     227------------------------------------------------------ 
     228"xxx" 
     22923      "ghf"  1 
     230"o" 
     231"geothermal_heat_flux" 
     232"geothermal_heat_flux" 
     233"m.w/m²" 
     234"Geothermal heat flux" 
    227235------------------------------------------------------ 
    228236"xxx" 

  • trunk/SOURCES/initial-0.3.f90

    r29 r34  
    7171  endif 
    7272 
    73 !cdc call sortie_ncdf_cat 
     73  call sortie_ncdf_cat 
    7474 
    7575  ! calcul de Hmx et Hmy -> shift=-1, dim=1 -> H(i-1,j) 

  • trunk/SOURCES/lect-eis.f90

    r29 r34  
    1818 
    1919  implicit none 
    20   character(len=17)  :: DIRNAMEINP 
     20  character(len=20)  :: DIRNAMEINP 
    2121  integer nx,ny,i,j,kdum 
    2222  real data(nx,ny) 
Note: See TracChangeset for help on using the changeset viewer.