Changeset 52 for trunk/libf/phylmd

Timestamp:

Jan 31, 2011, 12:42:42 AM (14 years ago)

Author:

aslmd

Message:

chantier principal du commit
--- version LMDZ5 qui fonctionne pour tests geantes
--- prochaine etape, tests sur GNOME

M libf/dyn3dpar/comconst.h
M libf/dyn3dpar/conf_planete.F90
ajout du flux de chaleur intrinseque: ihf
[par defaut il est nul]

M libf/dyn3dpar/gcm.F
changements cosmetiques
[pour diff plus efficace avec version non par]

M libf/dyn3dpar/iniacademic.F
possibilites de variations latitudinales
de temperature plus originales
[seulement pour planet_type.eq."giant"]

M libf/dyn3dpar/leapfrog_p.F

1. ajout d'une tendance causee par le flux de chaleur intrinseque

(seulement prise en compte si planet_type.eq."giant")

2. correction bugs problematiques a la compilation et au run

--> probleme dans les boucles (l'indice etait llm et non l)
--> ajout de SAVE pour les variables paralleles
--> correction des declarations de variables manquantes

M libf/dyn3dpar/calfis_p.F
correction d'une deuxieme parenthese manquante sur ALLOCATE(zteta(klon,llm))

M libf/phylmd/regr_lat_time_climoz_m.F90
erreur a la compilation avec FCM... il s'agit d'une routine terrestre
il y a visiblement un probleme avec o3_in
en attendant, les lignes sont commentees avec !AS

A deftanks/giant 8 fichiers
ajout de fichiers de configuration typiques pour les geantes gazeuses
[experimental pour le moment... on est loin de jupiter]

--> comparaisons entre un run ancien [avec LMDZ5-dev sur SVN ipsl sans cp var]
et run avec version sur ce SVN planeto donne des resultats similaires

pratique

A ioipsl
A ioipsl/compile_ioipsl.bash
A ioipsl/util 16 fichiers
script et utilitaire pour compiler IOIPSL de facon independante
il suffit d'executer ./compile_ioipsl.bash

M arch/arch-AMD64_CICLAD.path
si IOIPSL a ete compile avec la methode precedente, les bons
PATH sont definis dans ce fichier [le NETCDF est aussi OK]

M 000-README-svn
mise a jour options "svn status"

M mars/libf/phymars/meso_callkeys.h
mise a jour mineure du fichier
[ecri_phys etait defini mais pas dans la liste]

File:

• trunk/libf/phylmd/regr_lat_time_climoz_m.F90 (modified) (5 diffs)

trunk/libf/phylmd/regr_lat_time_climoz_m.F90

@@ -239 +239 @@
     end if
 
-    if (desc_lat) o3_in = o3_in(n_lat:1:-1, :, :, :)
-    if (desc_plev) o3_in = o3_in(:, n_plev:1:-1, :, :)
+    !!!! Aymeric; problem with compilation here.... pb with o3_in
+    !AS if (desc_lat) o3_in = o3_in(n_lat:1:-1, :, :, :)
+    !AS if (desc_plev) o3_in = o3_in(:, n_plev:1:-1, :, :)
 
     do m = 1, read_climoz
@@ -251 +252 @@
              ! Next to the south pole:
              j = 1
-             do while (o3_in(j, 1, l, m) == missing_value)
-                j = j + 1
-             end do
-             if (j > 1) o3_in(:j-1, :, l, m) = &
-                  spread(o3_in(j, :, l, m), dim=1, ncopies=j-1)
+!AS             do while (o3_in(j, 1, l, m) == missing_value)
+!AS                j = j + 1
+!AS             end do
+!AS             if (j > 1) o3_in(:j-1, :, l, m) = &
+!AS                  spread(o3_in(j, :, l, m), dim=1, ncopies=j-1)
              
              ! Next to the north pole:
              j = n_lat
-             do while (o3_in(j, 1, l, m) == missing_value)
-                j = j - 1
-             end do
-             if (j < n_lat) o3_in(j+1:, :, l, m) = &
-                  spread(o3_in(j, :, l, m), dim=1, ncopies=n_lat-j)
+!AS             do while (o3_in(j, 1, l, m) == missing_value)
+!AS                j = j - 1
+!AS             end do
+!AS             if (j < n_lat) o3_in(j+1:, :, l, m) = &
+!AS                  spread(o3_in(j, :, l, m), dim=1, ncopies=n_lat-j)
 
              ! Take care of missing values at high pressure:
@@ -272 +273 @@
                 ! missing values so the highest level has a valid value.
                 k = 2
-                do while  (o3_in(j, k, l, m) /= missing_value .and. k < n_plev)
-                   k = k + 1
-                end do
+!AS                do while  (o3_in(j, k, l, m) /= missing_value .and. k < n_plev)
+!AS                   k = k + 1
+!AS                end do
                 ! Replace missing values with the valid value at the
                 ! lowest level above missing values:
-                if (o3_in(j, k, l, m) == missing_value) &
-                     o3_in(j, k:n_plev, l, m) = o3_in(j, k-1, l, m)
+!AS                if (o3_in(j, k, l, m) == missing_value) &
+!AS                     o3_in(j, k:n_plev, l, m) = o3_in(j, k-1, l, m)
              end do
           end do
@@ -298 +299 @@
        print *, &
             "Found 12 months in ozone climatologies, assuming periodicity..."
-       o3_regr_lat(jjm+1:1:-1, :, 1:12, :) = regr1_step_av(o3_in, &
-            xs=sin(lat_in_edg), xt=sin((/- pi / 2, rlatv(jjm:1:-1), pi / 2/)))
+!AS       o3_regr_lat(jjm+1:1:-1, :, 1:12, :) = regr1_step_av(o3_in, &
+!AS            xs=sin(lat_in_edg), xt=sin((/- pi / 2, rlatv(jjm:1:-1), pi / 2/)))
        ! (invert order of indices in "o3_regr_lat" because "rlatu" is
        ! in descending order)
@@ -309 +310 @@
     else
        print *, "Using 14 months in ozone climatologies..."
-       o3_regr_lat(jjm+1:1:-1, :, :, :) = regr1_step_av(o3_in, &
-            xs=sin(lat_in_edg), xt=sin((/- pi / 2, rlatv(jjm:1:-1), pi / 2/)))
+!AS       o3_regr_lat(jjm+1:1:-1, :, :, :) = regr1_step_av(o3_in, &
+!AS            xs=sin(lat_in_edg), xt=sin((/- pi / 2, rlatv(jjm:1:-1), pi / 2/)))
        ! (invert order of indices in "o3_regr_lat" because "rlatu" is
        ! in descending order)