Changeset 2803 for LMDZ5/trunk/libf
- Timestamp:
- Feb 28, 2017, 10:05:52 PM (8 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
LMDZ5/trunk/libf/phylmd/rrtm/read_rsun_rrtm.F90
r2742 r2803 10 10 USE netcdf, ONLY: nf90_get_var, nf90_noerr, nf90_nowrite 11 11 12 USE phys_cal_mod, ONLY : day _cur, year_len12 USE phys_cal_mod, ONLY : days_elapsed, year_len 13 13 14 14 USE mod_phys_lmdz_mpi_data, ONLY: is_mpi_root … … 23 23 24 24 ! Input arguments 25 LOGICAL, INTENT(IN) 25 LOGICAL, INTENT(IN) :: debut 26 26 27 27 ! Local variables 28 INTEGER 28 INTEGER :: ncid, dimid, varid, ncerr, nbday 29 29 REAL, POINTER :: wlen(:), time(:) 30 REAL, ALLOCATABLE, SAVE, DIMENSION(:,:) 30 REAL, ALLOCATABLE, SAVE, DIMENSION(:,:) :: SSI_FRAC 31 31 !$OMP THREADPRIVATE(SSI_FRAC) 32 32 REAL, ALLOCATABLE, SAVE, DIMENSION(:) :: TSI(:) … … 90 90 ENDIF 91 91 92 !--only read at beginning of month 93 IF (debut.OR.day_cur.NE.day_pre) THEN 92 !--only read at beginning of day 93 !--day in year is provided as days_elapsed since the beginning of the year +1 94 IF (debut.OR.days_elapsed+1.NE.day_pre) THEN 94 95 95 !--keep memory of previous month96 day_pre=day _cur96 !--keep memory of previous day 97 day_pre=days_elapsed+1 97 98 98 99 !--copy 99 RSUN(1:NSW)=SSI_FRAC(:,day _cur)100 solaire=TSI(day _cur)100 RSUN(1:NSW)=SSI_FRAC(:,days_elapsed+1) 101 solaire=TSI(days_elapsed+1) 101 102 102 print *,'READ_RSUN_RRTM day=', day _cur,' solaire=', solaire, ' RSUN=', RSUN(1:NSW)103 print *,'READ_RSUN_RRTM day=', days_elapsed+1,' solaire=', solaire, ' RSUN=', RSUN(1:NSW) 103 104 104 105 ENDIF !--fin allocation
Note: See TracChangeset
for help on using the changeset viewer.