Changeset 237 for trunk/SOURCES/climat-forcage-insolation_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/climat-forcage-insolation_mod.f90 (modified) (2 diffs)

trunk/SOURCES/climat-forcage-insolation_mod.f90

@@ -12 +12 @@
 ! C. Dumas 07/2015
 
-use module3D_phy,only:nx,ny,S,H,flot,slv,Tann,Tjuly,Tmois,acc,coefbmshelf,ro,num_param,num_rep_42,dirnameinp,time,pi,dx,xlong,ylat,bsoc         
+use module3D_phy,only:nx,ny,S,H,flot,Tann,Tjuly,Tmois,acc,coefbmshelf,sealevel_2d,ro,num_param,num_rep_42,dirnameinp,time,pi,dx,xlong,ylat,bsoc         
 use netcdf
 use io_netcdf_grisli
@@ -637 +637 @@
 do j=1,ny
    do i=1,nx
-      Zs(i,j)=max(slv(i,j),S(i,j))
+      Zs(i,j)=max(sealevel_2d(i,j),S(i,j))
       do ictr=CO2SnapMin,CO2SnapMax,1
          do igtr=1,gtr