Changeset 220

Timestamp:

12/05/18 14:32:29 (6 years ago)

Author:

aquiquet

Message:

Further development of the new climat_forcage_mod: possibility to run transient OR equilibrium simulations using the snapshots provided (integer typerun).

File:

• trunk/SOURCES/climat_forcage_mod.f90 (modified) (10 diffs)

trunk/SOURCES/climat_forcage_mod.f90

@@ -38 +38 @@
   real,dimension(nx,ny) :: Zs        !< surface topography above sea level
 
+  integer :: typerun ! 1 for steady state using snap 1, 
+                     ! 2 for steady state using snap 2 ...
+                     ! 0 for transient
 
 ! S0 topo de ref correspondant a la climato est lu dans lect_topo 
@@ -188 +191 @@
 
 
-  namelist/clim_forcage/clim_ref_file,ttr,forcage_file1,forcage_file2,rappact,mincoefbmelt,maxcoefbmelt,filforc,lapserate
+  namelist/clim_forcage/clim_ref_file,ttr,forcage_file1,forcage_file2,typerun,rappact,mincoefbmelt,maxcoefbmelt,filforc,lapserate
 
   rewind(num_param)        ! pour revenir au debut du fichier param_list.dat
@@ -204 +207 @@
   write(num_rep_42,'(A,A)') 'forcage_file1 = ', forcage_file1
   write(num_rep_42,'(A,A)') 'forcage_file2 = ', forcage_file2
+  write(num_rep_42,'(A,A)') 'typerun = ', typerun
   write(num_rep_42,'(A,A)') 'rappact = ', rappact
   write(num_rep_42,*) 'mincoefbmelt = ', mincoefbmelt
@@ -234 +238 @@
   real,dimension(nx,ny,12) :: prmois ! precip sur topo Zs
   
-  
-  
-  if(time.le.tdate(1)) then   ! time avant le debut du fichier forcage
+  real (kind=kind(0.d0)) ::  timeloc
+  
+  if (typerun.eq.0) then
+     timeloc = time
+  else if ( (typerun .gt. 0) .and. (typerun .le. ubound(ttr,dim=1)) ) then
+     timeloc = ttr(typerun)
+  else
+     write(*,*) "Unknown type of simulation for climat_forcage, abort"
+     STOP
+  end if
+
+  if(timeloc.le.tdate(1)) then   ! time avant le debut du fichier forcage
      sealevel=spert(1)
      coeft=alphat(1)
@@ -242 +255 @@
      ift=1
 
-  else if (time.ge.tdate(nft)) then     ! time apres la fin du fichier forcage
+  else if (timeloc.ge.tdate(nft)) then     ! time apres la fin du fichier forcage
      sealevel=spert(nft)
      coeft=alphat(nft)
@@ -253 +266 @@
      do i=ift,nft-1  ! interpolation entre i et i+1
 
-        if((time.ge.tdate(i)).and.(time.lt.tdate(i+1))) then  ! cas general
+        if((timeloc.ge.tdate(i)).and.(timeloc.lt.tdate(i+1))) then  ! cas general
            sealevel=spert(i)+(spert(i+1)-spert(i))*    &
-                (time-tdate(i))/(tdate(i+1)-tdate(i))
+                (timeloc-tdate(i))/(tdate(i+1)-tdate(i))
 
            coeft=alphat(i)+(alphat(i+1)-alphat(i))*      &
-                (time-tdate(i))/(tdate(i+1)-tdate(i))
+                (timeloc-tdate(i))/(tdate(i+1)-tdate(i))
 
            coefp=alphap(i)+(alphap(i+1)-alphap(i))*      &
-                (time-tdate(i))/(tdate(i+1)-tdate(i))
+                (timeloc-tdate(i))/(tdate(i+1)-tdate(i))
                 
            ift=i
@@ -272 +285 @@
   slv(:,:)=sealevel
 
-  if(time.le.ttr(1)) then   ! time en dehors des limites du fichier forcage
+  if(timeloc.le.ttr(1)) then   ! time en dehors des limites du fichier forcage
     coeftime=0.
 !~     do j=1,ny
@@ -282 +295 @@
 !~     end do
     itr=1
-  else if (time.ge.ttr(ntr)) then        ! mis a 0
+  else if (timeloc.ge.ttr(ntr)) then        ! mis a 0
     coeftime=1.
     itr=ntr
@@ -295 +308 @@
      itr=1
      do l=itr,ntr-1    ! interpolation entre i et i+1
-        if((time.ge.ttr(l)).and.(time.lt.ttr(l+1))) then  ! test tranche
-           coeftime= (time-ttr(l))/(ttr(l+1)-ttr(l))
+        if((timeloc.ge.ttr(l)).and.(timeloc.lt.ttr(l+1))) then  ! test tranche
+           coeftime= (timeloc-ttr(l))/(ttr(l+1)-ttr(l))
 !           do j=1,ny
 !              do i=1,nx
@@ -345 +358 @@
 
 
-print*,'Tmois',Tmois(142,14,1)
-print*,'prmois',prmois(142,14,1)
   !$OMP WORKSHARE
   Tann=sum(Tmois,dim=3)/12.