Changeset 290

Timestamp:

11/18/19 17:08:48 (4 years ago)

Author:

dumas

Message:

ISMIP forcing : grad_smb is not interpolated between years

File:

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

trunk/SOURCES/Ant16_files/climat_InitMIP_years_perturb_mod.f90

@@ -90 +90 @@
 
   namelist/clim_smb_T_gen/smb_file,coef_smb_unit,temp_annual_file
+! version Aurelien Groenland avec file_lapse :
+!  namelist/smb_anom_initMIP/file_smb_anom,file_lapse,coef_smb_anom_unit,nb_snap,time_depart_snaps,massb_time
+! version utilisee pour Antarctique :  
   namelist/smb_anom_initMIP/file_smb_anom,file_lapse,coef_smb_anom_unit,nb_snap,time_depart_snaps,massb_time
   
@@ -353 +356 @@
      bm_time(:,:) = bm_0(:,:)
      tann_time(:,:) = ta0(:,:)
+     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 : 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)
+        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 
@@ -375 +382 @@
            tann_time(:,:) = ta0(:,:) + tann_anom(:,:,k)
            if (massb_time == 3 ) then
-              grad_smb_time(:,:) =   grad_smb(:,:,k) + (grad_smb(:,:,k+1)-grad_smb(:,:,k)) *   &
-                (time-time_snap(k))/(time_snap(k+1)-time_snap(k))
-              grad_tann_time(:,:) = grad_tann(:,:,k) + (grad_tann(:,:,k+1)-grad_tann(:,:,k)) *   &
-                (time-time_snap(k))/(time_snap(k+1)-time_snap(k))
+              grad_smb_time(:,:) =   grad_smb(:,:,k)
+              grad_tann_time(:,:) = grad_tann(:,:,k)
            endif
            exit