Changeset 20 for trunk/SOURCES

Timestamp:

09/14/15 16:53:18 (9 years ago)

Author:

dumas

Message:

climat-forcage-insolation_mod.f90 : Lecture du nombre de ligne du fichier CO2 a partir de la premiere ligne du fichier de forcage (plus de modification a faire en dur dans le code)

File:

• trunk/SOURCES/climat-forcage-insolation_mod.f90 (modified) (4 diffs)

trunk/SOURCES/climat-forcage-insolation_mod.f90

@@ -26 +26 @@
 integer,parameter :: mois=12
 integer,parameter :: ntr=2      ! nb de snapshots selon les paramètres orbitaux
-integer,parameter :: gtr=3      ! nb de snapshots selon l'état de la calotte
+integer,parameter :: gtr=4      ! nb de snapshots selon l'état de la calotte
 integer,parameter :: ctr=4      ! nb de snapshots selon le CO2
 real,dimension(ctr) :: Seuil_haut
@@ -34 +34 @@
 character(len=12) :: ice_int='yes'         !Choose between yes and no
 character(len=12) :: interpolation_ice='nonlinear'   !Choose between linear and nonlinear
-integer :: rayon_interpolation_ice=300  !(In km) Choose a multiple of nx.
+integer :: rayon_interpolation_ice=60  !(In km) Choose a multiple of nx.
 !real :: poids_nb_points=0.5
 !real :: poids_proximity=0.5
 character(len=15) :: CO2_function='variable'   !Choose between variable and constant
-integer,parameter :: nb_slope_dis_CO2=6 !10              !Number of slope discontinuities + 1 in the text file containing the CO2 function (nbr de ligne fichier CO2)
+integer :: nb_slope_dis_CO2  ! nbr de ligne fichier CO2
 real :: CO2_value=980.
 character(len=12) :: interpolation_CO2='logarithmic'   !Choose between linear and logarithmic
@@ -103 +103 @@
 
 !interpolation sur le CO2
-double precision,dimension(nb_slope_dis_CO2) :: CO2_time,CO2_val
+double precision,allocatable :: CO2_time(:), CO2_val(:) ! depend du nbr de points du fichier CO2
 real,dimension(nx,ny,mois) :: Tm_time_fin_3             ! temperature mensuelle au temps time (non corrige altitude)
 real,dimension(nx,ny,mois) :: Pm_time_fin_3             ! precipitation mensuelle au temps time
@@ -389 +389 @@
 case('variable')
    print *,'The CO2 will vary with time'
-
    open(75,file=trim(dirnameinp)//'forcing/'//trim(co2_file))
+   read(75,*) nb_slope_dis_CO2
+   allocate(CO2_time(nb_slope_dis_CO2),stat=error)
+   if (error.ne.0) then
+      print *,'error: could not allocate memory for array CO2_time'
+      stop
+   endif
+   allocate(CO2_val(nb_slope_dis_CO2),stat=error)
+   if (error.ne.0) then
+      print *,'error: could not allocate memory for array CO2_val'
+      stop
+   endif
    do w=1,nb_slope_dis_CO2
       read(75,*) CO2_time(w),CO2_val(w)
    end do
    close(75)
-   
    !print *,'CO2_time = ',CO2_time(:),'CO2_val = ',CO2_val(:)
-
 end select