Project

General

Profile

Feature #4981

Photolysis at polar caps

Added by Twan van Noije almost 2 years ago. Updated about 1 year ago.

Status:
Closed
Priority:
Low
Assignee:
-
Category:
-
Target version:
Start date:
10/26/2017
Due date:
% Done:

100%


Description

In the current TM5 code, the photolysis rates are set to zero at solar zenith angles higher than 85 degrees, because the parameterization cannot handle this region. It recommended to apply some sort of interpolation. Probably a linear interpolation in cos(SZA), reaching zero at either 90 degrees or 92-93 will do.

This is currently implemented and tested in C-IFS. The TM5 fix can follow the same approach.

no2_L29_18ut_old_photolysis.png View (370 KB) Philippe Le Sager, 07/20/2018 03:13 PM

no2_L30_18ut_new_photolysis.png View (383 KB) Philippe Le Sager, 07/20/2018 03:13 PM

History

#1 Updated by Twan van Noije almost 2 years ago

  • Description updated (diff)

#2 Updated by Philippe Le Sager about 1 year ago

  • Tracker changed from Bug to Feature
  • Status changed from New to In Progress
  • % Done changed from 0 to 20

H. Eskes wrote:

Vincent gave me a simple patch of the photolysis, existing of a linear extrapolation from sza=85 to sza=94. This provides a reasonable match with the photolysis in the model of Prather. I made a test run of one day, and it makesquite an impact near sunrise, and looks more realistic. See images (stratospheric layer).

#3 Updated by Philippe Le Sager about 1 year ago

  • Target version set to TM5-MP 2.1

#4 Updated by Philippe Le Sager about 1 year ago

Here are the plots from Henk. Current code:

and new code:

#5 Updated by Philippe Le Sager about 1 year ago

Created a dedicated branch to test the code: photo94sza. First test run with run time checks crashed in pdump with an overflow. Traceback"

forrtl: error (72): floating overflow
Image              PC                Routine            Line        Source             
libintlc.so.5      00002AAAAC30ACF1  tbk_trace_stack_i     Unknown  Unknown
libintlc.so.5      00002AAAAC308E2B  tbk_string_stack_     Unknown  Unknown
tm5mp.x            00000000015C4E64  Unknown               Unknown  Unknown
tm5mp.x            00000000015C4C76  Unknown               Unknown  Unknown
tm5mp.x            000000000156C879  Unknown               Unknown  Unknown
tm5mp.x            000000000157129C  Unknown               Unknown  Unknown
libpthread-2.11.3  00002AAAAD8A1850  Unknown               Unknown  Unknown
tm5mp.x            0000000000CD2FE4  mdf_mp_mdf_put_va       22829  mdf.F90
tm5mp.x            0000000001462232  user_output_pdump        4107  user_output_pdump.F90
tm5mp.x            00000000014906DE  user_output_pdump        1122  user_output_pdump.F90
tm5mp.x            0000000001302415  user_output_mp_us         968  user_output.F90
tm5mp.x            0000000000E6AB0B  modelintegration_         821  modelIntegration.F90
tm5mp.x            0000000000E68F95  modelintegration_         540  modelIntegration.F90
tm5mp.x            0000000001016452  tm5_mp_tm5_model_        1204  tm5.F90
tm5mp.x            00000000010C6616  MAIN__                     63  tracer.F90
tm5mp.x            000000000041021E  Unknown               Unknown  Unknown
libc-2.11.3.so     00002AAAADD47C36  __libc_start_main     Unknown  Unknown
tm5mp.x            0000000000410129  Unknown               Unknown  Unknown
where l. 4107 of pdump is the call to MDF_put_var:
      ! temperature (3d)
      call MDF_Put_Var( RF%ncid, RF%varid_temp, RF%data3d_t(:,:,:,:), status, start=(/i1,j1,1,1/), count=(/imr,jmr,lmr,RF%n_rec/) )
      IF_NOTOK_MDF(fid=RF%ncid)      

#6 Updated by Philippe Le Sager about 1 year ago

  • Status changed from In Progress to Resolved
  • % Done changed from 20 to 100

PDUMP overflow bug is not related to the new photolysis (discussed in #11761). The branch has been reintegrated in the trunk in r911.

#7 Updated by Philippe Le Sager about 1 year ago

  • Status changed from Resolved to Closed

Also available in: Atom PDF