Custom Query (96 matches)

Bug remonté des tests avec dynamico-aquaplanète de Yann. A vérifier, mais il semble qu'il manque dans physiq_mod.F90 les initialisations de mass_solu_aero() et mass_solu_aero_pi()...

Proposition de la modification:

Ligne 3301+ :

ELSE

tausum_aero(:,:,:) = 0.

+ mass_solu_aero(:,:)=0 + mass_solu_aero_pi(:,:)=0

IF (iflag_rrtm .EQ. 0) THEN !--old radiation

tau_aero(:,:,:,:) = 1.e-15 piz_aero(:,:,:,:) = 1. cg_aero(:,:,:,:) = 0.

ELSE

tau_aero_sw_rrtm(:,:,:,:) = 1.e-15 tau_aero_lw_rrtm(:,:,:,:) = 1.e-15 piz_aero_sw_rrtm(:,:,:,:) = 1.0 cg_aero_sw_rrtm(:,:,:,:) = 0.0

ENDIF

ENDIF

Rapport de possible bug dans yamada4 identifié par Yann dans le prototype dynamico-aquaplanète:

Produit d'un svn diff dans corrections dans yamada4.F90:

@@ -246,7 +246,7 @@
    IF (iflag_pbl==6) THEN

      DO k = 2, klev
-      q2(:, k) = l(:, k)**2*zz(:, k)
+      q2(:ngrid, k) = l(:ngrid, k)**2*zz(:ngrid, k)
      END DO


@@ -339,12 +339,14 @@
      ! print*,'Schema mixte D'
      ! print*,'Longueur ',l(:,:)
      DO k = 2, klev - 1
-      l(:, k) = max(l(:,k), 1.)
-      km(:, k) = l(:, k)*sqrt(q2(:,k))*sm(:, k)
-      q2(:, k) = q2(:, k) + dt*km(:, k)*m2(:, k)*(1.-rif(:,k))
-      q2(:, k) = min(max(q2(:,k),1.E-10), 1.E4)
-      q2(:, k) = 1./(1./sqrt(q2(:,k))+dt/(2*l(:,k)*b1))
-      q2(:, k) = q2(:, k)*q2(:, k)
+      DO ig = 1, ngrid
+        l(ig, k) = max(l(ig,k), 1.)
+        km(ig, k) = l(ig, k)*sqrt(q2(ig,k))*sm(ig, k)
+        q2(ig, k) = q2(ig, k) + dt*km(ig, k)*m2(ig, k)*(1.-rif(ig,k))
+        q2(ig, k) = min(max(q2(ig,k),1.E-10), 1.E4)
+        q2(ig, k) = 1./(1./sqrt(q2(ig,k))+dt/(2*l(ig,k)*b1))
+        q2(ig, k) = q2(ig, k)*q2(ig, k)
+      END DO
      END DO


@@ -376,7 +378,7 @@
    ! Transport diffusif vertical de la TKE.
    IF (iflag_pbl>=12) THEN
      ! print*,'YAMADA VDIF'
-    q2(:, 1) = q2(:, 2)
+    q2(1:ngrid, 1) = q2(1:ngrid, 2)
      CALL vdif_q2(dt, g, rconst, ngrid, plev, temp, kq, q2)
    END IF

With XIOS the model crashes on Ada with XIOS (but not without), with the following error message:

terminate called after throwing an instance of 'xios::CNetCdfException'
  what():  Error when calling function nc_enddef(ncId)
NetCDF: HDF error
Unable to end define mode of this file, given its id: 65536
forrtl: error (76): Abort trap signal
Image              PC                Routine            Line        Source             
libifcoremt.so.5   00002ADA7C66B3EF  for__signal_handl     Unknown  Unknown
libpthread-2.12.s  00002ADA7A7F27E0  Unknown               Unknown  Unknown
libc-2.12.so       00002ADA7EEED495  gsignal               Unknown  Unknown
libc-2.12.so       00002ADA7EEEEC75  abort                 Unknown  Unknown
libstdc++.so.6.0.  00002ADA7AA90BCD  _ZN9__gnu_cxx27__     Unknown  Unknown
libstdc++.so.6.0.  00002ADA7AA8EB86  Unknown               Unknown  Unknown
libstdc++.so.6.0.  00002ADA7AA8EBD1  Unknown               Unknown  Unknown
libstdc++.so.6.0.  00002ADA7AA8EDE8  Unknown               Unknown  Unknown
...
gcm_48x36x19_phyl  0000000000D33951  iophy_mp_histwrit         974  iophy.f90
gcm_48x36x19_phyl  0000000002EE3365  phys_output_write         458  phys_output_write_mod.f90
gcm_48x36x19_phyl  00000000018A0037  physiq_mod_mp_phy        4581  physiq_mod.f90
gcm_48x36x19_phyl  00000000013D9905  callphysiq_mod_mp          70  callphysiq_mod.f90
gcm_48x36x19_phyl  0000000001397F34  calfis_loc_               718  calfis_loc.f
gcm_48x36x19_phyl  00000000009B2B81  call_calfis_mod_m         205  call_calfis_mod.f90
gcm_48x36x19_phyl  0000000000596751  leapfrog_loc_             796  leapfrog_loc.f
gcm_48x36x19_phyl  000000000042AC4B  MAIN__                    444  gcm.f90

Problem found on Ada using LMDZ r3543 (but also occurs for other recent revisions) and XIO revision 1659.