Changeset 213

Timestamp:

10/02/18 18:18:58 (6 years ago)

Author:

dumas

Message:

back to previous version of bmelt-seuil-profondeur_mod.f90

File:

• trunk/SOURCES/bmelt-seuil-profondeur_mod.f90 (modified) (10 diffs)

trunk/SOURCES/bmelt-seuil-profondeur_mod.f90

@@ -22 +22 @@
 
   use module3d_phy
-  use netcdf
-  use io_netcdf_grisli
+
 
   implicit none
@@ -34 +33 @@
   real :: depth_talus                   !< profondeur de transition 
   real,dimension(nx,ny) ::  bmshelf     !< tableau fusion basale sous shelf
-  
-  
-! Pour le calcul du bmelt moyen par ice-shelf
-  integer,dimension(nx,ny) :: typeshelf !< Type de shelf Ronne->1 Ross ->2 ....
-  integer, parameter :: nb_reg = 18
-  real,dimension (nb_reg) :: bmelt_regions, bmgrz_regions
-  real :: bmelt_talus,bmgrz_talus       !< Au dela du talus continental
-  real :: bmelt_coef                    !< Coef pour corriger les valeurs std de bmelt (lu dans namelist mais pas utilisé ici)
-  integer,dimension(nb_reg):: i_bmelt   !< nbr de passage dans la boucle pour calcul de la moyenne du bmelt
-  integer,dimension(nb_reg):: i_bmgrz   !< nbr de passage dans la boucle pour calcul de la moyenne du bmgrz
-  
-  integer :: flag_dist                  !< Pas utilise ici mais pour compatibilite lecture namelist
-  
+
 
 contains
@@ -61 +48 @@
     ! Cette routine fait l'initialisation pour la fusion basale.
     ! Elle est appelée par inputfile-vec-0.5.f90
-    
-    ! Pour le calcul du bmelt moyen par ice-shelf
-    real*8, dimension(:,:),   pointer    :: tab               !< tableau 2d real pointer
-    character(len=100) :: file_number_shelves ! fichier avec les zones ice-shelves
-    character(len=100) :: file_bmelt_anom
-    real,dimension(nx,ny) :: dist_talus   ! distance au talus
-    character(len=100) :: file_dist_talus ! pas utilise ici mais pour compatibilite lecture namelist
-    real :: bmelt_dist0                   ! pas utilise ici mais pour compatibilite lecture namelist
-    
+
+
     namelist/bmelt_seuil/bm_grz,bmshelf_plateau,bmshelf_abysses,depth_talus
-    
-    namelist/bmelt_ant_reg_initmip/bmelt_regions,bmgrz_regions,  &
-      bmelt_talus,bmgrz_talus,bmelt_coef,file_number_shelves,   &
-      flag_dist,file_dist_talus,bmelt_dist0
 
     if (itracebug.eq.1)  call tracebug('entree dans init_bmelt de bmelt_seuil_prof')
@@ -111 +87 @@
 
     debug_3D(:,:,34)=bmshelf(:,:)
-    
-! Si on fait l'inversion on doit connaitre les regions pour calculer le taux de fonte par region    
-    if (ibmelt_inv.eq.1) then
-      rewind(num_param)        ! pour revenir au debut du fichier param_list.dat
-      read(num_param,bmelt_ant_reg_initmip)
-      ! formats pour les ecritures dans 42
-      write(num_rep_42,428)'!___________________________________________________________'
-      write(num_rep_42,428) '&bmelt-ant-regions-initmip                                 '
-      write(num_rep_42,*)
-      write(num_rep_42,'(A,18(f0.2,","))') 'bmelt_regions      = ', bmelt_regions(:)
-      write(num_rep_42,'(A,18(f0.2,","))') 'bmgrz_regions      = ', bmgrz_regions(:)
-      write(num_rep_42,*) 'bmelt_talus     = ', bmelt_talus
-      write(num_rep_42,*) 'bmgrz_talus     = ', bmgrz_talus
-      write(num_rep_42,*) 'bmelt_coef      = ', bmelt_coef
-      write(num_rep_42,*) 'file_numer_shelves = ', file_number_shelves
-      write(num_rep_42,*) 'flag_dist = ', flag_dist
-      write(num_rep_42,*) 'file_dist_talus = ', file_dist_talus
-      write(num_rep_42,*) 'bmelt_dist0 = ', bmelt_dist0
-      write(num_rep_42,*)'/'
-      write(num_rep_42,428) '! bmelt_regions    :  fonte basale sous shelves pour les 18 regions initMIP'
-      write(num_rep_42,428) '! bmgrz_regions    :  fonte basale grounding zone pour les 18 regions initMIP'
-      write(num_rep_42,428) '! bmelt_talus & bmgrz_talus  :  fonte basale apres talus cont'
-      write(num_rep_42,428) '! bmelt_coef                 :  coef fonte (1 pour conserver val)'
-      write(num_rep_42,428) '! file_numer_ice-shelves     : fichier zones ice shelves'
-!~     write(num_rep_42,428) '! flag_dist : flag pour bmelt fnct distance talus'
-!~     write(num_rep_42,428) '! file_dist_talus : fichier de distance talus'
-!~     write(num_rep_42,428) '! bmelt_dist0 : coef de bmelt au talus'
-      write(num_rep_42,*)
-      ! lecture du fichier contenant les types de shelves
-      ! 18 regions sur la base des bassins fournis par le projet ISMIP6-initmip antarctica
-      typeshelf(:,:)=100
-      file_number_shelves=TRIM(DIRNAMEINP)//trim(file_number_shelves)
-      call Read_Ncdf_var('z',file_number_shelves,tab)
-      typeshelf(:,:)  = tab(:,:)
-    endif
-    bmelt_regions(:)=0.
-    bmgrz_regions(:)=0.
-    i_bmelt(:)=0
-    i_bmgrz(:)=0
-
     return
   end subroutine init_bmelt
@@ -169 +105 @@
     integer :: ngr           ! nombre de voisins flottants
     real    :: coef_talus    ! pour ne pas changer la fusion au dessus de l'ocean profond
-    
-    ! Pour le calcul du bmelt moyen par ice-shelf
-    integer,dimension(nb_reg) :: nbr_pts_shelf, nbr_pts_grz ! nombre de points shelf et grz par zone
-    real,dimension(nb_reg) :: bmelt_regions_time, bmgrz_regions_time ! Pour le caclul du bemlt/bmgrz par zone
-    integer :: nregion ! pour calcul du bmelt sur toutes les regions
 
     if (itracebug.eq.1)  call tracebug('entree dans bmeltshelf de bmelt_seuil_prof')
-    
-    nbr_pts_shelf(:)=0
-    nbr_pts_grz(:)=0
-    bmelt_regions_time(:)=0.
-    bmgrz_regions_time(:)=0.
-    
+
     where (Bsoc0(:,:).lt.depth_talus)
         bmshelf(:,:)=bmshelf_abysses
@@ -206 +132 @@
 
             if (fbm(i,j)) then 
-              bmelt(i,j)=coef_talus*bmgrz(i,j)
+                bmelt(i,j)=coef_talus*bmgrz(i,j)
             endif
+
+
 
             ! ATTENTION le bloc suivant sert a determiner la fusion basale d'equilibre
@@ -216 +144 @@
                 !     bmelt(i,j)=corrbmelt(i,j)
                 corrbmelt(i,j)=corrbmelt(i,j)+hdot(i,j)*0.85
-                bmelt(i,j)=bmelt(i,j)+corrbmelt(i,j)
-                if (H(i,j).gt.50.) then ! on ne prend en compte que les points avec H > 50m
-                  if (fbm(i,j)) then ! points grounding line
-                    bmgrz_regions_time(typeshelf(i,j)) = bmgrz_regions_time(typeshelf(i,j)) + bmelt(i,j)
-                    nbr_pts_grz(typeshelf(i,j)) = nbr_pts_grz(typeshelf(i,j)) + 1
-                  else ! points shelf
-                    bmelt_regions_time(typeshelf(i,j)) = bmelt_regions_time(typeshelf(i,j)) + bmelt(i,j)
-                    nbr_pts_shelf(typeshelf(i,j)) = nbr_pts_shelf(typeshelf(i,j)) + 1
-                  endif
-                endif
+                bmelt(i,j)=bmelt(i,j)+corrbmelt(i,j)     
             endif
+
 
         else                   ! point posé, on compte le nombre de voisins flottants
@@ -237 +157 @@
             end if
 
+
             !   la fusion des points limites est une combinaison entre valeur posée et valeur flottante
             !   en fonction du nombre de points flottants
@@ -243 +164 @@
 
             if ((igrdline.eq.1).and.(ibmelt_inv.eq.1)) then
-                corrbmelt(i,j)=corrbmelt(i,j)+hdot(i,j)*0.85
+                corrbmelt(i,j)=corrbmelt(i,j)+hdot(i,j)  !*0.85
                 bmelt(i,j)=bmelt(i,j)+corrbmelt(i,j)
             endif
@@ -250 +171 @@
     end do
 
-    if ((ibmelt_inv.eq.1).and.(time.gt.10.)) then ! on ne prend pas en compte les 30 premiers pas de temps
-      do nregion=1,nb_reg
-        if (nbr_pts_shelf(nregion).gt.1) then
-          bmelt_regions(nregion) = bmelt_regions_time(nregion)/nbr_pts_shelf(nregion) + bmelt_regions(nregion)
-          i_bmelt(nregion) = i_bmelt(nregion)+1
-        endif
-        if (nbr_pts_grz(nregion).ge.1) then
-          bmgrz_regions(nregion) = bmgrz_regions_time(nregion)/nbr_pts_grz(nregion) + bmgrz_regions(nregion)
-          i_bmgrz(nregion) = i_bmgrz(nregion)+1
-        endif
-      enddo
-      print*,'time=',time
-      print*,'bmgrz des ',nb_reg,' regions'
-      print*,'Ross 15 16'
-      print*,bmgrz_regions_time(15)/nbr_pts_grz(15),nbr_pts_grz(15),bmgrz_regions_time(16)/nbr_pts_grz(16),nbr_pts_grz(16)
-      print*,'Ronny 14 18'
-      print*,bmgrz_regions_time(14)/nbr_pts_grz(14),nbr_pts_grz(14),bmgrz_regions_time(18)/nbr_pts_grz(18),nbr_pts_grz(18)
-      print*,'PIG 17 et Amery 5'
-      print*,bmgrz_regions_time(17)/nbr_pts_grz(17),nbr_pts_grz(17),bmgrz_regions_time(5)/nbr_pts_grz(5),nbr_pts_grz(5)
-      
-      if (time.eq.tend) then ! a le fin du run : calcul du bmelt et bmgrz moyen par secteur
-        do nregion=1,nb_reg
-          if (i_bmelt(nregion).gt.0) then
-            bmelt_regions(nregion)=bmelt_regions(nregion)/i_bmelt(nregion)
-          else
-            bmelt_regions(nregion)=0.
-          endif
-          if (i_bmgrz(nregion).gt.0) then
-            bmgrz_regions(nregion)=bmgrz_regions(nregion)/i_bmgrz(nregion)
-          else
-            bmgrz_regions(nregion)=0.
-          endif
-        enddo
-        print*,'Fin de la simule, time=',time
-        print*,'calcul des bmelt et bmgrz par inversion :'
-        print*,'bmelt des ',nb_reg,' regions'
-        print*,bmelt_regions(:)
-        print*,'bmgrz des ',nb_reg,' regions'
-        print*,bmgrz_regions(:)
-      endif
-    endif
-          
-    if ((igrdline.eq.1).and.(ibmelt_inv.eq.0)) then
-      where (i_Hp(:,:).eq.1)
+if ((igrdline.eq.1).and.(ibmelt_inv.eq.0)) then
+    where (i_Hp(:,:).eq.1)
         bmelt(:,:)=0.    ! hdot donne alors le -bmelt
-      endwhere
-    endif
-    
+    endwhere
+endif
+
 return
 end subroutine  bmeltshelf