Index: trunk/MESOSCALE/LMD_MM_MARS/SRC/POSTPROC/api.F90 =================================================================== --- trunk/MESOSCALE/LMD_MM_MARS/SRC/POSTPROC/api.F90 (revision 179) +++ trunk/MESOSCALE/LMD_MM_MARS/SRC/POSTPROC/api.F90 (revision 186) @@ -35,4 +35,68 @@ INCLUDE 'netcdf.inc' + ! + ! VARIABLES + ! + CHARACTER (LEN=500) :: path_to_input + CHARACTER (LEN=500) :: path_to_output + CHARACTER (LEN=500) :: input_name + CHARACTER (LEN=500) :: output_name + CHARACTER (LEN=20) :: process + CHARACTER (LEN=2000) :: fields + REAL, DIMENSION(299) :: interp_levels + INTEGER :: interp_method=1 + INTEGER :: extrapolate=0 + LOGICAL :: debug=.FALSE. + LOGICAL :: unstagger_grid=.FALSE. + LOGICAL :: bit64=.FALSE. + LOGICAL :: oldvar=.TRUE. + + INTEGER :: funit,ios + LOGICAL :: is_used + + ! + ! NAMELISTS + ! + NAMELIST /io/ path_to_input, input_name, path_to_output, output_name, & + process, fields, debug, bit64, oldvar + NAMELIST /interp_in/ interp_levels, interp_method, extrapolate, unstagger_grid + + ! + ! DEFAULT VALUES for VARIABLES + ! + path_to_input = './' + path_to_output = './' + output_name = ' ' + interp_levels = -99999. + process = 'all' + + + ! + ! READ NAMELIST + ! + DO funit=10,100 + INQUIRE(unit=funit, opened=is_used) + IF (.not. is_used) EXIT + END DO + OPEN(funit,file='namelist.api',status='old',form='formatted',iostat=ios) + IF ( ios /= 0 ) STOP "ERROR opening namelist.api" + READ(funit,io) + READ(funit,interp_in) + CLOSE(funit) + + !!! MAIN CALL + CALL api_main ( path_to_input, input_name, path_to_output, output_name, & + process, fields, debug, bit64, oldvar, & + interp_levels, interp_method, extrapolate, unstagger_grid, -99999. ) + + END PROGRAM api + + SUBROUTINE api_main ( path_to_input, input_name, path_to_output, output_name, & + process, fields, debug, bit64, oldvar, & + interp_levels, interp_method, extrapolate, unstagger_grid, onelevel ) + + IMPLICIT NONE + INCLUDE 'netcdf.inc' + !! !! EARTH CONSTANTS @@ -62,6 +126,6 @@ CHARACTER, ALLOCATABLE, DIMENSION(:,:,:,:) :: text CHARACTER (LEN=31),ALLOCATABLE, DIMENSION(:) :: dnamei, dnamej - CHARACTER(LEN=250),ALLOCATABLE, DIMENSION(:) :: input_file_names - CHARACTER(LEN=250),ALLOCATABLE, DIMENSION(:) :: output_file_names + CHARACTER(LEN=500),ALLOCATABLE, DIMENSION(:) :: input_file_names + CHARACTER(LEN=500),ALLOCATABLE, DIMENSION(:) :: output_file_names DOUBLE PRECISION, ALLOCATABLE, DIMENSION(:,:,:,:) :: ddata1, ddata2 REAL, ALLOCATABLE, DIMENSION(:,:,:,:) :: data1, data2, data3 @@ -79,16 +143,16 @@ INTEGER, DIMENSION(4) :: dims_in, dims_out INTEGER, DIMENSION(6) :: ishape, jshape - CHARACTER (LEN=80) :: cval + CHARACTER (LEN=500) :: cval !80 CHARACTER (LEN=31) :: cname, test_dim_name - CHARACTER (LEN=80) :: input_file, output_file, att_text - CHARACTER (LEN=250) :: path_to_input - CHARACTER (LEN=250) :: path_to_output - CHARACTER (LEN=250) :: input_name - CHARACTER (LEN=250) :: output_name, tmp_name + CHARACTER (LEN=500) :: input_file, output_file, att_text !80 + CHARACTER (LEN=500) :: path_to_input + CHARACTER (LEN=500) :: path_to_output + CHARACTER (LEN=500) :: input_name + CHARACTER (LEN=500) :: output_name, tmp_name CHARACTER (LEN=10) :: option CHARACTER (LEN=132) :: command CHARACTER (LEN=20) :: process, dummy CHARACTER (LEN=2000) :: fields, process_these_fields - REAL, DIMENSION(299) :: interp_levels + REAL, DIMENSION(299) :: interp_levels REAL :: rval REAL :: MISSING=1.e36 @@ -96,6 +160,6 @@ INTEGER :: map_proj INTEGER :: LINLOG = 1 - INTEGER :: interp_method=1 - INTEGER :: extrapolate=0 + INTEGER :: interp_method!=1 + INTEGER :: extrapolate!=0 INTEGER :: ncid, mcid, rcode INTEGER :: idm, ndims, nvars, natt, ngatts @@ -111,40 +175,17 @@ INTEGER :: kk LOGICAL :: is_used - LOGICAL :: debug=.FALSE. + LOGICAL :: debug!=.FALSE. LOGICAL :: interpolate=.FALSE. - LOGICAL :: unstagger_grid=.FALSE. + LOGICAL :: unstagger_grid!=.FALSE. LOGICAL :: fix_meta_stag=.FALSE. - LOGICAL :: bit64=.FALSE. + LOGICAL :: bit64!=.FALSE. LOGICAL :: first=.TRUE. - LOGICAL :: oldvar=.FALSE. - - ! - ! NAMELISTS - ! - NAMELIST /io/ path_to_input, input_name, path_to_output, output_name, & - process, fields, debug, bit64, oldvar - NAMELIST /interp_in/ interp_levels, interp_method, extrapolate, unstagger_grid - - ! - ! DEFAULT VALUES for VARIABLES - ! - path_to_input = './' - path_to_output = './' - output_name = ' ' - interp_levels = -99999. - process = 'all' - - ! - ! READ NAMELIST - ! - DO funit=10,100 - INQUIRE(unit=funit, opened=is_used) - IF (.not. is_used) EXIT - END DO - OPEN(funit,file='namelist.api',status='old',form='formatted',iostat=ios) - IF ( ios /= 0 ) STOP "ERROR opening namelist.api" - READ(funit,io) - READ(funit,interp_in) - CLOSE(funit) + LOGICAL :: oldvar!=.FALSE. + + REAL :: onelevel + if ( onelevel .ne. -99999. ) then + interp_levels(1) = onelevel + interp_levels(2:) = -99999. + endif ! @@ -152,11 +193,11 @@ ! lent = len_trim(path_to_input) - IF ( path_to_input(lent:lent) /= "/" ) THEN - path_to_input = TRIM(path_to_input)//"/" - ENDIF - lent = len_trim(path_to_output) - IF ( path_to_output(lent:lent) /= "/" ) THEN - path_to_output = TRIM(path_to_output)//"/" - ENDIF + !IF ( path_to_input(lent:lent) /= "/" ) THEN + ! path_to_input = TRIM(path_to_input)//"/" + !ENDIF + !lent = len_trim(path_to_output) + !IF ( path_to_output(lent:lent) /= "/" ) THEN + ! path_to_output = TRIM(path_to_output)//"/" + !ENDIF input_name = TRIM(path_to_input)//TRIM(input_name) ! @@ -526,5 +567,5 @@ write(6,*) " Data will be output on unstaggered grid " do kk = 1, times_in_file - IF ( DEBUG ) print *, kk + !IF ( DEBUG ) print *, kk IF (oldvar) THEN interm1(1:iweg-1,:,:) = ( u(1:iweg-1,:,:,kk) + u(2:iweg,:,:,kk) ) * .5 @@ -673,5 +714,5 @@ ENDDO deallocate (data1) - PRINT *, pres_field(10,10,:,1) + !PRINT *, pres_field(10,10,:,1) ENDIF @@ -949,5 +990,5 @@ DO ii = 1,4 dims_out(ii) = dvalj(jshape(ii)) - print *, dims_out(ii) + !print *, dims_out(ii) ENDDO !!! NB: what follows is useful because we'd like diagnostics for each history timestep @@ -1134,5 +1175,6 @@ write(6,*) "##########################################" - END PROGRAM api +END SUBROUTINE +! END PROGRAM api !--------------------------------------------------------------------- !--------------------------------------------------------------------- @@ -1573,5 +1615,4 @@ real, parameter :: RAD_PER_DEG = PI/180. -! print *, 'map ', map_proj IF ( map_proj .ge. 3 ) THEN ! No need to rotate !PRINT *, 'NO NEED TO ROTATE !!!! equivalent to output U,V with unstagger_grid' @@ -1579,5 +1620,4 @@ VVVmet(:,:,:) = VVV ELSE - !END IF cone = 1. ! PS @@ -1606,7 +1646,4 @@ END DO -!print *, longi(10,10) -!print *, lati(10,10) - DO i = 1, west_east_dim @@ -1624,6 +1661,4 @@ UUUmet(:,:,k) = VVV(:,:,k)*sin(alpha) + UUU(:,:,k)*cos(alpha) VVVmet(:,:,k) = VVV(:,:,k)*cos(alpha) - UUU(:,:,k)*sin(alpha) -!print *,UUU(10,10,k), UUUmet(10,10,k) -!print *,VVV(10,10,k), VVVmet(10,10,k) END DO END IF Index: trunk/MESOSCALE/LMD_MM_MARS/SRC/POSTPROC/comp_api =================================================================== --- trunk/MESOSCALE/LMD_MM_MARS/SRC/POSTPROC/comp_api (revision 179) +++ trunk/MESOSCALE/LMD_MM_MARS/SRC/POSTPROC/comp_api (revision 186) @@ -5,8 +5,8 @@ # -pgf90 api.F90 -L$NETCDF/lib -lnetcdf -lm -I$NETCDF/include -Mfree -o api +#pgf90 api.F90 -L$NETCDF/lib -lnetcdf -lm -I$NETCDF/include -Mfree -o api -#g95 api.F90 -L$NETCDF/lib -lnetcdf -lm -I$NETCDF/include -Wall -Wno=112,141,137,155 -fno-second-underscore -ffree-form -o api +g95 api.F90 -L$NETCDF/lib -lnetcdf -lm -I$NETCDF/include -Wall -Wno=112,141,137,155 -fno-second-underscore -ffree-form -o api #g95 api.F90 -L$NETCDF/lib -lnetcdf -lm -I$NETCDF/include -Wall -Wno=112,141,137,155 -ffree-form -o api #pgf90 -mcmodel=medium -Mlarge_arrays api.F90 -L$NETCDF/lib -lnetcdf -lm -I$NETCDF/include -Mfree -o p_interp Index: trunk/MESOSCALE/LMD_MM_MARS/SRC/POSTPROC/namelist.api =================================================================== --- trunk/MESOSCALE/LMD_MM_MARS/SRC/POSTPROC/namelist.api (revision 179) +++ trunk/MESOSCALE/LMD_MM_MARS/SRC/POSTPROC/namelist.api (revision 186) @@ -1,63 +1,12 @@ &io -! path_to_input = './dossier/ALT_POLAR_top0.1_sponge_smooth_topo/' !! ne doit pas etre trop long -! path_to_output = './smooth' - path_to_input = './dossier/ALT_POLAR_top0.1_sponge/' + path_to_input = './' path_to_output = './' -! input_name = 'wrfout_d01_2024*' - input_name = 'wrfout_d01_2024-07-03_06:00:00' -! path_to_input = './dossier_test/' -! input_name = 'wrfout_d01_2024-01-05_02:00:00' -! input_name = 'nca_wrfout_d01_2024-01-05_02:00:00' + input_name = 'wrfout_d01_9999-09-09_09:00:00' process = 'list' !! list fields required in "fields" (available tk, tpot, GHT) -! fields = 'W,SWDOWNZ,TAU_DUST,TSURF,XLONG,XLAT,HGT' -! fields = 'W' -! fields = 'U,V' !! beware, not only tk ! only one 3D field at least fields = 'tk,W,uvmet' - debug = .TRUE. / - process = 'all' !! process all fields in file - debug = .TRUE. - bit64 = .TRUE. ! bit64 = .TRUE. !! ne pas utiliser ? &interp_in -! interp_levels = 1000.,987.5,975.,962.5,950.,937.5,925., -! 900.,875.,850.,825.,800.,750.,700.,650., -! 600.,550.,500.,450.,400.,350.,300.,250., -! 225.,200.,175.,150.,137.5,125.,112.5,100., -! 87.5,75.,62.5,50.,37.5,25.,12.5, -!interp_levels = 5., 4.9, 4.8, 4.7, 4.6, 4.5, 4.4, 4.3, 4.2, 4.1, 4.0, -! 3.9, 3.8, 3.7, 3.6, 3.5, 3.4, 3.3, 3.2, 3.1, 3.0, -! 2.9, 2.8, 2.7, 2.6, 2.5, 2.4, 2.3, 2.2, 2.1, 2.0, -! 1.9, 1.8, 1.7, 1.6, 1.5, 1.4, 1.3, 1.2, 1.1, 1.0, -! 0.9, 0.8, 0.7, 0.6, 0.5, 0.4, 0.3, 0.2, 0.1, -! interp_levels = 5., 3., 1., 0.1 -! interp_method = 2, - interp_method = 3 !! INTERPOLATION: ALTITUDE ABOVE MOLA AREOID (km) -! interp_method = 4 -! interp_levels = 20. -! interp_levels = -09., -08., -07., -06., -05., -04., -03., -02., -01., 00., -! 01., 02., 03., 04., 05., 06., 07., 08., 09., -! interp_levels = 01., 02., 03., 04., 05., 06., 07., 08., 09., -! 10., 11., 12., 13., 14., 15., 16., 17., 18., 19., -! 20., 21., 22., 23., 24., 25., 26., 27., 28., 29., -!! 30., 31., 32., 33., 34., 35., 36., 37., 38., 39., -!! 40., 41., 42., 43., 44., 45., 46., 47., 48., 49. -! interp_levels = 01., 01.5, 02., 02.5, 03., 03.5, 04., 04.5, 05., 05.5, 06., 06.5, 07., 07.5, 08., 08.5, 09., 09.5 -! 10., 10.5, 11., 11.5, 12., 12.5, 13., 13.5, 14., 14.5, 15., 15.5, 16., 16.5, 17., 17.5, 18., 18.5, 19., 19.5 -! 20., 20.5, 21., 21.5, 22., 22.5, 23., 23.5, 24., 24.5, 25., 25.5, 26., 26.5, 27., 27.5, 28., 28.5, 29., 29.5 - interp_levels = 20.0 -! interp_levels = 01.00,01.25,01.50,01.75,02.00,02.25,02.50,02.75,03.00,03.25,03.50,03.75,04.00,04.25,04.50,04.75,05.00, -! 05.25,05.50,05.75,06.00,06.25,06.50,06.75,07.00,07.25,07.50,07.75,08.00,08.25,08.50,08.75,09.00, -! 09.25,09.50,09.75,10.00,10.25,10.50,10.75,11.00,11.25,11.50,11.75,12.00,12.25,12.50,12.75,13.00, -! 13.25,13.50,13.75,14.00,14.25,14.50,14.75,15.00,15.25,15.50,15.75,16.00,16.25,16.50,16.75,17.00 -! interp_method = 2, -! interp_levels = 1., + interp_method = 4 + interp_levels = 0.050 / - extrapolate = 1, - interp_method = 2, - unstagger_grid = .TRUE. !! pb !! non c OK - - extrapolate = 0 ;; set values below ground and above model top to missing (default) - extrapolate = 1 ;; extrapolate below ground, and set above model top to model top values - interp_method = 1 ;; linear in p interpolation (default) - interp_method = 2 ;; linear in log p interpolation