Changeset 237 for trunk/SOURCES/Hudson_files/climat-hudson_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/Hudson_files/climat-hudson_mod.f90 (modified) (2 diffs)

trunk/SOURCES/Hudson_files/climat-hudson_mod.f90

@@ -102 +102 @@
 !jalv: test pour que le bilan de masse soit 0 sur l'ocean plus profond de 600 m:
 !      c'est le bon critere??:
-where (sealevel-BSOC(:,:).gt.600) acc(:,:)=0.
-where (sealevel-BSOC(:,:).gt.600) abl(:,:)=0.
+where (sealevel_2d(:,:)-BSOC(:,:).gt.600) acc(:,:)=0.
+where (sealevel_2d(:,:)-BSOC(:,:).gt.600) abl(:,:)=0.
 
-where ((sealevel-BSOC(:,:).gt.100).and.(dist_heino(:,:).gt.2500.)) abl(:,:)=-1.
-where ((sealevel-BSOC(:,:).gt.100).and.(dist_heino(:,:).gt.2600.)) abl(:,:)=-5.
-where ((sealevel-BSOC(:,:).gt.100).and.(dist_heino(:,:).gt.2650.)) abl(:,:)=-10.
-where ((sealevel-BSOC(:,:).gt.100).and.(dist_heino(:,:).gt.2700.)) abl(:,:)=-20.
+where ((sealevel_2d(:,:)-BSOC(:,:).gt.100).and.(dist_heino(:,:).gt.2500.)) abl(:,:)=-1.
+where ((sealevel_2d(:,:)-BSOC(:,:).gt.100).and.(dist_heino(:,:).gt.2600.)) abl(:,:)=-5.
+where ((sealevel_2d(:,:)-BSOC(:,:).gt.100).and.(dist_heino(:,:).gt.2650.)) abl(:,:)=-10.
+where ((sealevel_2d(:,:)-BSOC(:,:).gt.100).and.(dist_heino(:,:).gt.2700.)) abl(:,:)=-20.
 
 bm(:,:)=acc(:,:)+abl(:,:)
@@ -135 +135 @@
 !      c'est le bon critere??:
 !where (BSOC(:,:).lt.B(:,:)) Tann(:,:)=0.
-where (sealevel-BSOC(:,:).gt.600) Tann(:,:)=0. 
+where (sealevel_2d(:,:)-BSOC(:,:).gt.600) Tann(:,:)=0.