Changeset 288

Timestamp:

11/15/19 17:32:31 (5 years ago)

Author:

dumas

Message:

ISMIP forcing : smb and Tann without interpolation between years

File:

• trunk/SOURCES/Ant16_files/climat_InitMIP_years_perturb_mod.f90 (modified) (6 diffs)

trunk/SOURCES/Ant16_files/climat_InitMIP_years_perturb_mod.f90

@@ -62 +62 @@
 integer                              :: nb_snap           !> nombre de snapshots
 
+real,dimension(nx,ny)                :: smb_anom_last20years  !> anomalie de smb des 20 dernieres annees du forcage
+real,dimension(nx,ny)                :: tann_anom_last20years !> anomalie de Tann des 20 dernieres annees du forcage
 
 contains
@@ -79 +81 @@
 
   integer            :: err                ! recuperation d'erreur
-  integer            :: i
+  integer            :: i,k
 
   !aurel for Tafor:
@@ -281 +283 @@
          
       endif
-      
+! calcul de la moyenne des 20 dernieres annees du forcage
+    do k=nb_snap,nb_snap-19,-1
+      smb_anom_last20years(:,:) = smb_anom_last20years(:,:) + smb_anom(:,:,k)
+      tann_anom_last20years(:,:) = tann_anom_last20years(:,:) + tann_anom(:,:,k)
+    enddo
+      smb_anom_last20years(:,:) = smb_anom_last20years(:,:)/20.
+      tann_anom_last20years(:,:) = tann_anom_last20years(:,:)/20.
   endif     
 
@@ -331 +339 @@
 
   implicit none
-  integer             :: k, k_snap               ! pour calculer les indices de temps
-  real :: time_prec
+  integer             :: k                ! pour calculer les indices de temps
   real,dimension(nx,ny) :: bm_time, tann_time ! pour calcul Bm et Tann
   real,dimension(nx,ny) :: grad_smb_time, grad_tann_time ! pour calcul Bm et Tann avec correction verticale
@@ -346 +353 @@
      bm_time(:,:) = bm_0(:,:)
      tann_time(:,:) = ta0(:,:)
-     k_snap       = 1
-!     S_ref(:,:)   = S(:,:)                   ! du coup sera la surface de reference avant le forcage
-     time_prec    = time
+  else if (time.ge.time_snap(nb_snap)) then  ! time apres le forcage : constant = derniere année nb_snap
+     bm_time(:,:) =  bm_0(:,:) + smb_anom (:,:,nb_snap)
+     tann_time(:,:) = ta0(:,:) + tann_anom(:,:,nb_snap)
      if (massb_time == 3 ) then
         grad_smb_time(:,:) =   grad_smb (:,:,1)
         grad_tann_time(:,:) = grad_tann(:,:,1)
      endif
-  else if (time.ge.time_snap(nb_snap)) then  ! time apres le forcage
-     bm_time(:,:) =  bm_0(:,:) + smb_anom (:,:,nb_snap)
-     tann_time(:,:) = ta0(:,:) + tann_anom(:,:,nb_snap)
-     if (abs(time-time_prec-1.).lt.dt) then   ! 
-        time_prec = time_prec + 1
-     endif
-     if (massb_time == 3 ) then
-        grad_smb_time(:,:) =   grad_smb (:,:,nb_snap)
-        grad_tann_time(:,:) = grad_tann(:,:,nb_snap)
-     endif
+!~   else if (time.ge.time_snap(nb_snap)) then  ! time apres le forcage : constant moy 20 dernieres années 
+!~      bm_time(:,:) =  bm_0(:,:) + smb_anom_last20years(:,:)
+!~      tann_time(:,:) = ta0(:,:) + tann_anom_last20years(:,:)
   else                                       ! cas general
      do k = 1 , nb_snap-1
         if((time.ge.time_snap(k)).and.(time.lt.time_snap(k+1))) then ! entre k et k+1 
-           bm_time(:,:) = bm_0(:,:) + smb_anom(:,:,k) + (smb_anom(:,:,k+1)-smb_anom(:,:,k)) *   &
-                (time-time_snap(k))/(time_snap(k+1)-time_snap(k))
-           tann_time(:,:) = ta0(:,:) + tann_anom(:,:,k) + (tann_anom(:,:,k+1)-tann_anom(:,:,k)) *   &
-                (time-time_snap(k))/(time_snap(k+1)-time_snap(k))
-! exactement sur le snapshot et avec un ecart 1 an par rapport au precedent stockage
-!           write(6,*) 'time,tests',k,time,time-time_snap(k),time-time_prec-1.
-           if ((abs(time-time_snap(k)).le.dt).and.(abs(time-time_prec-1.).lt.dt)) then    
-              k_snap    = k
-              time_prec = time_snap(k)     ! time_prec est le temps du precedent
-           endif
+!cdc version avec interpolation lineaire entre 2 années        
+!~            bm_time(:,:) = bm_0(:,:) + smb_anom(:,:,k) + (smb_anom(:,:,k+1)-smb_anom(:,:,k)) *   &
+!~                 (time-time_snap(k))/(time_snap(k+1)-time_snap(k))
+!~            tann_time(:,:) = ta0(:,:) + tann_anom(:,:,k) + (tann_anom(:,:,k+1)-tann_anom(:,:,k)) *   &
+!~                 (time-time_snap(k))/(time_snap(k+1)-time_snap(k))
+!cdc version ISMIP : forcage annee n pendant un an sans interpolation.
+           bm_time(:,:) = bm_0(:,:) + smb_anom(:,:,k)
+           tann_time(:,:) = ta0(:,:) + tann_anom(:,:,k)
            if (massb_time == 3 ) then
               grad_smb_time(:,:) =   grad_smb(:,:,k) + (grad_smb(:,:,k+1)-grad_smb(:,:,k)) *   &
@@ -396 +394 @@
   endif
 
-!~ ! garde les 10 dernieres annees et calcule la moyenne
-!~   if (icum.eq.1) then                           ! stockage dans le tableau bm_ref_10
-!~      do k = 9,1,-1
-!~         bm_ref_10(:,:,k+1) = bm_ref_10(:,:,k)   ! on decale tous les elements
-!~      end do
-!~      bm_ref_10(:,:,k+1) = bm(:,:)               ! le plus recent est en position 1
-!~      i_moy              = i_moy +1              ! compte combien il y en a pour la moyenne
-!~      i_moy              = min(i_moy,10)
-!~      bm_ref_t(:,:)   = 0.
-!~      do k = 1,i_moy
-!~         bm_ref_t(:,:) = bm_ref_t(:,:) + bm_ref_10(:,:,k)
-!~      end do
-!~      bm_ref_t(:,:)    = bm_ref_t(:,:)/i_moy
-!~      write(6,898) time, time_prec, icum, i_moy
-!~ 898  format('cumul pour gradient  ',2(f0.3,1x),2(i0,1x))
-!~   end if
 end subroutine massb_ISMIP_RCM