Changeset 75 for trunk/SOURCES/Ant40_files/lect-anteis_mod.f90

Timestamp:

06/28/16 18:04:12 (8 years ago)

Author:

dumas

Message:

OpenMP parallelization in remplimat-shelves-tabTu.f90 and simplified compilation options

File:

• trunk/SOURCES/Ant40_files/lect-anteis_mod.f90 (modified) (4 diffs)

trunk/SOURCES/Ant40_files/lect-anteis_mod.f90

@@ -37 +37 @@
 
     integer :: ios
+    ! pour les lectures ncdf
+    real*8, dimension(:,:),   pointer    :: tab               !< tableau 2d real pointer
 
     namelist/topo_file/topo_ref,topo_dep,grid_topo,ghf_fich
@@ -101 +103 @@
     ! voir init_iso et a avoir une surface de reference pour les temperatures
     ! lecture adaptee aux fichiers ZBL.dat ou netcdf ou grd
-    call lect_input(1,'Bsoc',1,Bsoc0,topo_ref,file_ncdf)    ! socle
-    call lect_input(1,'S',1,S0,topo_ref,file_ncdf)          ! surface
-    call lect_input(1,'H',1,H0,topo_ref,file_ncdf)          ! epaisseur
+!    call lect_input(1,'Bsoc',1,Bsoc0,topo_ref,file_ncdf)    ! socle
+!    call lect_input(1,'S',1,S0,topo_ref,file_ncdf)          ! surface
+!    call lect_input(1,'H',1,H0,topo_ref,file_ncdf)          ! epaisseur
+
+! lecture pour eviter plantage avec compile -O0
+     call Read_Ncdf_var('Bsoc',topo_ref,tab)
+     Bsoc0(:,:)  = tab(:,:)
+     call Read_Ncdf_var('S',topo_ref,tab)
+     S0(:,:)  = tab(:,:)
+     call Read_Ncdf_var('H',topo_ref,tab)
+     H0(:,:)  = tab(:,:)
 !cdc correction point pole sud :
 !    S0(71,71)=(S0(71,70)+S0(71,72)+S0(70,71)+S0(72,71))/4.
@@ -127 +137 @@
     !     ---------------------------
     ! lecture adaptee aux fichiers ZBL.dat ou netcdf ou grd
-    call lect_input(1,'Bsoc',1,Bsoc,topo_dep,file_ncdf)    ! socle
-    call lect_input(1,'S',1,S,topo_dep,file_ncdf)          ! surface
-    call lect_input(1,'H',1,H,topo_dep,file_ncdf)          ! epaisseur
+!    call lect_input(1,'Bsoc',1,Bsoc,topo_dep,file_ncdf)    ! socle
+!    call lect_input(1,'S',1,S,topo_dep,file_ncdf)          ! surface
+!    call lect_input(1,'H',1,H,topo_dep,file_ncdf)          ! epaisseur
+    ! lecture pour eviter plantage avec compile -O0
+     call Read_Ncdf_var('Bsoc',topo_dep,tab)
+     Bsoc(:,:)  = tab(:,:)
+     call Read_Ncdf_var('S',topo_dep,tab)
+     S(:,:)  = tab(:,:)
+     call Read_Ncdf_var('H',topo_dep,tab)
+     H(:,:)  = tab(:,:)
+    
 !    S(71,71)=(S(71,70)+S(71,72)+S(70,71)+S(72,71))/4.
 !    Bsoc(71,71)=(Bsoc(71,70)+Bsoc(71,72)+Bsoc(70,71)+Bsoc(72,71))/4.
@@ -242 +260 @@
 !    close(88)
 
-    call lect_input(1,'ghf',1,ghf,ghf_fich,file_ncdf)
+!    call lect_input(1,'ghf',1,ghf,ghf_fich,file_ncdf)
+    ! pour eviter plantage -O0
+    call Read_Ncdf_var('ghf',ghf_fich,tab)
+    ghf(:,:)  = tab(:,:)
 
     ! pour passer les flux des mW/m2 au J/m2/an