Custom Query (96 matches)

Sur Jean Zay, avec les options de debug le bench plante (en séquentiel) pour la r3572 dans pbl_surface lors du calcul de tke_x(i,k,is_ave):

       IF (iflag_split .eq.0) THEN
        wake_dltke(:,:,nsrf) = 0.
        DO k = 1, klev+1
           DO j = 1, knon
              i = ni(j)
!jyg<
!!              tke(i,k,nsrf)    = ytke(j,k)
!!              tke(i,k,is_ave) = tke(i,k,is_ave) + ytke(j,k)*ypct(j)
              tke_x(i,k,nsrf)    = ytke(j,k)
              tke_x(i,k,is_ave) = tke_x(i,k,is_ave) + ytke(j,k)*ypct(j)
!>jyg
           ENDDO
        ENDDO

pour i=14, k=1, car ytke(j,k)=nan (probablement pas/mal initialisé)

To better match with other models (ORCHIDEE, NEMO, ...), it would make sense to have the arch.path follow the same rules, i.e. definitions using 3 sets of keywords :

SOMETHING_LIBDIR ="-L....."
SOMETHING_LIB="-l..."
SOMETHING_INCDIR="-I...."

Of course thins implies adapting the makelmdz[_fcm] scripts accordingly.

An important benefit of having matching arch.path is that the file can then be shared by all models when compiling a combination (e.g. LMDZ+DYNAMICO+XIOS+IOIPSL)

Error reported by A. Cozic (02/12/2020):

Lors que je veux tourner le modèle LMDZ couplé avec INCA (158 traceurs) en 144x142x79 je semble dépasser une limite de taille autorisée par netcdf pour les fichiers en double précision. Lors de la création du fichier restart par dynredem0 j'ai le message d'erreur suivant

NetCDF: One or more variable sizes violate format constraints

Je pense cela vient bien du nombre de traceurs et des dimensions du tableau de concentration, car
* le même nombre de traceurs en 144x142x39 ne pose pas de problème
* idem cela passe bien si je travaille avec moins de traceurs en 144x142x79
* si je coupe le fichier restart en deux en mettant les concentrations de traceurs dans un fichier à part (restart_trac.nc) alors tout se passe sans soucis. 

Problem most likeky related to the use of:

NF90_CREATE(filename,NF90_CLOBBER,nid)

in the dynamics, whereas to enable larger files ("64bit_offset" mode),

NF90_CREATE(filename,IOR(NF90_CLOBBER,NF90_64BIT_OFFSET),nid)

would be necessary (and is moreover already used in the physics).