Changeset 905

Timestamp:

Mar 18, 2013, 11:58:34 AM (12 years ago)

Author:

emillour

Message:

Generic GCM:

• removed the "-static" ifort compilation option (problematic on Gnome cluster) and added some additional debug options
• fixed bug in newstart about initialization of albedo and thermal inertia.

EM

Location:

trunk/LMDZ.GENERIC

Files:

• README (modified) (1 diff)
• libf/dyn3d/newstart.F (modified) (2 diffs)
• makegcm_ifort (modified) (2 diffs)

trunk/LMDZ.GENERIC/README

@@ -915 +915 @@
 - corrected a bug in bilinearbig (found with debug mode)
   not harmful but possible -1 index at first call
+
+== 18/03/2013 == EM
+- removed the "-static" ifort compilation option (problematic onGnome cluster)
+  and added some additional debug options
+- fixed bug in newstart about initialization of albedo and thermal inertia.
+
+

trunk/LMDZ.GENERIC/libf/dyn3d/newstart.F

@@ -328 +328 @@
      .        cloudfrac,totalfrac,hice)
 
-        ! copy albedo and soil thermal inertia
+        ! copy albedo and soil thermal inertia on (local) physics grid
         do i=1,ngridmx
           albfi(i) = albedodat(i)
@@ -335 +335 @@
           enddo
         ! build a surfithfi(:) using 1st layer of ithfi(:), which might
-        ! be neede later on if reinitializing soil thermal inertia
+        ! be needed later on if reinitializing soil thermal inertia
           surfithfi(i)=ithfi(i,1)
         enddo
-
+        ! also copy albedo and soil thermal inertia on (local) dynamics grid
+        ! so that options below can manipulate either (but must then ensure
+        ! to correctly recast things on physics grid)
+        call gr_fi_dyn(1,ngridmx,iip1,jjp1,albfi,alb)
+        call gr_fi_dyn(nsoilmx,ngridmx,iip1,jjp1,ithfi,ith)
+        call gr_fi_dyn(1,ngridmx,iip1,jjp1,surfithfi,surfith)
       
       endif

trunk/LMDZ.GENERIC/makegcm_ifort

@@ -188 +188 @@
 else if $LINUX then
    #NB: on gnome -O3 ==> NaNs ...
-   set optim=" -O2 -fp-model precise -ip -mkl=sequential -align all -static "
-   set optim90=" -O2 -fp-model precise -ip -mkl=sequential -align all -static "
-   set optimtru90=" -O2 -fp-model precise -ip -mkl=sequential -align all -static "
+   set optim=" -O2 -fp-model precise -ip -mkl=sequential -align all "
+   set optim90=" -O2 -fp-model precise -ip -mkl=sequential -align all "
+   set optimtru90=" -O2 -fp-model precise -ip -mkl=sequential -align all "
    if ( `hostname` == ciclad1.ipsl.jussieu.fr ) then
        echo "YOU ARE ON CICLAD CLUSTER"
@@ -390 +390 @@
            set optim90="$optim90"" -G1 "
         else if $LINUX then
-           set optim=" -g -no-ftz -traceback -ftrapuv -fp-stack-check -check"
-           set optim90=" -g -no-ftz -traceback -ftrapuv -fp-stack-check -check"
-           set optimtru90=" -g -no-ftz -traceback -ftrapuv -fp-stack-check -check"
+           set optim=" -g -fpe-all=0 -no-ftz -traceback -ftrapuv -fp-stack-check -check all -debug"
+           set optim90=" -g -fpe-all=0 -no-ftz -traceback -ftrapuv -fp-stack-check -check all -debug"
+           set optimtru90=" -g -fpe-all=0 -no-ftz -traceback -ftrapuv -fp-stack-check -check all -debug"
         else 
            echo "pas d option debug predefinie pour cette machine"