| 102 | * '''thermosoil.f90 / thermosoil_finalize''' (2015/08/31) |
| 103 | The gtemp variable was not saved in the restart file.[[BR]] |
| 104 | The following correction should be applied: |
| 105 | {{{ |
| 106 | SUBROUTINE thermosoil_finalize(kjit, kjpindex, rest_id, temp_sol_new, gtemp, soilcap, soilflx) !@FM_150831 |
| 107 | !! 0. Variable and parameter declaration |
| 108 | !! 0.1 Input variables |
| 109 | INTEGER(i_std), INTENT(in) :: kjit !! Time step number (unitless) |
| 110 | INTEGER(i_std), INTENT(in) :: kjpindex !! Domain size (unitless) |
| 111 | INTEGER(i_std),INTENT (in) :: rest_id !! Restart file identifier(unitless) |
| 112 | REAL(r_std),DIMENSION (kjpindex), INTENT (in) :: temp_sol_new !! Surface temperature at the present time-step, |
| 113 | REAL(r_std),DIMENSION (kjpindex),INTENT(in) :: gtemp !! First soil layer temperature !@FM_150831 |
| 114 | ... |
| 115 | var_name= 'temp_sol_beg' |
| 116 | CALL restput_p(rest_id, var_name, nbp_glo, 1, 1, kjit, temp_sol_beg, 'scatter', nbp_glo, index_g) |
| 117 | |
| 118 | CALL restput_p(rest_id, 'gtemp', nbp_glo, 1, 1, kjit, gtemp, 'scatter', nbp_glo, index_g) !@FM_150831 |
| 119 | ... |
| 120 | }}} |
| 121 | |
| 122 | with the corresponding call in sechiba_finalize: |
| 123 | |
| 124 | {{{ |
| 125 | !! 5. Call soil thermodynamic to write restart files |
| 126 | !@FM_150831 CALL thermosoil_finalize(kjit, kjpindex, rest_id, temp_sol_new, soilcap, soilflx) |
| 127 | CALL thermosoil_finalize(kjit, kjpindex, rest_id, temp_sol_new, gtemp, soilcap, soilflx) !@FM_150831 |
| 128 | }}} |
| 129 | |
| 130 | The MICT thermosoil_finalize is different from the Trunk version so probably other variables have to be checked. |