Changeset 69 for trunk/SOURCES

Timestamp:

06/03/16 17:14:46 (8 years ago)

Author:

dumas

Message:

bmelt-ant-regions_mod.f90 read netcdf file

File:

• trunk/SOURCES/Ant40_files/bmelt-ant-regions_mod.f90 (modified) (7 diffs)

trunk/SOURCES/Ant40_files/bmelt-ant-regions_mod.f90

@@ -15 +15 @@
 
   use module3D_phy
-
+  use netcdf
+  use io_netcdf_grisli
+  
   implicit none
 
@@ -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,  & 
@@ -98 +102 @@
     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)
+!    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
@@ -119 +125 @@
        do i=1,nx
 
-          talus:    if (Bsoc(i,j).LT.2500.) 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         
@@ -199 +205 @@
        enddo
     enddo bms_init
+    
+    print*,'33,107bmlet',bmshelf(33,107),typeshelf(33,107),bsoc(33,107)
+        print*,'33,120bmlet',bmshelf(33,120),typeshelf(33,120),bsoc(33,120)
 
     bmshelf(:,:)=bmshelf(:,:)*bmelt_coef
@@ -219 +228 @@
        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