Changeset 77 for branches/iLoveclim/SOURCES/Ant40_files/bmelt-ant-regions_mod.f90

Timestamp:

07/05/16 11:20:52 (8 years ago)

Author:

dumas

Message:

Merge branche iLOVECLIM sur rev 76

Location:

branches/iLoveclim

Files:

• . (modified) (1 prop)
• SOURCES/Ant40_files/bmelt-ant-regions_mod.f90 (modified) (9 diffs)

branches/iLoveclim/SOURCES/Ant40_files/bmelt-ant-regions_mod.f90

@@ -14 +14 @@
 MODULE  BMELT_ANT_REGIONS         ! cat juillet 2005
 
-  USE module3D_phy
-
-
+  use module3D_phy
+  use netcdf
+  use io_netcdf_grisli
+  
   implicit none
 
@@ -26 +27 @@
   !< utilises pour moduler la fusion sous le shelf
   integer, dimension(10) :: region      !< pour écrire dans le fichier param
-  character(len=30),dimension(10) :: regname !< nom des régions 
+  character(len=30),dimension(10) :: regname !< nom des régions
+  character(len=100) :: file_number_shelves ! fichier avec les zones ice-shelves
   real :: bsupshelf
   integer :: grd                        !< pour une sortie 
@@ -38 +40 @@
   real :: bmelt_talus,bmgrz_talus       !< Au dela du talus continental
   real :: bmelt_coef                    !< Coef pour corriger les valeurs std de bmelt (lu dans namelist)
+  real, parameter :: depth_talus=-2500. !< Profondeur du talus continental
 
 CONTAINS
@@ -46 +49 @@
     ! Cette routine fait l'initialisation pour la fusion basale.
     ! Elle est appelée par inputfile-vec-0.5.f90
-
+        
+    real*8, dimension(:,:),   pointer    :: tab               !< tableau 2d real pointer
 
     namelist/bmelt_ant_reg/bmelt_Ross,bmgrz_Ross,bmelt_FRis,bmgrz_FRis,  & 
          bmelt_Amery,bmgrz_Amery,bmelt_PIG,bmgrz_PIG,  &
          bmelt_Pen,bmgrz_Pen,bmelt_other,bmgrz_other,  &       
-         bmelt_talus,bmgrz_talus,bmelt_coef      
+         bmelt_talus,bmgrz_talus,bmelt_coef,file_number_shelves      
 
 
@@ -78 +82 @@
     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,*)'/'                            
     write(num_rep_42,428) '! bmelt_Ross & bmgrz_Ross    :  fonte basale secteur Ross'
@@ -87 +92 @@
     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,*)
 
 
-    !    ecriture dans le fichier parametres
-    write(num_rep_42,*)'fusion basale sous les ice shelves : bmelt-ant-regions_mod'
-    write(num_rep_42,*)'----------------------------------------------------------'
-
-
-    ! lecture du fichier contenant les distances au talu continental (m)
-    open(88,file=TRIM(DIRNAMEINP)//'distance_talu-40km.xy')
-
-    do j=1,ny
-       do i=1,nx
-          read(88,'(i3,1x,i3,1x,f10.2)') k,k,dist_talu(i,j)
-       enddo
-    enddo
-    close(88)
-
-    ! les bords de la grille sont mis a dist_talus=0 de force
-
-
-    dist_talu(1,:)=0.   !a gauche
-    dist_talu(2,:)=0.
-
-    dist_talu(nx-1,:)=0. !  a droite 
-    dist_talu(nx,:)=0.  
-
-    dist_talu(:,1)=0.    ! en bas   
-    dist_talu(:,2)=0.
-
-    dist_talu(:,ny-1)=0.    ! en haut
-    dist_talu(:,ny)=0.
-
-
-    debug_3D(:,:,32)=dist_talu(:,:)
 
     ! lecture du fichier contenant les types de shelves Ronne-Flichner
     !  ->1, Ross -> 2 , Amery -> 3, PIG-> 4, les petits shelves au dessus de PIG -> 5
-
-    open(88,file=TRIM(DIRNAMEINP)//'numer-ice-shelves-juil07.dat',status='OLD')
     typeshelf(:,:)=100
-    do k=1,nx*ny
-       read(88,*,end=36) i,j,typeshelf(i,j)
-    end do
-36  close(88)
-
-
+    file_number_shelves=TRIM(DIRNAMEINP)//trim(file_number_shelves)
+!    call lect_input(1,'z',1,typeshelf,file_number_shelves,file_ncdf)
+        call Read_Ncdf_var('z',file_number_shelves,tab)
+        typeshelf(:,:)  = tab(:,:)
+!    open(88,file=TRIM(DIRNAMEINP)//'numer-ice-shelves-juil07.dat',status='OLD')
+!    open(88,file=TRIM(file_number_shelves),status='OLD')
+!       do j=1,ny
+!               do i=1,nx
+!                       read(88,*) typeshelf(i,j)
+!       enddo   
+!       enddo
+!       close(88)
 
     region(:)=0
@@ -147 +125 @@
        do i=1,nx
 
-          talus:    if (dist_talu(i,j).GT.5000.) then
+          talus:    if (Bsoc(i,j).GT.depth_talus) then
              if (nint(typeshelf(i,j)).eq.1) then ! Ronne-Filchner FRis
                 bmshelf(i,j)=bmelt_FRis         
@@ -227 +205 @@
        enddo
     enddo bms_init
+    
 
     bmshelf(:,:)=bmshelf(:,:)*bmelt_coef
@@ -247 +226 @@
        do i=2,nx-1
 
-          talus_nochange : if (typeshelf(i,j).lt.10) then 
+          talus_nochange : if (Bsoc(i,j).GT.depth_talus) then 
              coef_talus = coefbmshelf
           else
-             coef_talus = 1         ! pas de changement au dela du talus continental
+             coef_talus = 1         ! pas de changement au dela du talus continental               
           endif talus_nochange