From: Solmon Fabien <fsolmon@ictp.it>
To: R Climate <rclimate@yahoo.com>
Cc: "regcnet@lists.ictp.it" <regcnet@lists.ictp.it>
Sent: Saturday, January 12, 2013 10:35:16 AM
Subject: Re: [RegCNET] Dust Run
Hello !
I think you need to rerun terrain with the option ltexture = true. After
that you can check that you have the texture fields ( needed by dust
emission) in your DOMAIN000.nc file.
Also you need to perform the chbc step to generate dust boundary
conditions if you want them, or set ichebdy=0 if you want to relax toward
zero concentration
at the boundary ( in which case the chbc step is not
necessary).
Finally the rdtsemfac might also be adjusted.
Usually I use rdstemfac=1 (default) when I use the new option
ichdustemd = 2 ( which I think is better at least over the sahara). With
the option ichdustemd =1 , the emission more intense and I remember using
a rdstemfac = 0.7 to get more realistic emissions. But again all of this
domain dependant ..
Fabien
> Hello Dr. Giulian and Regcnet,
> I'm running a dust simulation with the RegCM4.3.5.1 and I put the
> following selection in the name-list:
>
> ltexture = .true.,
>
> chemsimtype = 'DUST ',
> ichsolver = 1,
> ichsursrc = 1,
> ichdrdepo = 1,
> ichebdy = 1,
> ichcumtra = 1,
> ichremlsc = 1,
> ichremcvc = 1,
>
ichdustemd = 1,
> ichdiag = 0,
> idirect = 1,
> rdstemfac = 1.0,
>
>
> Then, when I run the simulation I got the following:
>
> Constant hor. diff. coef. = 0.360000E+05 m^2 s-1
> Maximumt hor. diff. coef. = 0.375000E+06 m^2 s-1
> Calculate solar declination angle at 2009010100
> JDay 0.00 solar declination angle = -23.05862917 degrees
> solar TSI irradiance = 1367.0000 W/m^2
> READY IC DATA for 2009-01-01 00:00:00 UTC
> SEARCH BC data for 2009010106
> READY BC
from 2009010100 to 2009010106
> # Ozone profiles #
> 0.496E-07 0.496E-07 0.496E-07 0.496E-07 0.496E-07 0.519E-07
> 0.527E-07
> 0.574E-07 0.623E-07 0.710E-07 0.846E-07 0.106E-06 0.148E-06
> 0.223E-06
> 0.223E-06 0.409E-06 0.793E-06 0.141E-05 0.317E-05 0.00
> Calling inidust
> No such file or
> directory
>
--------------------------------------------------------------------------
> MPI_ABORT was invoked on rank 7 in communicator MPI_COMM_WORLD
> with errorcode 1.
>
> NOTE: invoking MPI_ABORT causes Open MPI to kill all MPI processes.
> You may or may not see output from other processes, depending on
> exactly when Open MPI kills them.
>
>
> Any help is welcome
>
>
> Thanks
> Rahul
> Manger of Climate Division_______________________________________________
> RegCNET mailing list
>
RegCNET@lists.ictp.it>
https://lists.ictp.it/mailman/listinfo.cgi/regcnet--
+39 040 2240206 (office) | Earth System Physics
+39 040 2240449 (fax) | ICTP
(UNESCO-IAEA)
esp.ictp.it | Strada Costiera 11
| 34014 Trieste, Italy