correct the attributes on p2l output variable (#1395)

Author: mickaelaccensi

model/src/w3ounfmetamd.F90

@@ -3639,14 +3639,14 @@ SUBROUTINE DEFAULT_META()
     ! IFI=6, IFJ=9, P2L
     META => GROUP(6)%FIELD(9)%META
     ! Information for spectral microseismic generation data (2nd file)
-    META(1)%FSC = 0.0004
     META(1)%VARNM='p2l'
     META(1)%VARNL='base ten logarithm of power spectral density of equivalent surface pressure'
     !META(1)%VARNS='base_ten_logarithm_of_power_spectral_density_of_equivalent_surface_pressure'
     META(1)%VARNS=''
     META(1)%VARNG='base_ten_logarithm_of_power_spectral_density_of_equivalent_surface_pressure'
-    IF (NCVARTYPE.EQ.2) THEN
+    IF (NCVARTYPE.LE.3) THEN
       META(1)%UNITS='log10(Pa2 m2 s+1E-12)'
+      META(1)%FSC = 0.0004
       META(1)%VMIN = -12.
       META(1)%VMAX = 12.
     ELSE

model/src/w3profsmd.F90

@@ -677,7 +677,7 @@ SUBROUTINE W3XYPFSN2 ( ISP, C, LCALC, RD10, RD20, DT, AC)
         U(IP) = MAX(0.d0,U(IP)-DTSI(IP)*ST(IP)*(1-IOBPA(IP)))*DBLE(IOBPD(ITH,IP))
 
 #ifdef W3_REF1
-        WRITE(10111,*) REFPARS(3), IOBPD(ITH,IP), IOBPA(IP), U(IP),  AC(IP)
+        ! WRITE(10111,*) REFPARS(3), IOBPD(ITH,IP), IOBPA(IP), U(IP),  AC(IP)
         IF (REFPARS(3).LT.0.5.AND.IOBPD(ITH,IP).EQ.0.AND.IOBPA(IP).EQ.0) THEN
           U(IP) = AC(IP) ! restores reflected boundary values
         ENDIF