Changeset 237 for trunk/SOURCES/Ant16_files/climat_InitMIP_years_perturb_mod.f90

Timestamp:

01/09/19 17:09:26 (6 years ago)

Author:

aquiquet

Message:

Sealevel is now treated as a 2D variable (sealevel_2d while sealevel remains the eustatic sea level), results should remain identical as sealevel_2d is equal to sealevel in this revision.

File:

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

trunk/SOURCES/Ant16_files/climat_InitMIP_years_perturb_mod.f90

@@ -16 +16 @@
 
 
-use module3d_phy,only: nx,ny,S,S0,H,Bsoc,acc,abl,BM,Tann,Tjuly,Ts,time,dt,num_param,num_rep_42,num_forc,dirforcage,dirnameinp,tafor,time,sealevel,coefbmshelf
+use module3d_phy,only: nx,ny,S,S0,H,Bsoc,acc,abl,BM,Tann,Tjuly,Ts,time,dt,num_param,num_rep_42,num_forc,dirforcage,dirnameinp,tafor,time,sealevel,sealevel_2d,coefbmshelf
 use netcdf
 use io_netcdf_grisli
@@ -285 +285 @@
 100 continue
 
+  sealevel_2d(:,:) = sealevel
+
   Tann (:,:) = Ta0 (:,:) + T_lapse_rate * (S(:,:)-S0(:,:)) +Tafor
   Ts(:,:)    = Tann(:,:)